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ABSTRACT 


A  steady  flow  analysis  is  conducted  for  a  NACA  0012  airfoil  in  low  Reynolds 
number  flows  ranging  from  540,000  to  1,000,000.  Emphasis  is  placed  on  prediction  and 
location  of  the  separation  bubble.  Computational  methods  include  the  direct  boundary 
layer  method,  the  viscous-inviscid  interaction  method,  and  the  time-averaged  Navier- 
Stokes  method.  Characteristic  trends  in  skin  friction  coefficient,  displacement  thickness, 
and  boundary  layer  velocity  profiles  with  increasing  angle  of  attack  are  observed. 
Computational  results  are  compared  to  each  other  and  to  experimental  photographs 
visualizing  the  density  flowfield  using  Point  Diffraction  Interferometry.  Both  the 
viscous-inviscid  method  and  the  Navier-Stokes  method  failed  to  accurately  represent 
leading  edge  separation  bubbles.  The  direct  boundary  layer  method,  usually  considered 
of  very  limited  usefulness  due  to  a  singularity  in  the  underlying  equations  at  separation, 
is  shown  to  exhibit  unexpected  recovery  behavior  for  small  amounts  of  separation. 
Furthermore,  the  results  near  the  leading  edge,  where  separation  bubbles  were  computed, 
were  validated  by  the  experiment. 
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1.  INTRODUCTION 


While  there  are  many  reliable  solution  methods  for  high  Reynolds  number  flows, 
the  low  Reynolds  number  regime  currently  has  considerably  fewer  options.  This  is  due 
in  part  to  the  fact  that  high  Reynolds  number  flows  account  for  most  aeronautical 
applications  of  interest.  There  are,  however,  important  applications  involving  low 
Reynolds  numbers,  such  as  turbomachinery  blades.  The  more  likely  reason  for  the  lack 
of  reliable  codes  for  low  Reynolds  numbers  is  the  greater  difficulty  of  accurately 
representing  the  flow.  Most  methods  make  use  of  approximations  in  the  formulation  of 
their  underlying  equations  to  obtain  computational  solutions  in  a  reasonable  amount  of 
time.  These  approximations  often  become  less  and  less  accurate  as  the  Reynolds  number 
decreases.  An  obvious  question  which  may  arise  concerns  the  value  of  bothering  with 
such  seemingly  limited,  simplistic  codes  when  the  state  of  the  art  is  Navier-Stokes  (NS) 
solvers.  There  are  several  considerations  which  make  the  effort  worthwhile.  First,  not 
everyone  has  access  to  the  supercomputers  or  mini-supercomputers  that  are  necessary  for 
extensive  NS  solutions.  Even  if  these  computers  are  available,  user  time  may  be  limited. 
Preliminary  research  using  a  simpler  method  may  cut  down  the  amount  of  advanced 
calculations  needed  considerably,  thus  reducing  overall  cost.  Another  pertinent  factor 
is  the  time  savings.  A  design  team  for  a  new  aircraft  cannot  afford  months  of  detailed 
refinement  using  NS  solvers  when  a  proposal  deadline  looms  near.  Less  than  NS 
accuracy  is  certainly  acceptable,  especially  when  the  significant  gains  in  speed  and  cost 
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reduction  are  considered.  An  efficient  program  only  needs  to  calculate  to  the  level  of 
accuracy  necessary  to  accomplish  the  desired  goal.  To  this  end,  methods  such  as  the 
direct  boundary  layer  method  and  the  viscous-inviscid  interaction  method  offer 
opportunities  to  investigate  low  Reynolds  number  phenomena,  such  as  laminar  separation 
bubbles.  A  final  point  to  consider  is  that  all  methods,  including  NS  solvers,  are  really 
approximations  in  that  they  involve  empirical  models  for  transition  and  turbulence.  The 
suitability  of  these  models  to  low  Reynolds  number  flows  will  influence  the  results 
directly. 

This  investigation  explores  the  various  computational  methods,  comparing  them  to 
each  other  and  to  experimental  data.  A  panel  code  was  first  developed,  which  provided 
some  of  the  input  needed  for  a  direct  boundary  layer  code.  The  direct  boundary  layer 
code  was  studied  extensively  and  several  modifications  were  made  to  enable  further 
analysis  of  boundary  layer  profiles  and  transition  effects.  Similar  calculations  were  then 
performed  for  a  viscous-inviscid  interaction  code.  Experimental  interferograms,  obtained 
by  Point  Diffraction  Interferometry,  provided  a  reliable  reference  for  comparison. 
Finally,  a  Navier-Stokes  code  was  evaluated.  Emphasis  throughout  the  investigation  was 
placed  on  detection  and  location  of  laminar  separation  bubbles,  as  well  as  a  thorough 
consideration  of  transition  and  turbulence  models. 
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n  PANEL  CODES 


A.  THEORY 

In  potential  flow  theory,  the  flow  field  around  an  airfoil  may  be  represented  by  the 
velocity  potential.  Considering  contributions  from  the  freestream  flow  and  the  source 
and  vorticity  distributions,  the  total  potential  may  be  constructed: 


*=  4),+  4>v 


where 


(2.1) 


cosa  +y  sina) 


(2.2) 


The  source  distributions  (q)  vary  from  panel  to  panel,  while  the  vorticity  strength  (7)  is 
assumed  constant  for  all  panels.  The  value  of  representing  the  flow  past  an  airfoil  by 
surface  singularity  distributions  lies  in  the  fact  that  these  singularity  distributions 
automatically  satisfy  Laplace’s  equation,  the  governing  flow  equation  for  inviscid, 
incompressible  flow: 


<l>jB  +  ^jy=0  (2.3) 

Since  Laplace’s  equation  is  a  linear  homogeneous  second  order  partial  differential 
equation,  the  superposition  principle  used  in  Equation  2.1  is  valid.  The  boundary 
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conditions  include  flow  tangency  at  control  points  (midpoints  of  panels)  and  the  Kutta 
condition  at  the  trailing  edge,  requiring  equal  tangential  velocities  for  the  first  and  last 
panels.  By  evaluating  the  integrals  along  the  airfoil  surface,  the  potential  may  be 
determined  at  any  point  in  the  flow  field.  Each  point  is  defined  at  a  radius  (r)  and  angle 
(0)  from  a  chosen  reference  point  on  the  airfoil.  The  reference  point  in  this  study  is  the 
leading  edge. 

B.  PANEL  METHOD  GEOMETRY 

For  computational  purposes,  it  is  not  feasible  to  evaluate  every  point  in  the  flow 
field.  The  airfoil  is  represented  by  a  number  of  defined  points,  called  nodes.  More 


connecting  these  nodes  are 

the  panels.  There  are  (n)  panels  and  (n+l)  nodes,  with  the  first  and  last  node 
overlapping.  Figure  Figure  2.1  depicts  the  panel  geometry.  Numbering  starts  at  the 
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trailing  edge,  then  progresses  along  the  lower  surface,  leading  edge,  and  upper  surface, 
and  ends  at  the  trailing  edge.  The  unit  normal  vectors  (fi^)  are  perpendicular  to  the 
panels  and  directed  outward  from  the  airfoil  surface.  The  unit  tangential  vectors  (t,)  are 
parallel  to  the  panels  and  the  positive  direction  is  defined  with  increasing  numbering  (n 
to  n+ 1).  The  panels  may  vary  in  length,  with  the  exception  of  the  first  and  last  panels, 
which  must  be  equal  in  order  to  use  the  Kutta  condition  at  the  trailing  edge. 

C.  COMPUTER  CODE 

1.  Overview 

In  order  to  study  the  steady,  incompressible,  inviscid  flow  over  arbitrary 
airfoils,  a  panel  code  called  panel  was  developed.  The  required  input  consists  of  the 
number  of  nodes  on  the  airfoil  surface,  the  coordinates  of  the  nodes  referenced  from  the 
leading  edge,  and  the  angle  of  attack  in  degrees.  The  program  produces  normalized 
velocities  and  pressure  coefficients  at  each  control  point  as  output.  The  program  was 
later  modified  to  produce  an  output  file  compatible  with  bl2d,  a  direct  boundary  layer 
program  described  in  Chapter  2.  Additional  input  consists  of  Reynolds  number  and 
transition  information.  This  data  is  not  used  by  the  program  panel,  but  is  simply 
transferred  to  the  output  file  which  will  be  used  as  input  for  bl2d. 

2.  Influence  Coefficients 

The  use  of  influence  coefficients  leads  to  a  straightforward  procedure  for 
programming  the  equations.  An  influence  coefficient  is  defined  as  the  velocity  induced 
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at  a  field  point  by  a  unit  strength  singularity  distribution  on  one  panel.  For  the  two 
dimensional  steady  flow  problem,  the  following  influence  coefficients  are  needed: 


•  A“j;  normal  velocity  component  induced  at  the  i*  panel  control  point  by  unit 
source  distribution  on  the  j*  panel 


i4  =  .X[sm(0.-0^.)ln^  +cos(0j-0^)  i*j 

1/ 

1  .  . 

-j  .  -j 


(2.4) 


•  A*ij:  tangential  velocity  component  induced  at  the  i*  panel  control  point  by  unit 
strength  source  distribution  on  the  j*  panel 


A 


\  =  ^[sin(0,.-0j.) pj,.  -cos(0,.-0^.)ln^],  i*j 
=0  ,  i=j 


(2.5) 


•  B“ij:  normal  velocity  component  induced  at  the  i*  panel  control  point  by  unit 
strength  vorticity  distribution  on  the  j*  panel 


-sin(0re,)Pj/].  i*J 

=0  ,  i=j 


(2.6) 


•  B'jj:  tangential  velocity  component  induced  at  the  i'*’  panel  control  point  by  unit 
strength  vorticity  distribution  on  the  j*  panel 


B  = ±Ccos(e,-ej)  p, .  since, 

r^j 

~  2 


],  i*j 

.  H 


(2.7) 


where  the  geometrical  quantities,  depicted  in  Figure  2.2,  are  defined  by: 
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(2.8) 


ry  =  <Jixmi-x/  *(ymry/ 


xm,^ - 

'  2 


ymr 


0j  =arctan(^'**-^') 


J3n,-jcy„ 


xnif-Xj 


i  +  1  (xnii  ,  yiTii) 


i"'  panel 

Figure  2.2  Relationships  Betu'een  Geometrical  Quantities 
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The  formula  for  /Sy  may  be  verified  as  follows: 

a  +  +  X  +  /3y  =  180“  (triangle) 

X  +  b  +  =  180“  (supplementary  angles) 

Setting  these  equations  equal  to  each  other  and  eliminating  common  terms, 

0,  =  b-a 

Inspection  of  the  diagram  shows  that  angle  b  is  in  fact  the  arctangent  of  the  quantity  in 
parentheses  in  the  first  term  of  the  formula  for  0^.  Likewise,  angle  a  matches  the  second 
term. 

3.  Program  Description 

a.  Boundary  conditions 

The  first  boundary  condition  requires  flow  tangency  at  control  points: 

(K"),=0  ,  i  =  l,2,...n  a-9) 

In  terms  of  influence  coefficients  (with  V*  =  l), 


Z;.,  .  '->.2 . " 

The  second  boundary  condition  is  the  Kutta  condition,  which  states  that 
the  pressures  on  the  lower  and  upper  panels  at  the  trailing  edge  must  be  equal  if  the  flow 
is  to  leave  the  trailing  edge  smoothly.  Using  a  form  of  Bernoulli’s  equation. 
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(2.11) 


the  pressure  equilibrium  also  implies  equal  velocities  for  incompressible  flow.  Since  the 
normal  velocities  are  taken  to  be  zero,  the  boundary  condition  may  now  be  stated  as: 

(V%-~iVX  (2.12) 

where  the  negative  sign  is  strictly  due  to  the  adopted  convention  of  positive  tangential 
velocities  in  the  direction  of  increasing  node  numbering.  Since  the  flow  is  positive  to 
the  right  (as  shown  in  Figure  2.1),  the  panels  downstream  of  the  front  stagnation  point 
will  have  negative  values  for  computational  purposes  only.  It  is  important  to  note  that 
not  all  the  lower  surface  panels  have  a  reversed  sign,  only  those  downstream  from  the 
stagnation  point.  This  is  especially  significant  for  non-symmetrical  airfoils  or  any  airfoil 
at  an  angle  of  attack. 

In  terms  of  influence  coefficients,  the  normalized  equation  becomes: 


-5^.. 

r;.,  S'.)*"*'*-®.) 


h.  Solution  procedure 

Equations  2.10  and  2.13  represent  a  linear  algebraic  system  of  (n+1) 
equations  and  (n+1)  unknowns.  The  unknowns  are  the  source  strengths  which  vary  from 
panel  to  panel  (qi...qj  and  the  vorticity  strength  y. 
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Expanding  and  rearranging  Equation  2.10  for  an  example  airfoil  of 


n=73  nodes  and  panels  results  in: 


♦...+  B",73)  =  -sin(a-0,) 

^"21^1  ^"22^2 "21  +^"22  B\73)  =  -sm(a-02)  (2.14) 

"73,1 "73^92  Y  "73,1  ®  "73^  + ...  +  B  "73^73)  =  -  sin(a  -073) 


The  equations  now  readily  lend  themselves  to  solution  in  matrix  form.  Recasting  with 
a  simpler  notation,  the  A“ij  terms  (coefficients  of  q^)  may  be  renamed  a^j  and  the  sum  of 
all  B“j  terms  in  parentheses  (coefficients  of  7)  renamed  ain+,,  where  i=l,2,...,n  and 
j  =  l,2,...,n.  The  terms  on  the  right  sides  of  the  equations  may  be  renamed  bj. 

The  (n+l)**  equation,  or  in  this  example,  the  74*  equation,  comes  from 
Equation  2.13  in  a  similar  manner: 


...  *(B',,„*B'„,„)].  a.i» 

cos(a-0i)-cos(a-073) 

The  coefficients  of  qj  may  be  renamed  a74j.  All  of  the  B*  terms  in  the  brackets  together 
form  the  coefficient  of  7,  now  renamed  a74  74.  The  entire  right  side  of  the  equation 
constitutes  the  new  term  b74. 
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Finally  expressing  this  system  in  a  concise  matrix  form  for  the  general 


case, 


.  ^1^*1 

^1  ^  ^  .  ^*1  ^2  ^2 

^.1  ^,1  ^33  . — •  ^3  ^  h  (2.16) 

^11,1  ®/i;2  ^<1,3  .  9h 

-  Y. 


This  system  is  solved  in  the  program  using  a  Gaussian  Elimination  subroutine. 

With  the  values  of  the  and  y  known,  the  velocity  at  each  panel 
control  point  may  be  calculated: 

»  =  (2.17) 

The  total  velocity  is  equal  to  the  tangential  velocity  due  to  taking  the  normal  velocity  to 
be  zero. 

c.  Numerical  techniques 

Although  programming  most  of  the  described  procedure  is  a  relatively 
straightforward  task,  there  are  a  number  of  potential  pitfalls.  Different  programming 
languages  each  have  their  own  special  rules  and  format,  but  the  following  hints  for 
FORTRAN,  a  commonly  used  engineering  language,  apply  to  many  others  as  well. 

All  angles  entered  into  the  equations,  such  as  a,  /S,  and  6,  must  have 
values  in  radians.  The  angle  of  attack  (a)  is  an  input  parameter  that  should  also  be 
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converted  relative  to  the  zero  lift  line  of  the  airfoil  in  order  for  the  resulting  velocities 
to  match  published  values.  The  angle  6  is  used  in  calculating  the  influence  coefficients. 
6  is  the  angle  of  a  panel  from  the  positive  x  axis,  counter-clockwise  positive.  The 
formula  given  for  d,  involving  taking  an  arctangent,  will  produce  the  correct  physical 
angle  if  interpreted  correctly.  Many  programming  languages  use  -t/2  to  t/2  as  the 
default  range  for  the  standard  inverse  tangent  function,  which  uses  only  one  value  for  an 
argument.  The  function  cannot  determine  whether  a  negative  sign  was  in  the  numerator 
or  denominator.  A  problem  arises  when  an  angle  is  actually  in  the  2“*  or  3"*  quadrant 
because  the  function  will  assign  values  from  the  1“  or  4*  quadrant.  The  values  for  6 
must  be  in  the  range  -r<d<x  to  work  properly.  This  may  be  accomplished  by  using 
the  ATAN2  function  in  FORTRAN,  which  accepts  both  a  numerator  and  denominator 
as  arguments  and  assigns  quadrants  correctly. 

The  angle  jS  may  be  calculated  from  two  inverse  tangents,  as  presented 
in  the  formula.  However,  a  more  efficient  algorithm  may  be  used  for  computer 
calculation.  Recalling  that  0  =  b-a,  it  follows  that: 
tan  /3  =  tan(b-a) 

=  Sigftba) 

cos(b-a) 

=  sin  b  cos  a  -  cos  b  sin  a 

cos  b  cos  a  +  sin  b  sin  a  (2.18) 

The  sines  and  cosines  for  the  angles  a  and  b  may  be  easily  determined  from  the  geometry 
of  Figure  2.2.  For  example,  sin  b  =  (yntj  -  yj+,)/hypotenuse  b.  Noting  that  all  the 
hypotenuse  values  may  be  cancelled  out  in  Equation  2.18,  the  arctangent  of  the  angle 
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may  be  found  with  simply  x  and  y  differences.  The  ATAN2  function  in  FORTRAN 
should  also  be  used  to  calculate  /3. 

D.  RESULTS 

1.  Eppler  E585  AirfoU 

The  first  investigation  was  conducted  for  an  Eppler  E585  airfoil  (shown  in 
Figure  2.3),  with  n=71  nodes  (or  panels).  It  is  one  of  a  series  of  airfoils  without  flaps 
designed  for  the  Reynolds  number  range  of  sailplanes,  about  100,000  to  500,000.  The 
angle  relative  to  the  zero  lift  line  is 
5.53°.  If  not  compensated  for,  the  results 
would  all  be  shifted  by  this  amount. 

Figure  2.4  shows  the  results  of  the 
program  panel  for  a  =  3,  7,  11°.  The 
published  data  in  Figure  2.5  [Ref.  IJ 
compares  favorably.  The  velocities  match 
the  given  distribution  well  except  for 
slight  deviations  at  the  trailing  edge.  This 
difference  can  be  attributed  to  the  higher 
order  panel  method  used  in  Eppler’s  computations.  The  panels  are  defined  by  third 
degree  polynomials  whose  coefficients  are  computed  by  a  spline  program.  Additional 
points  are  splined  in  as  needed.  For  the  best  precision  from  a  panel  method,  the  steps 
must  be  smaller  near  the  leading  and  trailing  edges.  The  simpler  method  of  connecting 
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Figure  2.4  Velocity  Distribution  Computed  by  PANEL 


just  the  given  airfoil  coordinates  naturally  leads  to  less  accurate  results  in  the  critical 
trailing  edge  area.  Experimentation  with  linearly  interpolated  extra  points  verified  that 
a  better  resolution  could  be  obtained  in  this  manner.  For  most  cases,  the  straight  line 
segment  panels  using  only  the  given  coordinates  produce  quite  sufficient  resolution  with 
a  great  advantage  in  computational  speed. 

2.  NACA  0012  Airfoil 

This  is  a  well  known  airfoil  (shown  in  Figure  2.6)  used  for  many  studies  and 
comparisons,  elementary  to  advanced,  as  well  as  practical  applications.  It  is  one  of  the 
original  4-digit  series  of  1932, 
where  the  first  two  digits  indicate 
camber  amount  and  location  (00  is 
symmetrical)  and  the  last  two 
digits  indicate  maximum  thickness 
in  percent  mean  aerodynamic 
chord  (12%).  Figures  2.7  and  2.8 
show  a  comparison  of  program 
results  and  those  provided  by  Anderson  [Ref.  2]  for  the  NACA  0012  airfoil  at  9°  angle 
of  attack.  The  pressure  distributions  both  come  to  a  suction  peak  of  -5.2.  It  should  be 
noted  that  Figure  2.8  differs  slightly  from  the  original  reference  plot  in  that  the  error  in 
the  decimal  place  of  the  ordinate  values  has  been  corrected. 
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Figure  2.8  Pressure  CoefTicient  Provided  by  Anderson 
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E.  USER’S  GUIDE 


These  detailed  instructions  are  accurate  at  the  time  of  writing,  oriented  toward  the 
Advanced  Computation  Laboratory  (rm.  136)  of  the  Naval  Postgraduate  School. 
However,  due  to  the  dynamic  nature  of  ariy  computer  environment,  some  instructions 
may  change  over  a  period  of  time.  If  difficulty  is  encountered,  check  with  the  manager 
for  changes  in: 

1 .  the  account  location  of  the  files 

2.  plotting  package  availability  or  procedures 

3.  computer  informal  addresses  .i.e.  madmax,  suzqt,  indigol,  etc.) 

In  all  instruction  lines,  the  prompt  which  appears  on  the  screen  is  in  normal  type,  while 
literal  user  input  is  in  boldface.  User  input  which  is  general  and  requires  the  actual 
word  instead  will  be  in  italics.  For  example,  filename  could  be  replaced  by  vel.dat. 
Although  not  specifically  stated  on  each  line,  a  carriage  return  (ENTER)  is  assumed. 
As  many  commands  are  case  sensitive,  type  each  line  in  the  format  shown. 

1.  Stardent 

a.  Program  operation 

Using  the  Stardent  terminal,  login  and  change  to  the  directory  to  be 
used  for  the  panel  code,  making  a  new  one  if  desired  (i.e.  mkdir  paneldir,  cd 
paneldir).  Then  type; 

>cp  /alpha/acc/mune/panel  panel 
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to  copy  the  compiled  program  from  the  account  where  it  is  stored.  Also  copy  a  sample 
input  file  for  a  NACA  0012  airfoil: 

>  cp  /alpha/accrname/points.dat  points.dat 

These  two  files  are  all  that  is  needed  to  begin.  To  run  the  program,  type; 

>  panel 

Answer  the  prompted  questions,  using  100  points  and  0  for  the  transition  code.  When 
finished,  the  program  will  respond  with: 

Calculations  complete,  output  files  are: 
vel.dat,  cp.dat,  bl2d.dat 

The  first  two  files  are  simply  x-y  type  column  data  for  plotting  the  velocity  and  pressure 
distributions.  The  last  file  is  created  for  use  with  another  program,  bl2d,  described  in 
the  next  chapter. 

b.  Plotting  procedures 

To  view  the  resulting  velocity  distribution,  the  file  vel.dat  can  be  used 
wiw.  any  standard  plotting  package.  On  the  Stardent,  gnuplot  can  be  used.  First  type; 

>  xterm  -t  & 

to  open  up  a  tektronics  window  for  plotting  that  will  operate  in  background.  When  it 
appears,  move  the  cursor  there  and  type: 

>  gnuplot 

Note  that  it  is  somewhat  "messy"  to  work  in  this  window  because  mistakes  do  not 
disappear  from  the  screen  with  a  backspace.  The  computer  will  take  the  overwritten 
characters  as  the  input,  however.  Control-C  may  be  pressed  instead  to  simply  type  the 
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line  over  or  type  clear  to  erase  the  whole  screen  when  needed.  At  the  gnuplot  prompt, 
type 

gnuplot  >  set  term  tek40xx 

to  set  the  proper  terminal  type  for  plotting.  To  plot  the  x-y  format  data  file  vel.dat, 
type: 

gnuplot  >  plot  "vel.dat" 

For  options  in  gnuplot  such  as  scaling  axes  and  plotting  multiple  data  sets,  see  Appendix 
B. 

2.  IRIS 

a.  Using  a  Stardent  window 

To  simply  operate  the  program  panel  from  an  IRIS  terminal,  login  and 
open  up  a  Stardent  window  using  the  following  procedures.  First,  click  the  right  mouse 
button  inside  the  original  console  window  and  select  Clone.  Move  the  cursor  to  the  new 
window  and  change  to  the  desired  directory  on  the  IRIS,  making  a  new  one  (mkdir)  if 
desired.  Then  type: 

>  telnet  suzqt 

Login  to  the  Stardent  and  follow  the  instructions  in  the  previous  section  to  operate  the 
program.  Do  not  use  the  same  plotting  procedures  when  complete,  however.  Keeping 
both  windows  open,  move  the  cursor  back  to  the  IRIS  window  and  type: 

>  rep  sazqU/alph^/ loginnamel directory hcLdsA  vel.dat 
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This  command  remote  copies  from  one  system  to  the  other.  The  loginname  on  the 
command  line  is  the  name  of  the  account,  usually  the  user’s  last  name.  The  directory 
is  the  one  created  on  the  Stardent  where  the  program  was  run.  This  command  will  only 
work  if  an  account  is  held  on  both  the  Stardent  and  IRIS  under  the  same  loginname.  If 
this  is  not  the  case,  files  can  be  transferred  using  the  file  transfer  utility  ftp,  described 
in  Appendix  B. 

b.  Using  the  IRIS 

It  is  also  possible  to  do  all  calculations  directly  on  the  IRIS.  This  may 
be  more  useful  when  a  user  holds  an  account  on  the  IRIS  but  not  on  the  Stardent.  After 
obtaining  a  copy  of  the  source  code  panel.f  from  the  Stardent  using  the  ftp  procedures 
(Appendix  B),  compile  it  for  the  IRIS: 

>  f77  -03  -o  panel  panel.f 

Program  operation  is  as  described  for  the  Stardent. 

c.  Plotting  procedures 

The  plotting  package  usually  used  on  the  IRIS  is  XYPLOT.  At  the 
prompt,  from  the  directory  with  the  plotting  data  (such  as  vel.dat)  type: 

>  xyplot 

Answer  the  questions  that  follow: 

Name  of  1st  input  file?  vel.dat 

Name  of  2nd  input  file?  (press  ENTER  since  only  1  file  to  plot) 

Default  configuration  file?  (ENTER,  since  none  specified  yet) 
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A  new  window  will  pop  up  with  the  plot.  The  mouse  can  be  used  in  this  window  to 
change  the  appearance  of  the  plot  in  many  ways.  The  plotting  program  is  very  user- 
friendly  and  can  be  operated  with  little  prior  instruction. 

Gnuplot  is  also  available  on  the  NFS  IRIS.  It  is  not  necessary  to  open 
a  special  window  for  the  plot  because  one  will  be  created  automatically  when  the 
program  is  invoked. 
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ra.  DIRECT  BOUNDARY  LAYER  CODE 


A.  THEORY 

The  two-dimensional  flow  around  an  airfoil  may  be  represented  by  the  simultaneous 
solution  of  the  continuity  equation, 

(3.1) 

dt  dx  dy 

the  components  of  the  momentum  equation, 

dpu  ^  dpu^  ^  dpuv  __dp  ^  ^  0^2) 

dt  dx  dy  dx  dx  dy  ^ 

dpv  ^dpuv ^  dpv\_  dp ^  >  Oj) 

dt  dx  dy  coc  dx  dy  ^ 

and  the  energy  equation,  often  collectively  referred  to  as  the  Navier-Stokes  equations, 
although  technically  this  name  applies  only  to  the  momentum  equation  applied  to  a 
Newtonian  fluid.  Equations  3. 1-3.3  are  expressed  in  general  form  for  unsteady, 
compressible,  viscous  flows  with  body  forces.  In  practice,  however,  such  a  complete 
solution  is  not  usually  feasible,  or  even  possible,  for  many  applications.  The 
development  of  the  Thin  Shear  Layer  (TSL)  equations,  sometimes  referred  to  as  the 
boundary  layer  equations,  enables  a  computationally  practical  scheme  for  solving  the 
flow. 
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Considering  the  steady,  incompressible,  flow  around  an  airfoil  with  no  body  forces, 
some  initial  simplifications  may  be  made.  The  energy  equation  becomes  decoupled  from 
the  other  equations  and  is  no  longer  needed  in  the  solution.  The  continuity  equation 
immediately  reduces  to 

(3.4) 

dx  dy 

Eliminating  the  time  dependent  term,  expanding  the  partial  derivatives,  subtracting  u 
times  the  continuity  equation,  and  dividing  through  by  p  in  the  x  momentum  equation 
yields 

=  (3.5) 

dx  dy  p  dx  p  dx  pc^ 

with  a  similar  expression  for  the  y  component.  The  equation  is  further  developed  by 
using  the  assumption  of  a  Newtonian  fluid,  in  which  stress  is  proportional  to  rate  of 
strain; 


o„=2p(-^) 

dc 


,du  dv. 


(3.6) 


(3.7) 


where  p  is  the  viscosity.  In  a  constant-property  flow,  p  may  be  taken  outside  the 
derivative  after  substituting  Equations  3.6  and  3.7  into  Equation  3.5,  and  may  be 
rewritten  in  terms  of  the  kinematic  viscosity,  vsp/p.  Thus,  the  x-component  momentum 
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equation  for  a  Newtonian  fluid  in  constant-property,  steady,  two-dimensional,  laminar 
flow  is 


du  du  I  dp 
dOc  dy  p  dx 


dx^  dy^ 


(3.8) 


To  incorporate  the  effects  of  a  turbulent  flow,  all  instantaneous  flow  quantities  are 
replaced  by  a  mean  term  plus  a  fluctuating  part,  i.e.  u  =  u  +  u'.  Expanding,  eliminating 
zero  products,  and  rearranging  the  equation  gives 

=  -1^  +v(-^  (3.9) 

dx  dy  p  dx  dx^  dy^  dx  dy 

where  the  overbars  on  the  mean  velocity  components  u  and  v  are  omitted  for  simplicity. 
The  extra  turbulent  stress  terms  in  Equation  3.9,  as  compared  to  Equation  3.8,  are  often 
called  the  Reynolds  stresses. 

To  obtain  the  TSL  equations,  an  order  of  magnitude  analysis  is  applied  with  the 
assumption  of  5/1  <  1 .  In  other  words,  the  boundary  layer  (of  thickness  6)  is  very  small 
compared  to  the  characteristic  length  of  the  body.  Using  the  following  order  of 
magnitude  approximations 


K~« 


€ 


dy  5 


du 

dx"  I 


(3.10) 


where  the  subscript  e  refers  to  the  edge  of  the  boundary  layer,  the  first  term  in  the 
parentheses  and  the  first  Reynolds  stress  in  Equation  3.9  may  be  neglected.  When 
compared  term  by  term  to  the  x-component  equation  with  the  assumption  u>v,  the  y- 
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component  of  the  momentum  equation  reduces  to  the  approximation  of  constant  pressure 
in  the  normal  direction.  Summarizing,  the  two-dimensional,  incompressible,  steady 
boundary-layer  equations  for  both  laminar  and  turbulent  flows  are: 


a*  dy 


(3.11) 


du  du  I  dp  du'v’  ,, 

dx  dy  g  dx  dy'^  dy 

(3.13) 

dy 

Note  that  these  equations  are  for  a  surface  coincident  with  the  x-axis.  In  order  to  use 
the  equations  for  an  airfoil,  the  airfoil  surface  must  be  "unwrapped"  onto  the  x-axis.  The 
usual  x/c  and  y/c  coordinates  which  define  the  airfoil  must  be  transformed  to  a  surface 
coordinate. 

The  applicable  boundary  conditions  on  the  surface  of  a  solid  airfoil  are 

y=0  u=0  v=0 

and  at  the  outer  edge  of  the  boundary  layer, 

y=5  u=Ue(x) 
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B.  COMPUTER  CODE:  BL2D 


1.  Overview 

The  program  bl2d,  developed  by  T.  Cebeci,  provides  a  solution  to  the 
boundary  layer  equations.  The  same  restrictions  apply  to  the  program  as  to  the  equations 
on  which  it  is  based:  two-dimensional,  steady,  incompressible,  viscous  flow.  The 
program  accepts  input  of  Reynolds  number  and  prescribed  transition  locations,  as  well 
as  panel  coordinate  and  velocity  information  computed  by  a  separate  routine.  Output  is 
generated  for  many  features  of  the  resulting  boundary  layer,  including  skin  friction 
coefficient  and  displacement  thickness.  Run  time  is  less  than  one  minute  on  a  Stardent 
computer  and  less  than  two  minutes  on  a  personal  computer  (PC). 

2.  Models 

a.  Turbulence  Model 

In  order  to  use  equation  3.12,  an  expression  must  be  found  for  the 
Reynolds  shear  stress  term.  Since  it  is  not  feasible  to  attempt  calculating  the  actual 
value,  empirical  models  are  usually  used.  One  such  model  is  the  eddy-viscosity  concept: 


(3.14) 


where  is  an  empirical  term  called  the  turbulent  eddy  viscosity.  Another  model  is  the 
mixing-length  concept  first  proposed  by  Prandtl  in  1925: 
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■p«V=p/^^ 


diM 

du 

dy 

dy 

(3.15) 


where  1  is  the  empirically  determined  mixing  length.  Although  both  of  these  models 
were  originally  derived  based  on  erroneous  physical  arguments,  they  have  nevertheless 
produced  remarkably  successful  results  for  many  applications.  A  third  model,  which 
incorporates  the  eddy-viscosity  concept,  is  the  Cebeci-Smith  (CS)  model,  in  which  the 
viscous  region  is  divided  into  an  inner  layer  and  an  outer  layer,  each  with  its  own 
formula  for 


where 


=  0.16^[l-e 


-(4) ,  , 
1  ^ 


0.0168 


(3.16) 


(3.17) 


u 

X 

1  -exp[-G(jc-ar^)  J 

1200(«.j  ^  'ir 

and  the  Falkner-Skan  variables  ij,  and  f  are  used.  The  term  is  a  factor  which 
models  the  length  of  the  transition  region,  explained  further  in  the  next  section.  This 
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turbulence  model  is  used  in  the  program.  Using  Equation  3.14  allows  rewriting  the 
momentum  equation  for  a  turbulent  flow,  Equation  3.12,  in  the  same  form  as  a  laminar 
flow: 

u^+v— (b— )  (3.18) 

dx  ^  p  dx  dy  dy 

where  b =»>+€„.  Thus,  the  only  computational  difference  between  laminar  and  turbulent 
boundary  layers  is  the  addition  of  the  turbulent  eddy  viscosity. 

b.  Transition  Model 

The  determination  of  the  location  of  transition  from  laminar  to  turbulent 
flow  is  one  of  the  most  critical  factors  in  the  success  of  many  computational  efforts  to 
predict  or  reproduce  physical  phenomena.  Yet,  even  in  today’s  age  of  supercomputers, 
this  area  of  research  remains  widely  neglected.  The  most  advanced  Navier-Stokes 
solvers  often  ignore  the  issue  entirely,  arbitrarily  declaring  the  entire  flow  to  be 
turbulent.  Others  make  the  effort  of  computing  a  transition  point,  at  which  the  flow 
instantaneously  changes  from  laminar  to  turbulent.  These  may  be  quite  reasonable 
approximations  for  many  applications,  especially  when  the  Reynolds  number  is  high. 
However,  there  are  still  a  number  of  important  flows  that  require  more  accuracy.  Until 
the  transition  mechanism  and  the  many,  varied  factors  which  can  affect  it  are  more  fully 
understood,  programmers  must  rely  on  the  traditional  engineering  approach  of  modeling. 

The  program  incorporates  a  transition  model  determined  by  Chen  and 
Thyson,  utilizing  a  transition  range  rather  than  a  point.  This  range  is  neither  fully 
laminar  nor  fully  turbulent.  It  is  a  region  of  iniermittency,  in  which  turbulent  spots 
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gradually  appear  with  progression  in  the  streamwise  direction.  This  feature  has  been 
shown  to  be  essential  for  low  Reynolds  number  flows.  The  convergence  of  the  entire 
boundary  layer  solution  is  very  sensitive  to  transition  related  factors,  such  as  the  input 
value  supplied  for  the  start  of  transition.  An  experimental  value  may  not  have  been 
measured  accurately  and  an  empirically  calculated  value  may  deviate  from  the  actual 
onset  of  transition.  If  the  code  does  not  run  and  all  other  input  has  been  verified  to  be 
correct,  a  solution  may  often  be  obtained  by  experimenting  with  slight  deviations  in  the 
transition  location  for  the  upper  surface  specified  as  input.  The  lower  surface  value  does 
not  exhibit  the  same  sensitivity. 

3.  Numerical  Techniques 

a.  Transformation  of  Airfoil  Coordinates 

The  x/c  and  y/c  coordinates  of  the  airfoil  are  supplied  as  part  of  the 
input,  called  xc  and  yc  in  the  program.  Starting  from  the  stagnation  point,  the  program 
redefines  these  coordinates  into  a  single  parameter  corresponding  to  a  surface  distance: 

Thus,  the  variable  x  used  internally  by  the  program  in  the  boundary  layer  equations  is 
this  surface  coordinate.  The  values  are  printed  under  the  heading  (S)  in  the  output. 

b.  Transformation  of  Variables 

The  well  known  Falkner-Skan  transformation  is  used  to  transform  the 


variable  y: 


(3.20) 


where  y  is  the  normal  coordinate  along  which  the  thickness  of  the  boundary  layer  is 
measured.  The  dimensionless  similarity  variable  eliminates  the  growth  of  the  boundary 
layer  in  laminar  flow  and  reduces  it  in  turbulent  flow.  This  enables  larger  steps  in  the 
streamwise  direction  and  improves  computational  efficiency.  The  x  transformation  is 
simply  a  scaling  by  the  reference  length,  usually  the  chord  for  a  airfoil,  so  that  ^=x/c. 
Since  the  surface  distance  x  is  used,  these  will  not  be  the  same  as  the  input  x/c 
coordinates. 

The  dimensionless  stream  function  f(x,Tj)  is  defined  by 

=  (3.21) 

Equations  3.11  and  3.18  and  the  boundary  conditions  may  be  rewritten  in  terms  of  the 
new  variables: 


-iff  .  { if^  -fK) 

I  0^ 


TJ  =0 


/=o 

/=i 


(3.22) 


where  ij*  is  the  transformed  boundary  layer  thickness  corresponding  to  5  and  Rl  is  the 
Reynolds  number  based  on  reference  velocity  u„  and  reference  length,  the  chord  for  an 
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airfoil.  The  prime  denotes  differentiation  with  respect  to  tj.  The  dimensionless  pressure 
gradient  parameter  m  is  defined  by 


m 


5  d(ujuj 


(3.23) 


uju,  dl 

The  velocity  components  v  and  v  are  related  to  the  dimensionless  stream  function  by 


jLAjiTx^^^r^ 


(3.24) 


c.  Keller  Box  Method 

Equation  3.22,  a  second-order  partial  differential  equation,  may  be 
solved  by  various  numerical  methods  such  as  the  Crank-Nicholson  or  Keller  Box 
methods.  The  latter  method,  depicted  in  Figure  3.1,  has  proven  to  be  efficient  for 
boundary  layer  calculations.  The  Keller  Box  method  first  requires  reformulating  higher 


Figure  3.1  Grid  Box  for  Centered-Difference  Approximations 
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order  equations  into  a  set  of  first  order  equations.  At  each  rectangular  grid  section,  these 
equations  are  approximated  using  centered-difference  derivatives,  averaging  values  at  the 
center  of  the  "box".  The  truncation  error  is  second  order.  The  resulting  implicit, 
nonlinear  difference  equations  are  linearized  by  Newton’s  method  and  solved  by  a  block 
elimination  method. 

Using  the  following  definitions, 

/'  =  «  tt'=v  (3.25) 

Equation  3.22  may  be  expressed  as  a  first  order  system: 

t1=0  M=0  f^fjx)  ,  M  =  1 

Note  that  the  u  and  v  in  Equations  3.25  and  3.26  are  not  the  velocity  components.  They 
are  two  new,  arbitrarily  selected  names  of  variables  for  the  expression  of  the  first  order 
system.  Since  the  program  was  coded  using  these  particular  variable  names,  the  present 
numerical  discussion  will  use  them  for  consistency. 

Referring  to  Figure  3.1,  the  grid  points  may  be  described  as 

=  n  =  1.2,...,N  5°=0 

^r'\j-\*K  i  =  l,  =  0  (3.27) 

where  the  capital  letters  N  and  J  are  the  maximum  number  of  grid  points  used  in  the 
streamwise  and  normal  directions,  respectively.  The  superscript  n  is  not  an  exponent, 
but  a  counter  just  like  the  subscript  j.  This  upper  and  lower  notation  allows  using  both 
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counters  on  the  same  variable  when  needed.  Considering  one  box  of  the  grid,  the  finite 
difference  approximations  of  Equation  3.25  may  be  written  for  midpoint  of  the  right  side, 
segment  P,P2,  using  centered  difference  derivatives; 


Jfll 

J  ~fj~t  _  _ n 


h: 


2 


^  *1 

hj  2 


(3.28) 


Equation  3.26  may  be  approximated  in  the  same  manner  for  the  midpoint  of  the  box  by 
centering  first  in  one  direction  then  the  other.  The  resulting  finite  difference  equation 
is 


T  ~ ^y-i';-i)  *  *1  )y-i4  'y-vt)  ■  ^y-v6 


where 


r;:^,  -  -L;:i .  -m 

2 
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W"  +  l  n  n  n 

ai  = - +  a  a2=w  +a" 


The  boundary  conditions  of  Equation  3.26  are  rewritten  at  as 


/o"=/w  «o’"0  «y"  =  l 


(3.30) 
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d.  Newton's  Method 


Equations  3.28,  3.29,  and  3.30  comprise  a  set  of  3J+3  equations  and 
3J+3  unknowns  (fj",  Uj“,  Vj",  where  j=  0,  1,  2,  J),  with  If  *,  Uj"*',  and  Vj“  '  known. 

Newton’s  method  is  applied  to  linearize  this  system.  The  method  assumes  that  an 
{q)proximate  solution  is  known,  either  from  the  preceding  iteration  cycle  or  from  the 
previous  streamwise  station.  Then  small  unknown  quantities  are  added  to  the 
approximate  solution.  Using  the  arbitrary  iteration  variable  i  (the  superscript  n  omitted 
for  clarity), 


uP^=Uj-*6Uj-  vj*‘=v/  +  8vy  (3.31) 

with  i=0  corresponding  to  known  values  at  the  previous  streamwise  station  (P  ‘),  these 
expressions  may  be  substituted  into  Equations  3.28  and  3.29  for  the  unknowns.  After 
dropping  higher  order  terms  of  8,  a  linear  system  of  equations  results: 


8«;  -  (6v/  +  5v/.,)  =  (r^j 

(s,)yav/  +  =  ir^)j 


(3.32) 


where  the  right  hand  sides  are 


~  ^-vi "  ~  *1  "  *2(“  ^y-vfc 

'v 


and  the  coefficients  are 
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(^6)y=-«2“/-i 


The  boundary  conditions  of  Equation  30  become 

5/o=0  6m<;=0  6«/=0 

These  equations  may  be  easily  identified  in  the  subroutine  COEF  of  bl2d.  Since  they 
may  be  arranged  into  a  block  tridiagonal  structure  in  matrix-vector  form,  the  subroutine 
SOLVE  uses  the  efficient  block  elimination  method  to  solve  for  the  small  5  quantities. 
The  iteration  of  Newton’s  method  continues  until  the  small  quantities  are  small  enough 
to  be  neglected. 

4.  Program  Modification  for  Boundary  Layer  Profiles 

At  each  station  along  the  airfoil  surface,  the  program  calculates  the  u 
velocities  for  each  value  of  ij  in  the  grid.  Eta  is  the  coordinate  in  the  normal  direction 
representing  the  transformed  boundary  layer  thickness.  In  order  to  retrieve  the  physical 
boundary  layer  thickness,  an  inverse  transformation  is  required: 
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(3.34) 


Since  the  kinematic  viscosity  shows  up  only  indirectly  in  the  non-dimensional 
form  of  the  Reynolds  number,  the  actual  equation  used  is 


y=Ti 


(3.35) 


where  RLspu^xZ/i,  or  equivalently,  UjX/t',  and  x  and  u^  are  used  in  the  non-dimensional 
forms  of  x/c  and  Uj/u*.  The  value  of  x  used  here  is  the  surface  distance. 

Plotting  the  shape  of  the  actual  velocity  profile  at  a  given  station  requires  the 
station  number  (NX),  the  corresponding  x  coordinate,  the  u  velocities,  and  the 
corresponding  y  values.  All  of  these  values  are  provided  by  the  original  program  or 
Equation  3.35.  Velocity  profiles  may  be  plotted  at  this  point,  but  only  the  shape  will  be 
revealed.  To  visualize  the  growth  of  the  boundary  layer,  the  height  of  the  boundary 
layer  at  each  station  is  needed.  This  may  be  determined  by  finding  where  the  u  velocity 
has  reached  freestream  velocity,  indicating  the  edge  of  the  boundary  layer. 
Computationally,  this  is  accomplished  by  allowing  u  to  reach  0.995  of  u.,  the  freestream 
velocity.  Even  though  the  remaining  u  values  in  the  grid  above  this  height  will  still  be 
calculated  by  the  program,  no  more  values  are  written  to  the  plotting  output  file. 
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5.  Program  Modification  for  Estimating  Transition  Location 

The  original  program  uses  input  values  to  specify  the  onset  of  transition.  The 
transitional  flow  region  is  then  calculated  using  the  Chen-Thyson  model,  shown  as  in 
Equation  3.17.  In  order  to  provide  an  initial  estimate  for  the  transition  location  when 
no  other  method  of  determination  is  available,  a  modification  using  Michel’s  criterion 
was  incorporated: 

p  =1.174(1+^^)^,  (3.36) 

Re  ' 

where  R^  =  \ij9/u  is  the  Reynolds  number  based  on  momentum  thickness  at  transition  and 
Re^r  is  the  Reynolds  number  based  on  the  transition  location. 

C.  RESULTS 

Studies  were  conducted  to: 

•  validate  the  program  by  comparing  to  known  data 

•  determine  the  effect  of  prescribed  onset  of  transition 

•  investigate  the  possible  occurrence  of  zero  or  negative  skin  friction  to  indicate  the 
start  of  a  laminar  separation  bubble  before  breakdown  of  the  code 

•  analyze  boundary  layer  velocity  profiles  on  the  airfoil  upper  surface 

•  evaluate  the  effects  of  changes  in  Reynolds  number 

•  assess  the  validity  of  obtaining  an  unsteady  boundary  layer  solution  by  extracting 
steady  velocities  from  unsteady  pressure  distributions 
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All  studies  presented  are  for  a  NACA  0012  airfoil  defined  by  100  points.  The  Reynolds 
number  is  540,000,  except  where  noted  in  the  validation  study  and  the  Reynolds  number 
effect  study. 

1.  Program  Validation 

In  order  to  ascertain  that  the  results  of  the  program  could  be  considered 
reliable,  an  initial  test  case  was  run  to  compare  with  published  data  [Ref.  3].  The  test 
conditions  were  an  AOA  of  0°  at  a  Reynolds  number  of  6,000,000.  The  first  results 
obtained  were  in  the  expected  range  but  the  curves  were  not  smooth.  TTie  problem  was 
traced  to  a  very  small  discrepancy  in  the  original  airfoil  coordinate  input  file  provided 
with  the  program  which  was  not  noticeable  when  the  file  was  checked  by  plotting.  After 
generating  new  airfoil  coordinates  and  running  them  through  panel  to  get  new  velocities, 
smooth  boundary  layer  results  were  obtained.  Figure  3.2  shows  the  computed  skin 
friction  coefficient  and  displacement  thickness  and  Figure  3.3  shows  the  published 
results.  The  plots  exhibit  excellent  agreement. 

2.  Transition  Onset  Location 

In  the  next  study,  the  convergence  of  bl2d  at  a  lower  Reynolds  number  of 
540,000  was  investigated.  The  input  value  for  the  location  of  the  start  of  the  transition 
range  was  found  to  be  the  most  critical  factor.  Starting  at  an  angle  of  attack  of  0“, 
transition  values  obtained  from  incompbl  (Chapter  4)  were  initially  used.  The  program 
converged  to  a  reasonable  solution  as  determined  by  viewing  plots  of  the  various  output 
files.  At  2®  however,  the  program  would  not  converge  using  the  estimated  transition 
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value.  Experimentation  with  this  input  parameter  showed  that  moving  it  back  usually 
made  the  convergence  even  worse,  but  in  moving  it  forward,  a  point  was  reached  where 
a  solution  could  be  obtained.  Moreover,  this  solution  was  one  that  could  be  reasonably 
expected  based  on  the  previous  solution  at  0°.  Similar  experimentation  was  performed 
in  AOA  increments  of  2°  until  excessive  separation  at  high  angle  of  attack  caused  the 
code  to  break  down.  The  same  study  was  also  conducted  later  using  the  version  of  bl2d 
modiried  to  make  transition  estimates.  Table  3. 1  shows  a  comparison  between  the  values 
estimated  by  the  two  programs  and  the  value  actually  needed  for  the  program  to 
converge. 


Table  3.1  COMPARISON  OF  TRANSITION  ONSET 


AOA 


INCOMPBL 


BL2D 


convergence 


%  diff 


The  first  significant  point  to  note  is  that  the  transition  values  produced  by  the 


modified  version  of  bl2d  are  consistently  close  to  the  values  produced  by  incompbl. 


showing  that  the  criterion  has,  in  fact,  been  programmed  correctly.  The  differences  can 
be  attributed  to  several  factors: 

•  slightly  different  values  are  input  into  the  criterion  equation  for  each  program 

•  incompbl  outputs  the  x/c  value  of  the  nodal  point  nearest  the  calculated  transition 
onset,  as  opposed  to  the  actual  value 

•  if  the  onset  of  transition  is  located  inside  a  separation  bubble  by  the  initial 
calculation,  incompbl  arbitrarily  moves  it  to  the  start  of  the  bubble 

The  third  column  shows  how  far  forward  the  transition  point  was  moved  to  obtain 
convergence,  where  the  first  AOA  did  not  exhibit  sensitivity.  The  last  column  shows  the 
percent  difference  between  the  best  estimate  and  the  value  required  for  convergence. 
Most  were  fairly  close,  with  even  the  worst  case  less  than  30%  forward  of  the  first 
estimate.  This  sets  a  reasonable  bound  for  necessary  experimentation  with  the  transition 
location. 

3.  Laminar  Separation 

a.  Sfdn  Friction  Coefficient  and  Displacement  Thickness 

Figures  3.4-3. 10  show  the  progression  of  skin  friction  coefficient  (Cf) 
and  displacement  thickness  (5*)  as  the  AOA  ranges  from  0  to  12  degrees.  The  transition 
onset  may  be  observed  as  the  point  where  Cf  reaches  as  minimum  then  dramatically 
increases,  indicating  the  change  from  laminar  to  turbulent  flow.  The  transition  point 
moves  forward  as  the  angle  of  attack  increases.  The  minimum  value  of  Cf  decreases 
with  increasing  angle  of  attack.  When  the  Q  reaches  zero,  separation  is  indicated. 
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Cf,  delta 


Figure  3.4  BL2D;  C,  and  6*,  NACA  0012,  AOA=0%  Re=540,000 
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Figure  3.10  BL2D:  Cf  and  5*,  NACA  0012,  AOA  =  12%  Re=540,000 

Knowing  that  the  boundary  layer  equations  break  down  when  separation 
occurs,  the  anticipated  information  of  this  study  was  the  x/c  location  of  transition  for  as 
many  angles  of  attack  as  possible  before  the  Cf  became  negative,  where  it  was  assumed 
the  program  would  not  run.  A  beneficial  discovery  of  this  study  is  the  ability  of  the 
program  to  recover  from  mild  amounts  of  separation  with  meaningful  results.  At  4°,  the 
first  encounter  with  separation  may  be  observed  as  the  Cf  just  dips  below  zero.  The 
pattern  over  the  remaining  airfoil  surface  suggests  a  separation  "bubble"  after  which  the 
flow  reattaches,  as  opposed  to  near-stall  separation.  The  program  exhibited  this  recovery 
behavior  all  the  way  to  10°.  The  final  plot  at  12°  shows  that  even  though  a  solution  was 
produced,  convergence  was  not  attained  and  the  results  were  meaningless  due  to  the 
greater  amount  of  separation. 
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b.  Boundary  Layer  Velocity  Profiles 

To  complement  the  skin  friction  study  and  to  further  investigate  the 
pattern  of  laminar  separation,  the  program  was  modified  to  calculate  and  produce  plotting 
output  for  velocity  profiles  at  evenly  spaced  intervals  along  the  top  surface  of  the  airfoil. 
Results  are  shown  in  Figures  3. 1 1-3. 16  for  an  AOA  range  of  0  to  10  degrees.  Boundary 
layer  growth  is  evident  as  the  angle  of  attack  increases.  Furthermore,  the  region  most 
prone  to  separation,  as  revealed  by  the  point  of  inflection  in  the  velocity  profile,  moves 
forward  with  increasing  angle  of  attack,  confirming  the  indications  of  the  skin  friction 
plots. 
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Figure  3.11  BUD:  Velocity  Profiles,  NACA  0012,  AOA=0“,  Re=540,000 
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Figure  3.12  BUD:  Velocity  Profiles,  NACA  0012,  AOA=2%  Re=540, 
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Figure  3.13  BUD:  Velocity  Profiles,  NACA  0012,  AOA=4°,  Re=540,000 
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Airfoil  Upper  Surface  Station 


Figure  3.14  BL2D:  Velocity  ProHles,  NACA  0012,  AOA=6%  Re=540, 
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Figure  3.15  BUD:  Velocity  ProHies,  NACA  0012,  AOA=8%  Re=540,000 
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Figure  3.16  BL2D:  Velocity  Profiles,  NACA  0012,  AOA=10“,  Re=540,000 

4.  Reynolds  Number  Changes 

The  direct  boundary  layer  code  was  run  for  Reynolds  numbers  of  540,000, 
750,000,  and  1,000,000.  A  representative  sample  of  the  results  at  2°  is  presented  in 
Figure  3.17,  showing  a  comparison  of  the  skin  friction  coefficients.  As  the  Reynolds 
number  increases,  two  effects  may  be  observed.  First,  the  transition  point  moves 
forward,  so  there  is  more  turbulent  flow.  This  occurs  on  both  surfaces,  although  only 
the  upper  surface  is  plotted.  In  addition,  the  minimum  value  of  Cf  increases.  Both 
effects  suggest  that  separation  is  less  likely  to  occur  as  Reynolds  number  increases,  all 
other  conditions  being  constant. 
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NACA  0012  ACA=2  deg.  Variation 
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Figure  3.17  BL2D:  Effects  of  Reynolds  Number  Changes 
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5.  Unsteady  Boundary  Layers 

By  correlating  an  unsteady  pressure  distribution  with  a  similar  steady  one,  the 
question  arises  as  to  whether  the  steady  direct  boundary  layer  code  can  process  the  data 
in  the  form  of  a  steady  velocity  distribution  and  produce  results  that  correspond  to  the 
unsteady  case.  In  the  studies  of  Neace  [Ref.  4],  it  was  shown  that  an  unsteady  pressure 
distribution  at  a  certain  angle  of  attack  has  a  closely  matching  steady  pressure  distribution 
at  an  angle  of  attack  approximately  2°  lower.  For  example,  10. 14°  unsteady  corresponds 
to  8°  steady.  An  unsteady  panel  method  program  called  U2DIIF  was  used  to  calculate 
unsteady  pressure  distributions  for  ramp  motion.  With  the  hypothesis  that  the  boundary 
layer  characteristics,  steady  or  unsteady,  are  driven  by  the  pressure  distribution,  the 
pressures  were  converted  to  "steady"  velocities  for  the  unsteady  angles  of  attack  using 
Bernoulli’s  equation.  Steady  velocities  for  steady  angles  of  attack  were  computed  with 
the  program  panel.  In  Figure  3.18,  the  comparison  between  steady  and  unsteady 
velocity  distributions  reveals  excellent  agreement  on  the  upper  surface  and  fair  agreement 
on  the  lower  surface.  Since  velocity  is  the  main  input  to  bl2d,  and  the  upper  surface  is 
much  more  critical  (for  positive  angles  of  attack),  the  strong  correlation  suggests  that  the 
steady  boundary  layer  code  may  indeed  be  able  to  produce  results  for  the  unsteady  case. 
Many  attempts  were  made  to  obtain  such  results  using  all  of  the  methods  explained  in 
other  sections  to  facilitate  convergence.  No  solution  was  found  for  any  of  the  cases 
investigated,  however.  Apparently,  the  differences  on  the  lower  surface  had  more  effect 
than  anticipated.  Future  investigations  could  include  a  modified  velocity  input,  where 
the  unsteady  upper  surface  is  spliced  with  the  corresponding  steady  lower  surface. 
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Smoothness  should  be  ensured  in  the  connecting  areas.  A  compressibility  correction  may 
also  make  a  difference. 

D.  USER’S  GUTOE  TO  BUD 

1.  Output  from  the  Program  PANEL 

Three  files  are  output  from  the  panel  code  described  in  the  previous  chapter. 
The  file  vel.dat  is  simply  a  printout  of  x/c  and  velocity  relative  freestream  data,  mainly 
for  plotting  purposes  if  desired.  Even  though,  during  computation  in  the  program,  panels 
on  the  lower  surface  downstream  of  the  stagnation  point  were  taken  to  have  negative 
"directional"  velocities,  the  output  file  correctly  displays  the  positive  "physical" 
velocities.  The  file  cp.dat  lists  x/c  and  pressure  distribution  data.  The  other  output  file, 
called  bl2d.dat,  is  generated  to  be  compatible  with  the  program  blld  as  input.  An 
example  input  file  for  the  modified  program  is  included  in  Appendix  A. 

2.  Input  Description 

The  first  line  of  bl2d.dat  consists  of  Reynolds  number,  transition  location 
(x/c)  on  the  upper  surface,  and  transition  location  on  the  lower  surface.  The  next  line 
indicates  the  number  of  points  and  the  i  value  of  the  location  of  the  stagnation  point.  For 
the  modified  version,  an  additional  input  value  is  the  transition  flag.  This  indicates  if  the 
run  is  an  initial  estimate  (0)  or  the  transition  values  are  fairly  well  determined  and  valid 
boundary  layer  calculations  are  desired.  The  remaining  lines  are  identical  to  the  velocity 
output  file,  except  that  y/c  values  of  the  airfoil  coordinates  are  also  included.  Note  that 
these  velocities  correspond  to  a  particular  angle  of  attack,  the  one  that  was  specified 


53 


when  panel  was  run.  The  actual  value  of  the  angle  of  attack  is  not  listed  separately  in 
the  file. 

The  file  bl2d.dat  may  be  edited,  replacing  the  originally  specified  values  in 
the  first  few  lines  with  new  values  for  the  desired  flow.  The  Reynolds  number  is  based 
on  the  appropriate  reference  length,  usually  the  chord  for  an  airfoil.  The  transition 
locations  may  be  obtained  from  experimental  data  or  from  a  calculation  method.  In  the 
modified  version,  Michel’s  method  may  be  used  to  initially  estimate  the  location  of 
transition  onset.  If  no  information  is  known  before  using  the  program,  the  transition 
values  specified  for  the  first  run  should  be  large  enough  to  be  downstream  of  the  actual 
transition  points,  yet  not  so  large  that  the  program  will  not  converge.  The  arbitrary 
values  supplied  by  the  program  panel  are  0.8  for  the  upper  surface  and  0.999  for  the 
lower  surface.  If  the  transition  flag  at  the  end  of  the  next  line  in  the  input  file  is  set  to 
0,  a  run  of  the  program  will  show  estimates  for  transition  location  on  the  screen  but  these 
values  will  not  be  used  for  the  boundary  layer  calculations.  The  input  file  should  be 
edited,  replacing  the  initial  downstream  transition  values  with  the  estimates  and  changing 
the  transition  flag  to  1.  If  the  program  does  not  converge  with  these  values,  it  will  be 
necessary  to  experiment  with  slight  deviations  in  the  upper  surface  value.  Si..ce  the 
transition  value  calculated  is  only  an  estimate,  this  experimentation  procedure  is  the  rule 
rather  than  the  exception,  especially  at  higher  angles  of  attack. 

For  the  stagnation  location,  it  is  important  to  not  simply  take  the  i  value  from 
the  velocity  output  which  corresponds  to  the  velocity  closest  to  zero.  In  order  to  work 
properly,  bl2d  usually  requires  the  i  value  for  the  first  point  after  the  stagnation  point, 
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where  points  are  numbered  in  a  clockwise  direction  from  the  trailing  edge.  With  only 
positive  velocities  to  inspect,  it  is  impossible  to  tell  where  this  occurs.  The  program 
panel,  however,  uses  the  negative  directional  velocities  in  its  calculations  and 
automatically  determines  the  proper  i  value  to  send  to  the  output  file  bl2d.dat.  If 
convergence  cannot  be  obtained  by  varying  the  transition  location,  it  may  help  in  some 
cases  to  decrease  the  predetermined  i  value  of  the  stagnation  location  by  one. 

3.  Program  Operation 

When  all  necessary  values  have  been  changed,  rename  the  input  file 
appropriately  for  reference,  such  as  bl5.dat  for  an  AOA  of  5°.  Subsequent  runs  of  the 
program  panel  for  other  AOA’s  will  overwrite  the  output  file  bl2d.dat.  When  ready  to 
run  bl2d,  copy  the  desired  input  file  to  bl2d.dat,  the  required  input  file  name.  To  run 
the  program,  type; 

>bl2d 

Convergence  may  most  easily  be  observed  by  plotting  output  files  such  as  the  skin 
friction  coefficient.  Modifications  to  various  parameters  as  explained  in  the  previous 
section  may  facilitate  convergence.  In  some  cases,  however,  such  as  an  unusually  shaped 
airfoil,  a  highly  cambered  airfoil,  or  a  standard  airfoil  with  a  faulty  input  file, 
convergence  may  not  be  attainable.  Additionally,  all  airfoils  at  a  high  enough  angle  of 
attack  will  cause  the  program  to  break  down,  as  the  direct  boundary  layer  method  cannot 
handle  significant  separation. 
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4.  Output  Description 


The  output  to  the  screen  indicates  the  progress  of  the  program  as  it  runs.  In 
the  modified  program,  the  first  estimates  for  transition  locations  will  also  print  to  the 
screen  when  the  transition  flag  is  set  to  0.  The  remaining  output  in  this  case  is  not 
applicable  and  the  program  should  be  rerun  with  new  transition  values. 

The  standard  output  file  is  named  bl2d.out.  The  first  line  repeats  the  values 
of  Reynolds  number  and  transition  locations  that  were  supplied  as  input.  Next,  a 
summary  of  boundary  layer  solutions  is  presented  for  ISF=1,  the  upper  surface 
downstream  of  the  forward  stagnation  point.  Reading  across,  the  data  consists  of  station 
or  point  number  (NX),  x/c  value  (XC),  distance  over  the  airfoil  surface  starting  at  the 
forward  stagnation  point  for  that  AOA  (S),  dimensionless  wall  shear  parameter  (VW), 
skin  friction  coefficient  (CF),  displacement  thickness,  8*  (DLS),  and  6,  the  momentum 
thickness  (THT).  The  same  data  is  listed  for  ISF=2,  the  lower  surface  at  that  AOA. 
It  is  important  to  realize  that  an  otherwise  converged  solution  may  still  have  several 
highly  divergent  values  for  all  parameters  near  the  trailing  edge.  This  does  not  invalidate 
the  whole  solution.  Since  the  current  investigation  concerns  primarily  the  leading  edge 
area,  the  last  few  values  may  be  ignored. 

Several  modifications  were  made  to  the  program  to  allow  additional 
informative  output.  The  skin  friction  coefficient  and  displacement  thickness  are  printed 
in  the  output  files  cf.dat  and  dls.dat,  respectively,  with  x/c  values  for  plotting. 
Individual  boundary  layer  velocity  profiles  at  every  five  stations  along  the  airfoil  upper 
surface  are  printed  in  various  output  files  as  follows: 


56 


NX =5  FOR031.DAT 

NX =10  FOR032.DAT 

NX=15  FOR033.DAT 

NX=20  FOR034.DAT 

NX =25  FOR035.DAT 

NX =30  FOR036.DAT 

NX =35  FOR037.DAT 

NX =40  FOR038.DAT 

NX =45  FOR039.DAT 

etc. 

where  station  (NX)  1  is  the  first  point  after  the  stagnation  point  and  numbering  increases 
across  the  top  of  the  airfoil.  The  example  shown  is  for  a  100  point  airfoil  with  50  points 
across  the  upper  surface.  For  a  different  number  of  input  points,  the  number  of  output 
files  will  adjust  accordingly.  The  station  nearest  the  trailing  edge  is  not  output  because 
the  results  are  often  divergent.  A  comprehensive  summary  of  these  files  is 
simultaneously  stored  in  the  output  file  FOR060.DAT  for  ease  in  plotting  all  of  the 
profiles.  An  additional  output  file,  FOR055.DAT,  contains  plotting  data  for  drawing 
evenly  spaced  lines  across  the  velocity  profiles.  This  is  a  visual  effect  only  and  is  not 
necessary  if  not  desired. 

While  there  are  many  ways  to  display  this  data,  a  command  file  compatible 
with  gnuplot  called  profile  shows  the  data  to  best  advantage.  A  copy  of  this  file  is 
provided  in  Appendix  A.  It  automatically  plots  all  the  profiles  on  one  plot,  showing 
boundary  layer  growth  and  the  changing  slope  of  the  profiles,  indicating  when  separation 
occurs.  Using  gnuplot  in  the  directory  where  the  data  files  and  the  file  profile  reside, 
type: 

>loadl  ’’profile” 
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Modifications  may  be  made  to  the  resulting  plot  using  gnuplot  commands,  Appendix  B, 
or  editing  the  command  file  in  a  separate  window. 

5.  PC  Version 

To  increase  its  accessibility,  bl2d  was  also  converted  to  a  version  compatible 
with  personal  computers.  Additional  programs  included  on  the  Boundary  Layer 
Analysis  disk  are  a  PC  version  of  panel,  an  airfoil  point  generation  program  called 
airfoil,  a  PC  version  of  gnuplot,  and  the  command  plotting  file  profile.  Most  of  the 
instructions  are  the  same  for  this  version,  but  there  are  a  few  differences.  The  programs 
panel  and  bl2d  are  restricted  to  100  or  less  airfoil  points  due  to  array  limitations  in  the 
PC  FORTRAN  compiler.  An  automatic  rerun  feature  was  incorporated  into  bl2d  for  the 
case  of  initial  transition  estimates,  where  the  first  run  internally  restarts  using  the 
calculated  estimates  for  boundary  layer  calculations.  Finally,  the  velocity  profiles  are 
output  only  to  a  comprehensive  file  called  profilel  .dat  instead  of  FOR060.DAT.  There 
are  no  individual  velocity  profile  files.  The  horizontal  line  file  called  FOR055.DAT  in 
the  UNIX  version  is  called  profile2.dat  for  the  PC  version. 
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rv.  VISCOUS-IN  VISCID  INTERACTION  CODE 


A.  THEORY 

The  direct  boundary  layer  code,  described  in  the  previous  chapter,  calculates  a 
displacement  thickness  for  a  prescribed  pressure  distribution  (or  equivalently,  a  velocity 
distribution,  for  incompressible  flow).  As  the  name  implies,  it  is  a  direct  calculation 
involving  one  pass,  thus  very  little  time  is  required.  Another  method,  known  as  the 
inverse  boundary  layer  method,  calculates  a  pressure  (or  velocity)  distribution  for  a 
prescribed  displacement  thickness.  The  displacement  thickness  represents  an  "effective 
body"  as  far  as  the  flow  is  concerned.  Iteration  is  required  and  the  method  requires 
considerably  more  calculation  time.  The  main  advantage  of  the  inx  erse  method  is  the 
ability  to  calculate  through  regions  of  separation. 

The  simplest  viscous-inviscid  interaction  method  divides  the  flowfield  into  an  inner 
viscous  region  where  boundary  layer  calculations  are  performed  and  an  outer  inviscid 
region  where  potential  flow  analysis  prevails.  The  solutions  are  then  iterated  until  they 
match  along  the  dividing  line.  This  method  has  "weak"  interaction  because  the  only 
exchange  of  information  is  along  the  boundary. 

In  the  strong  interaction  method,  both  the  pressure  and  displacement  thickness  are 
treated  as  unknowns  and  are  solved  simultaneously  with  successive  sweeps  over  the 
airfoil.  The  external  boundary  condition  for  the  boundary  layer  equation  at  the  outer 
edge  of  the  viscous  region  is 
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u^ix)  =uj^{x) + 5  w^(x)  (4-1) 

where  Ue°(x)  is  the  inviscid  velocity  over  the  airfoil  and  5Uj(x)  is  the  perturbation  due  to 
viscous  effects,  expressed  as 

6u/x)=—f  ^[u(a)5*(a)]  —  da  (4.2) 

itJ  aa'-  ^  ^x-a 

‘s 

where  d(Ue6*)/da  is  the  blowing  velocity.  Equations  4.1  and  4.2  comprise  the 
interaction,  or  coupling  law.  The  interaction  takes  place  between  x,  and  x^.  The 
integral  term  is  known  as  the  Hilbert  integral,  in  which  the  displacement  effect  is 
modeled  by  source/ sink  distributions  using  potential  flow  theory. 

B.  COMPUTER  CODE 
1.  Overview 

The  viscous-inviscid  program  incorporates  a  self-contained  panel  code  based 
on  the  Hess-Smith  method,  a  boundary  layer  calculation  routine,  and  an  interaction 
scheme.  The  inviscid  panel  method  is  similar  to  the  program  panel  described  in  Chapter 
2;  however,  the  boundary  condition  of  zero  normal  velocity  on  the  surface  of  the  body 
is  replaced  by  a  blowing  velocity  determined  from  the  boundary  layer  calculations.  This 
blowing  velocity  is  used  to  represent  the  viscous  effects  on  the  inviscid  flow.  The 
boundary  layer  method  is  similar  to  the  program  bl2d  with  some  modifications  that 
account  for  the  presence  of  the  wake  and  for  low  Reynolds  number  flows.  The  viscous- 
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inviscid  code  used  for  the  investigation  was  developed  by  T.  Cebeci.  Run  time  on  the 
Stardent  computer  is  about  five  minutes. 

2.  Models 

a.  Turbulence  Model 

As  in  the  direct  boundary  layer  program,  the  eddy-viscosity  formulation 
of  Cebeci  and  Smith  is  used,  with  separate  formulas  for  the  inner  and  outer  regions. 
This  model  in  this  program  has  the  additional  features  of  low  Reynolds  number  effects 
and  a  wake  flow  model.  The  modified  equations  are  expressed  as: 


where 


(€«)o  =  «M*Y,rY 


du 
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and  the  transition  length  is  represented  by 


0^=213  Pogfe^ -4.7323] 


For  the  wake  flow,  the  eddy-viscosity  formulation  is 


(4.5) 


(4.6) 


where  (e,J*  specifies  the  eddy  viscosity  of  the  far  wake,  taken  as  the  maximum  of  the 
lower  and  upper  wake  eddy  viscosities: 


>1* 

(Ol-OOM  /  (».-»)<*)' 

(4.7) 

«• 

(€X  =  0064  J  (u^-u)dy 

ywim 


with  the  location  where  u=u„un. 

b.  Transition  Model 

The  program  uses  an  empirical  formula  called  Michel’s  criterion  to 
calculate  a  first  approximation  to  the  transition  location  on  the  upper  and  lower  airfoil 
surfaces.  It  is  expressed  as  a  relationship  between  the  Reynolds  numbers  based  on 
momentum  thickness  and  on  the  x  (surface)  location  of  transition: 


=1.174(1 (4.8) 

'  Re  ^ 

This  is  the  same  equation  used  in  the  modified  version  of  bl2d.  The  program  incompbl, 
however,  sometimes  adjusts  the  resulting  value.  If  an  area  of  separation  with  subsequent 
reattachment  is  calculated,  the  transition  onset  is  moved  to  the  beginning  of  the 
separation  bubble.  The  transition  location  that  is  printed  in  the  output  file  always 
corresponds  to  a  nodal  point.  In  many  cases,  the  program  will  produce  a  fairly  accurate 
result.  Sometimes,  however,  a  refinement  may  be  needed.  This  may  be  done  by 
experimenting  with  small  shifts  around  the  calculated  value,  as  described  in  Chapter  3 
for  the  direct  boundary  layer  program.  The  experimentation  may  be  implemented  by 
using  the  transition  specification  option,  which  overrides  calculation.  Input  options  are 
described  in  the  User’s  Guide  section. 

A  more  advanced  technique,  suggested  by  Cebeci,  is  the  e“  method, 
which  makes  use  of  linear  stability  theory.  A  separate  stability/transition  code 
incorporating  this  method  is  run  using  the  output  of  the  viscous/inviscid  code  (Michel’s 
criterion)  as  the  first  estimate.  The  new  value  is  supplied  as  input  to  the  first  code,  this 
time  overriding  Michel’s  criterion.  This  type  of  manual  iteration  continues  until 
convergence,  usually  within  three  to  four  cycles.  At  the  present  time,  determination  of 
all  required  input  values  and  analysis  of  the  output  require  significant  experience  and 
judgment,  prohibiting  a  programmed  link  between  the  two  codes  until  further  refinement 
is  accomplished.  The  e"  method  was  not  used  in  this  study. 
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As  in  the  program  bl2d,  the  program  incompbi  uses  the  Chen-Thyson 
transition  range  model  to  calculate  the  length  of  the  transition  region.  Equation  4.5 
shows  the  modified  version  of  this  model.  The  program  incorporates  an  improvement 
for  the  transition  length  parameter.  may  be  identified  by  reducing  the  (3/C^ 
term  to  the  form  (1/G.^),  with  G.^=CV3.  In  the  original  model,  the  constant  C  has  a 
recommended  value  of  60,  resulting  in  G.^„=1200.  Whereas  a  value  of  1200  may  work 
well  for  large  Reynolds  numbers,  values  from  20  to  80  have  been  shown  to  be  most 
successful  in  low  Reynolds  number  flows  where  separation  bubbles  exist  [Ref.  5].  The 
program  determines  an  appropriate  value  using  an  empirical  correlation  formula  in  the 
form  of  C^,  also  shown  in  Equation  4.5. 

3.  Numerical  Techniques 
a.  Hilbert  Integral 

Equation  4.1,  containing  the  Hilbert  integral,  may  be  approximated  in 
discretized  form  as 


=  « «V)  +  5^  “) 

where  u/fx)  corresponds  to  the  inviscid  velocity  distribution  which  contains  the 
displacement  thickness  effect  ($*)“  and  Cjj  is  a  matrix  of  interaction  coefficients  which  are 
functions  of  geometry  only. 
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b.  FLARE  Appmximation 

In  regions  of  recirculating  flow,  such  as  a  separation  bubble,  numerical 
stability  difficulties  may  be  encountered.  The  FLARE  approximation,  due  to  Fliigge- 
Lotz  and  Reyhner,  neglects  the  longitudinal  convective  term  u(du/dx)  in  the  region  of 
negative  u  velocity. 

C.  RESULTS 

Studies  were  conducted  to: 

•  Determine  the  effect  of  including  the  wake  in  the  calculations 

•  Investigate  the  possible  occurrence  of  negative  skin  friction  and  to  determine  its 
significance 

•  Compare  the  results  with  the  direct  boundary  layer  code 

Since  the  viscous-inviscid  method  has  the  ability  to  calculate  through  regions  of 
separation,  additional  information  was  anticipated  beyond  that  provided  by  the  direct 
boundary  layer  method. 

1.  Wake  Calculations 

Since  one  of  the  input  options  is  for  the  inclusion  of  wake  calculations, 
investigations  were  performed  to  determine  its  effect.  Calculations  may  be  limited  to  the 
airfoil  surface  only,  or  a  grid  extending  into  the  wake  region  may  also  be  used.  A 
representative  sample  of  the  results  is  shown  in  Figure  4.1,  depicting  the  skin  friction 
coefficient  for  three  angles  of  attack.  At  2®,  the  main  difference  is  a  movement  aft  of 
the  transition  point,  where  the  remainder  of  the  curve  maintains  its  original  shape.  At 
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6°,  the  transition  point  also  moves  forward;  however,  the  shape  of  the  turbulent  section 
is  slightly  altered.  At  10°,  the  transition  point  is  unchanged.  More  significant  changes 
occur  progressing  across  the  top  of  the  airfoil,  leading  to  a  delay  in  the  point  where  the 
skin  friction  falls  below  zero.  Thus,  the  primary  effect  of  wake  flow  is  to  reduce  flow 
separation  on  the  airfoil,  especially  important  at  higher  angles  of  attack.  This  will  allow 
calculations  to  continue  when  convergence  may  not  have  been  reached  otherwise.  These 
results  are  consistent  with  those  found  by  Cebeci  [Ref.  6].  Therefore,  the  wake  was  used 
in  all  further  studies  using  this  code. 

2.  Laminar  Separation 

a.  Comparison  with  Direct  Boundary  Layer  Method 

The  program  incompbl  was  run  for  the  same  conditions  as  the  previous 
study  with  bI2d,  a  NACA  00X2  airfoil  with  a  Reynolds  number  of  540,000.  The  angle 
of  attack  was  increased  in  two  degree  increments.  Figures  4.2  to  4.7  compare  the  skin 
friction  results  for  both  programs.  The  low  angles  of  attack  show  excellent  agreement 
in  the  laminar  region.  As  the  angle  of  attack  increases,  the  most  notable  difference  is 
the  absence  of  Cf  <  0  for  the  viscous-inviscid  method.  Since  the  direct  boundary  layer 
code  failed  to  converge  with  separation  greater  than  that  produced  at  10°,  no  comparison 
could  be  performed  beyond  this  point. 
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Figure  4.6  Cf  Comparison,  NACA  0012,  AOA=8°,  Re=540,000 
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b.  Investigation  of  Higher  Angles  of  Attack 


In  Figures  4.8  to  4. 14,  the  angle  of  attack  for  the  viscous-inviscid  code 
was  increased  further  in  search  of  separation,  indicated  by  Cf<0.  Finer  increments  were 
used  when  separation  appeared  to  be  imminent.  Separation  did  not  occur  until  13.3° 
and  was  gone  by  13.4°.  A  final  plot  at  14°  shows  that  the  remaining  trend  is  an  increase 
in  the  minimum  value  of  Cf.  The  amount  of  separation  at  13.3°  was  minuscule,  as  only 
one  data  point  fell  below  zero  with  a  value  of  -0.00065. 
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Figure  4.13  INCOMPBL:  C„  N  AC  A  0012,  AOA  =  13.4%  Re =540,000 
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D.  USER’S  GUIDE  TO  EVCOMPBL 

1.  Required  Files 

The  files  required  in  a  directory  to  run  this  program  are; 

•  incompbl*;  a  compiled,  executable  program  compatible  with  the  computer  being 
used  (arbitrary  name  if  starting  from  source  code) 

•  FOR001.DAT  (Stardent)  or  fort.l  (IRIS);  first  input  file  containing  airfoil  related 
data  (required  name) 

•  incompbl.dat;  second  input  file  containing  flow  and  calculation  information 
(arbitrary'  name) 

2.  UNIX  FORTRAN 

There  are  three  ways  to  get  data  from  external  input  files  into  a  FORTRAN 
program  on  a  UNIX  based  system.  The  first  uses  an  OPEN  statement  which  associates 
a  specific  input  file  name  with  a  file  number,  just  as  on  a  personal  computer.  This 
method  is  not  used  in  this  program.  The  second  method  has  READ  statements  using  file 
numbers  not  declared  by  an  OPEN  statement.  In  this  case,  a  default  file  name  such  as 
FOR(X)l.DAT  is  assumed  for  a  READ  statement  referencing  a  file  number  of  one,  for 
example.  The  program  uses  this  method  for  the  first  input  file  (hence  the  required  name) 
as  well  as  many  output  files,  which  the  user  can  modify  or  add  to  in  the  source  code 
before  compilation.  The  third  method,  which  may  be  used  only  once  in  a  program,  is 
to  specify  an  input  file  on  the  command  line  at  execution  time.  The  file  can  have  any 
name  chosen  by  the  user  as  long  as  that  file  has  data  in  the  correct  format  expected  by 
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the  program.  READ  statements  using  a  file  number  of  five  assume  interactive  input 
from  the  keyboard  during  execution  or  an  input  file  name  in  the  command  line. 

3.  Starting  from  the  Source  Code 

If  the  compiled  version  is  not  available  or  it  is  not  certain  which  source  code 
file  corresponds  to  the  compiled  file  that  exists,  the  source  code  is  the  best  starting  point. 
Otherwise,  proceed  with  execution  procedures  in  the  next  section. 

The  source  code  incompbi.f  may  be  obtained  from  either  the  Stardent  or  IRIS 
computer  and  transferred  from  one  to  the  other  as  needed  using  the  ftp  utility.  Once 
compiled,  however,  the  program  will  be  computer  specific,  even  if  the  executable 
versions  have  the  same  name. 

Modifications  may  be  made  to  the  source  code  first  if  desired.  For  example, 
a  WRITE  (32,*)  statement  could  be  inserted  to  output  specific  data  to  a  file  with  a 
default  name  of  FOR032.DAT  (Stardent)  or  fort. 32  (IRIS).  This  method  is  often  used 
to  quickly  obtain  files  of  x/c  vs.  velocity,  pressure,  or  other  parameters  for  plotting. 

To  compile  the  FORTRAN  source  code  for  the  Stardent,  type: 

>  fc  -02  -o  incompbl  incompbi.f 

The  term  fc  is  for  FORTRAN  Compile.  The  term -02  permits  vector  optimization.  The 
program  will  run  without  it,  but  a  warning  notice  will  be  issued  after  compilation.  Note 
the  use  of  the  letter  "O",  not  the  number  zero,  "0".  DO  NOT  use  the  option  -03 
instead.  This  option  is  for  parallel  processing,  a  feature  not  currently  incorporated  on 
the  NFS  Stardent.  The  program  may  appear  to  successfully  compile  and  run,  but  there 
will  usually  be  errors  in  the  output.  The  term  following  the  -o  is  the  name  of  the  output 
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executable  program.  Any  name  may  be  used.  The  program  name  will  be  displayed  with 
an  asterisk  (*)  following  it  in  the  directory  listing  to  indicate  that  it  is  an  executable 
program.  Finally,  the  source  code  or  codes  are  listed  in  order,  only  one  being  used  in 
this  case.  To  compile  on  the  IRIS,  type 
>  f77  -03  -0  incompbl  incompbl.f 

The  optimization  levels  have  different  meanings  on  the  IRIS  and  -03  is  the  correct 
parameter. 

During  the  compilation  process,  a  file  incompbl. o,  called  an  object  file,  is 
produced  appears  in  the  directory  listing.  This  file  is  not  needed  in  this  application  and 
may  be  deleted. 

4.  Input  File  Editing 

The  first  input  file  pertains  to  the  panel  method  part  of  the  program  and  is 
called  FOR001.DAT  (Stardent)  or  fort.l  (IRIS).  A  sample  file  is  included  in  Appendix 
A.  The  first  line  is  simply  a  number  telling  how  many  of  the  following  lines  are  for 
comments.  The  next  few  lines  contain  the  comments,  such  as  the  type  airfoil  being 
analyzed,  the  date  of  the  test,  or  any  other  information  useful  to  the  user.  The  next 
group  of  data  consists  of  the  angle  of  attack  (ALPI),  the  x/c  location  of  the  pivot  about 
which  the  airfoil  rotates  to  a  new  angle  of  attack  (PIVOT),  and  the  number  of  panels 
defining  the  lower  and  upper  surfaces  (NLOWER  and  NUPPER,  respectively).  Finally, 
the  x/c  and  y/c  coordinates  are  listed  in  separate  blocks,  with  the  order  starting  at  or  near 
the  trailing  edge,  proceeding  across  the  lower  surface,  then  the  upper  surface,  and  ending 
at  or  near  the  trailing  edge.  The  number  of  points  will  be  one  higher  than  the  number 
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of  panels,  even  if  the  first  and  last  points  coincide.  The  trailing  edge  point  simply  is 
listed  twice  in  this  case.  The  only  part  of  this  file  that  is  likely  to  get  changed  on  a 
routine  basis  is  the  angle  of  attack,  as  the  other  values  are  usually  fixed  for  a  given 
airfoil. 

The  other  input  file  pertains  to  the  boundary  layer  part  of  the  code.  This  file 
is  normally  called  incompbl.dat,  although  the  user  may  give  it  any  name  and  use  this 
name  on  the  command  line  at  the  time  of  execution.  A  sample  input  file  is  included  in 
Appendix  A.  IWAKE  is  the  viscous  wake  flow  flag.  A  zero  indicates  that  these  effects 
will  not  be  included,  while  a  one  indicates  that  they  will  be  included.  NXT  dictates  the 
number  of  chordwise  stations  on  the  body.  NW  sets  the  number  of  chordwise  stations 
in  the  wake.  ITREND  refers  to  the  number  of  calculation  cycles,  where  20  is  a  good 
starting  number,  and  30  or  40  may  be  needed.  ITR(l)  is  a  flag  for  the  transition  location 
specification  method  for  the  upper  surface.  A  zero  will  activate  a  calculation  using 
Michel’s  method,  and  a  one  indicates  that  the  location  will  be  provided  as  part  of  the 
input.  ITR(2),  for  the  lower  surface,  should  be  zero.  ISWPMX  is  the  number  of  sweeps 
in  each  cycle.  A  cycle  corresponds  to  the  calculation  of  inviscid  and  viscous  flow 
equations.  One  sweep  is  usually  sufficient  but,  in  some  cases,  it  may  be  necessary  to  use 
2  or  3  sweeps  in  one  cycle.  RL  is  the  Reynolds  number  based  on  chord  length. 
XCTR(l)  is  the  x/c  value  for  the  transition  location  on  the  upper  surface.  This  value  is 
only  used  if  rTR(l)  =  1;  otherwise,  it  will  be  ignored.  IP  is  the  print  flag,  which  should 
normally  be  set  to  one  to  obtain  output.  This  screen  output  can  be  redirected  to  a  file 
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for  later  use  by  including  the  proper  option  on  the  command  line  at  the  time  of 
execution.  If  IP=0,  the  standard  output  will  not  be  generated. 

5.  Program  Execution 

After  ensuring  the  appropriate  files  are  properly  edited  and  are  present  in  the 
directory,  type: 

>incompbl  <  incompbl.dat  >  incompbi.out 
Several  new  files  will  appear  in  the  directory  after  running  the  program.  Any  write 
statements  incorporated  into  the  program  using  the  default  numbering  format  will  produce 
files  such  as  FOR032.DAT  (Stardent)  for  fort.32  (IRIS).  The  output  file  incompbi.out 
will  have  a  comprehensive  summary  of  the  input  data,  and  output  data  such  as  Cp,  Co, 
Cl,  Cm,  and  boundary  layer  properties,  including  skin  friction  coefficient  and 
displacement  thickness. 
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V.  EXPERIMENT 


A.  BACKGROUND 

The  experimental  work  described  in  this  chapter  was  performed  as  part  of  a 
continuing  series  of  investigations  by  M.  S.  Chandrasekhara  and  L.  W.  Carr  in  the 
Compressible  Dynamic  Stall  Facility  (CDSF)  of  the  Fluid  Mechanics  Laboratory  (FML) 
at  NASA  Ames  Research  Center.  The  steady  flow  density  field  around  a  NACA  0012 
airfoil  at  a  Reynolds  number  of  540,000  was  photographed  using  the  technique  of  Point 
Diffraction  Interferometry  [Ref.  7]. 

A  sample  photograph,  called  an  interferogram,  is  shown  in  Figure  5.1.  The  bright 
and  dark  areas  emanating  from  the  airfoil  are  called  fringes,  and  each  one  represents  a 
line  of  constant  density.  The  stagnation  point  may  be  easily  identified  as  the  center  of 
the  smallest  fringe  on  the  lower  surface  (for  a  positive  angle  of  attack).  The  flow 
accelerates  around  the  leading  edge  of  the  airfoil.  The  example  shown  also  reveals  the 
presence  of  a  laminar  separation  bubble  just  aft  of  the  leading  edge,  distinguished  by  a 
characteristic  fringe  pattern. 

For  a  given  freestream  Mach  number,  the  fringes  may  also  be  correlated  with 
particular  Mach  numbers  and  pressures  in  isentropic  flow.  For  a  standard  interferometer 
and  two  dimensional  flow,  the  path  length  difference  APL  due  to  density  changes  can  be 
related  to  the  fringe  number  e; 
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where  n  is  the  refractive  index  of  the  signal  beam,  n„f  refers  to  the  reference  beam,  \ 
is  the  wave  length  of  light  used,  and  L  is  the  test  section  span.  An  integer  value  for  e 
results  in  a  bright  fringe,  while  an  integer  plus  one  half  corresponds  to  a  dark  fringe. 
Using  the  Gladstone-Dale  equation  [Ref.  8]  and  the  perfect  gas  equation.  Equation  5.1 
reduces  to 


p-p-=(^)(^)Me  (5.2) 

/Iq-I  L 

where  p  is  the  density  of  the  fringe,  p«  refers  to  freestream  conditions,  and  A  is  a 
constant  determined  from  the  experimental  parameters.  With  Xo=532  nm,  L=25  cm, 
(no-l)=2.733  X  lO"^  ,  and  total  or  stagnation  density  po=l,21  kg/m^  A=0.(X)9421 
kg/m^  Dividing  by  Po, 


Using  the  relationship 


p  _  P»^i4e 

Po  Po  Po 


1 

»-T 


(5.3) 


(5.4) 


withp=Pc„  and  M=M*,  pJpo  is  a  function  of  freestream  Mach  number  only  for  7=  1.4. 
For  the  present  case  of  M=0.3,  the  term  p*/po  is  constant  at  0.956.  Equation  5.3  may 
now  be  written  as 
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=0.956  +  0.007786€ 
Po 


(5.5) 


Thus,  quantitative  density  measurements  are  available  knowing  only  the  fringe  number. 
The  fringe  numbers  are  determined  by  identifying  the  fringe  number  of  the  stagnation 
fringe  as  the  highest  fringe  number,  and  counting  down  around  the  leading  edge  and 
upper  surface  of  the  airfoil.  With  p/po=l  Equation  5.5,  e  *6.  Knowing  the  densities, 
pressures  may  be  calculated  in  a  straightforward  manner: 


(5.6) 


Mach  numbers  for  each  fringe  may  be  calculated  using  Equation  5.4. 

The  photographs  for  various  angles  of  attack  are  identified  using  particular  settings 
on  the  equipment.  A  correlation  between  settings  and  angles  of  attack  is  presented  in 
Table  5.1. 


Table  5.1  CORRELATION  OF  EXPERIMENTAL  SETTINGS  AND  AOA 
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154 

4 
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9.5 
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57 

10 
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B.  IMAGE  PROCESSING 


1,  Scanning 

The  processing  task  began  with  a  set  of  experimental  interferograms  at  angles 
of  attack  ranging  from  0®  to  10®.  Each  was  scanned  into  an  IRIS  computer  using  the 
program  pixscan  at  the  Numerical  Aerodynamic  Simulation  Facility  at  NASA  Ames 
Research  Center.  Options  for  contrast  enhancement  (gamma  correction)  and  grayscale 
were  used.  A  photograph  of  the  airfoil  in  no-flow  conditions  was  also  scanned.  Ail  of 
the  photographs  were  carefully  placed  on  the  scanner  in  a  position  such  that  the  top 
surfaces  of  the  two  bottom  reference  triangles  made  a  horizontal  line  even  with  the  edge 
of  the  scanning  window.  The  resulting  scannal  images  were  transferred  via  ftp  to  the 
IRIS  computer  at  the  Naval  Postgraduate  School  for  further  processing. 

2.  Editing 

Using  a  program  called  pixedit,  the  original  flowfield  images  were  first  overlaid 
with  the  airfoil  image.  This  permitted  a  more  defined  surface  for  analysis  when  the 
diffraction  caused  by  the  interferometry  technique  distorted  the  airfoil  surface.  The 
images  were  then  cropped  closer  to  the  airfoil.  This  procedure  reduced  the  number  of 
pixels  in  the  image,  which  was  necessary  for  use  with  the  digitizing  program.  A  small 
section  near  the  leading  edge  was  also  selected  for  enlargement,  thus  giving  better  detail 
of  the  laminar  separation  bubble  region. 
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In  anticipation  of  the  digitization  process,  reference  marks  at  known  coordinates 
were  added  to  the  two  images.  The  original  photographs  provided  the  basis  for 
coordinate  system  definition.  The  three  dark  triangles  are  located  as  shown  in  Figure 
5.2.  The  distance  between  horizontal  or  vertical  faces  of  the  triangles  is  0.2  of  the 
chord.  The  vertical  faces  of  the  two  triangles  on 
the  right  are  located  at  0.25  chord.  Thus,  the 
coordinates  of  the  point  at  the  right  angle  of  the 
bottom  right  triangle  are  (0.25,  -0. 1).  A 
reference  mark  is  needed  in  the  top  left  corner 
for  use  with  the  digitization  program,  but  none  is 
provided.  Therefore,  a  mark  was  constructed 
with  several  applications  of  the  snap  new  image 
feature  of  pixedit.  Very  thin  horizontal  and  vertical  lines  were  initially  saved  as  separate 
images.  For  each  interferogram,  these  images  were  read  in  and  placed  to  be  even  with 
existing  reference  marks  of  the  airfoil  leading  edge  (vertical  line)  and  the  bottom  edge 
of  the  upper  triangle  (horizontal  line).  The  small  area  in  the  top  left  corner  where  the 
lines  crossed  was  saved  as  another  image  and  the  long  lines  were  deleted.  The 
coordinates  of  the  cross  mark  are  (0.0,  0.1).  Finally,  the  complete  image  was  saved  for 
further  processing.  A  similar  procedure  was  used  to  make  a  cross  mark  at  (0.05,  0.0) 
for  the  lower  right  reference  on  the  enlarged  leading  edge  images. 


Figure  5.2  Reference  Triangles 
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3.  Fringe  Tracing 

Both  of  the  new  images,  the  density  flowfield  and  the  enlarged  leading  edge, 
were  digitized  using  the  program  DigiCurv.  The  appropriate  lower  right  and  upper  left 
coordinates  were  provided  and  the  corresponding  points  on  the  image  were  selected  with 
the  mouse.  With  the  coordinate  system  thus  initialized,  each  curve  was  digitized 
separately.  The  centerline  of  the  dark  fringe  was  used  for  digitization.  A  representative 
number  of  points  were  chosen  for  each  fringe  using  the  left  button  of  ^he  mouse,  always 
starting  with  the  point  where  the  fringe  intersected  the  airfoil.  Depressing  the  right 
button  activated  a  menu  with  a  curve  fit  option.  The  program  calculated  a  best  fit  curve 
to  the  chosen  points  using  a  spline  routine  and  output  up  to  30  (default  value)  new  points 
to  describe  the  curve.  The  fitted  curve  was  also  displayed  on  the  screen  for  acceptance. 
In  most  cases,  the  computed  curve  had  outstanding  agreement  with  the  centerline  of  the 
dark  fringe,  even  when  a  sharp  comer  was  involved  near  the  end  of  a  separation  bubble. 
If  the  displayed  curve  needed  refinement,  the  backup  option  removed  points  one  by  one, 
more  closely  spaced  pwints  were  selected,  and  the  curve  was  refit.  After  all  dark  fringes 
were  digitized  for  an  image,  the  data  point  output  file  was  saved  for  plotting. 

The  digitized  data  can  be  viewed  using  any  plotting  routine  compatible  with 
sequential  data  separated  by  blank  lines.  The  program  xyplot  does  not  recognize  blank 
lines  and  is  not  a  good  choice.  The  program  gnuplot  will  properly  accept  data  in  this 
format.  The  digitized  plot  corresponding  to  the  interferogram  in  Figure  5. 1  is  shown  in 
Figure  5.3. 


86 


4.  PriatiAf 

To  ot>Uin  graphical  outpui  on  a  pouscnpt  Uucr  pnnurr,  the  files  must  be  in 
postacnpt  formal,  fhe  dtgitued  plots  are  in  this  format  and  can  be  easily  printed  using 
the  following  commands  in  the  Advanced  Compuiauon  L^ioratory  at  the  Naval 
PoMgraduate  School.  From  the  Siardent.  type 

>  rglp  JUtHomt 

From  the  IRIS,  type 

>  Ip  -dlaser  fiitmamt 

The  filenames  may  have  the  suffu  ps  to  indicate  that  they  are  postscript  files,  but  this 
IS  for  reference  only  and  is  not  required  by  the  printer 

The  images  outpui  from  pixrdit  are  not  directly  compatible  Nsith  the  posiscnpt 
printer.  Using  the  IRIS  computer,  they  must  first  be  converted  to  sgi  format: 

>  pixlsgi  inJUiiuumi  outfUeiumi.igi 

Again,  the  suffix  is  optional.  The  next  step  is  a  conversion  to  postscnpt  format: 

>  tops  infile namt. sgi  -p  98.0  >  outfiiename.ps 

All  suffixes  are  optional.  The  -p  98.0  option  preserves  the  proponions  of  the  original 
image;  otherwise,  it  may  be  distoned  when  stretched  to  fill  up  the  printed  page. 
The  >  symbol  redirects  the  screen  output  to  the  specified  output  file.  There  are  a 
number  of  other  options  available,  including  size  and  color.  For  a  full  description  of 
options,  type 

>  man  tops 

to  obtain  the  on  line  help  manual  entry  on  this  conversion  program. 
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C.  ANALYSIS 


TSc  pnnury  utfonnaiid*!  drufed  from  ihe  ii»ieftcro<nim»  w  as  t*!  the 

tfart  of  ihc  Uirunar  icfittmion  buti6k.  ii  raitfod  Fi^urr  5  4  «ii  mlaffod 
imafc  41  6'‘ .  vkhcfc  a  bubble  it  futi  formed  Uttnf  the  d«(^ttired  data  ftto  of  thu  ima^e. 
pfcctic  a/c  locaticwit  of  ihc  interwctiOAt  of  the  fnAftt  «kiih  the  airfoil  wrfacic  mat  be 
determined,  at  it  (he  firti  potni  difitired  for  each  fnnfe  Sumn^  laiih  ihc  tta^naiion 
fnnge,  ihe  firti  dark  fnnge  aiviayt  hat  a  number  of  >  5  (for  (he  pretem  cate  of 
M,  "OJ)  Subtequeni  dark  innget  have  numbv  4  5,  5.  etc  Negative  fnnge 
numbert  are  pottible.  Converting  pretturet  to  pretturc  coefficienn  (C^).  a  table  of 
fringe  numbert  and  corretponding  C/t  may  be  produced  and  utod  for  all  interferogramt 
with  the  tame  freetiream  Mach  number.  The  program  pres  (Appendix  A)  wat  wntteti 
to  process  the  data  as  described  and  output  a  plotting  file  of  x/c  vs.  C^.  Figure  S.S 
shows  the  digitized  image  at  6*  and  the  corresponding  C,  plot.  The  most  important 
feature  is  the  pressure  plateau,  which  reveals  the  existence  and  location  of  a  separation 
bubble.  The  bubble  starts  at  a  chord  location  of  O.OIS.  Interferograms  for  other  angles 
of  attack  were  analyzed  in  the  same  manner.  Figure  5.6  shows  a  comprehensive  plot  for 
an  angle  of  attack  range  from  6°  to  10°.  As  the  angle  of  attack  increases,  the  length  of 
the  bubble  increases,  and  the  starting  location  moves  slightly  forward. 


89 


Figure  5.4  Enlarged  Leading  Edge  Interferogram,  AOA=6° 
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Rgure  5.5  DigHized  Interferogram  and  Plo<,  AOA=6* 
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D.  COMPARISON  OF  RESULTS  TO  COMPUTATION 

Table  5.2  compares  the  experimental  results  to  those  obtained  by  bt2d  and 
incompbl.  The  direct  boundary  layer  code  shows  a  very  small  separation  bubble  at  4°; 
however,  none  is  present  in  the  experiment.  At  6®,  the  bubble’s  existence  is  correctly 
detected,  but  the  computational  location  is  slightly  aft  of  the  experimental  value.  The 
higher  angles  of  attack  show  excellent  agreement,  with  differences  of  only  0.3%  of 
chord.  In  contrast,  the  viscous-inviscid  code  failed  to  predict  any  separation  until  13.3°. 
Since  this  is  past  the  steady  stall  angle  of  12.4°,  there  is  no  experimental  bubble  to 
compare  with.  Even  if  stall  had  not  occurred,  the  trend  clearly  indicates  that  the  location 
would  be  in  great  error  as  well. 


Table  5.2  COMPARISON  OF  BUBBLE  START  LOCATIONS 


1  AOA  1  BL2D 

INCOMPBL 

EXPERIMENT 

0 

... 

— 

... 

2 

— 

— 

— 

4 

.245 

— 

... 

6 

.038 

— 

.015 

8 

.017 

— 

.014 

10 

.012 

— 

.009 

13.3 

NO  SOLN 

.024 

STALL 
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VI.  NAVIER-STOKES  CODE 


A.  OVERVIEW 

A  time-averaged  Navier-Stokes  (NS)  code  called  iis2  was  used  for  a  final  study. 
This  method  has  the  advantage  of  including  compressibility  effects.  White  at  a  Mach 
number  of  0.3  the  effects  are  small,  it  is  just  on  the  border  of  the  region  that  is  usually 
considered  acceptable  for  the  assumption  of  incompressibility.  A  disadvantage  of  this 
method,  as  mentioned  in  the  transition  discussion  of  Chapter  2.  is  the  lack  of  a  transition 
model.  Turbulent  flow  is  assumed  throughout  the  flowfield.  The  code  also  takes  three 
to  four  hours  to  run  on  a  Stardent  computer.  Nevertheless,  Navier  Stokes  codes  are 
often  regarded  as  the  best  compulation  method  currently  available  and  warrant 
consideration.  The  details  of  the  lime-averaged  NS  equations,  their  derivation,  and  their 
discretization  are  well  documented  elsewhere  fRef.  9J  and  will  not  be  reviewed  here. 
The  particular  code  used  for  this  investigation  was  developed  by  J.  A.  Ekaterinaris  of  the 
Navy-NASA  Joint  Institute  of  Aeronautics. 

B.  RESULTS 

The  code  was  run  for  a  NACA  0012  airfoil,  with  a  Mach  number  of  0.3  and  a 
Reynolds  number  of  540,000.  A  161x64  C-type  grid  with  a  very  fine  distribution  normal 
to  the  surface  in  the  viscous  region  was  used.  The  grid  is  shown  in  Figure  6. 1 .  The 
Baldwin-Lomax  model  was  used  for  turbulence  modeling.  The  program  was  run  to  20(X) 
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Figure  6.1  161x64  Viscous  Grid  for  NACA  0012  Airfoil 
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iterations,  where  the  residuals  had  dropped  two  orders  of  magnitude  Fii^uto  6  T  k<  e  4 
show  comparisons  of  the  computational  density  fields  to  the  etpcnmenul 
At  0°,  there  is  no  separation  bubble  and  the  afreemeni  appears  to  tw  a;  n'  tv 

experimental  bubble  first  appears  at  x/c  »0  015  Ail  of  the  computaitonal  omxify  c  unrs 
emanate  from  the  leading  edge  This  is  consisicni  with  the  lack  oi  «  iranMUon  mudri 
At  10°,  both  methods  show  a  targe  bubble  The  SS  bi^ible  is  muc!.  tunher  ait  there 
is  also  a  small  extra  bubble  near  the  leading  edge  It  appears  that  the  *s'u.te  erf  the  art 
method  is  not  always  best  for  represcntmi  reality  pafixuUr'y  k»«  ttr>m<«;<Ss  ngn.hrt 
flows.  The  author  is.  howeser.  currently  »o«kiftg  on  a  ^ervon  oi  the  pfe^gratn 
will  incorporate  transition  cakuUtions  and  preliminary  sMSa^asrs  that  rrsui'.i^ 


be  much  better 
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Figure  6.3  Comparison  of  NS  and  Experiment,  AOA=6“ 


Figure  i.4  Con^Mfisoii  ot  MS  siid  F^perimml.  AOA  » 10” 


vn.  CONCLUSIONS 


In  the  computational  investigations  of  a  NACA  0012  airfoil  in  low  Reynolds 
number  flows,  several  important  discoveries  have  been  made.  First,  even  though  Navier- 
Stokes  codes  are  the  most  advanced  computational  method  currently  available,  they  are 
not  always  appropriate.  Specifically,  a  transition  model  is  necessary  for  the  case  of  low 
Reynolds  number  flows.  Even  if  accurate  results  may  be  obtained,  the  method  is  not 
practical  for  many  applications  due  to  its  high  cost  in  time  and  money.  The  viscous- 
inviscid  method  seems  to  offer  very  advanced  calculations  at  a  very  inexpensive  price. 
The  problem  with  this  code  used  alone,  however,  is  that  it  simply  does  not  give  correct 
results  for  separation  bubbles  in  low  Reynolds  number  flows.  A  compressibility 
correction  in  a  future  version  may  alleviate  the  problem.  Using  a  stability/transition 
method  in  conjunction  with  the  code  will  certainly  provide  more  refined  transition 
estimates,  which  may  influence  separation  bubble  results.  Finally,  the  relatively  simple 
and  often  overlooked  direct  boundary  layer  method  can  provide  meaningful  information 
about  separation  bubbles  in  low  Reynolds  number  flows.  Given  that  the  code  is  also  very 
fast,  it  may  be  used  efficiently  in  the  design  stages  and  quality  assurance  checking  of 
many  aeronautical  af^licadons. 

Successful  or  not,  all  current  computational  methods  have  several  important 
limitations.  All  are  dependent  on  empirical  models  for  transition  onset,  transition  length, 
and  turbulence.  The  models  are  often  formulated  for  specific  parameter  ranges  outside 
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of  which  agreement  is  poor.  Until  there  is  sufficient  computing  power  available  to  solve 
the  full  Navier-Stokes  codes,  it  is  imperative  to  check  the  applicability  of  the  models  and 
the  assumptions  of  the  equations  on  which  a  method  is  based  before  using  a  program  and 
counting  on  the  results  to  be  reliable. 


101 


LIST  OF  REFERENCES 


1.  Eppler,  R.,  Airfoil  Design  and  Data,  Springer- Verlag,  1990. 

2.  Anderson,  J.  D.,  Fundamentals  of  Aerodynamics,  2d  ed.,  McGraw-Hill,  Inc., 
1991. 

3.  Cebeci,  T.,  and  Bradshaw,  P.,  Momeruum  Transfer  in  Boundary  Layers, 
Hemisphere  Publishing  Corporation,  1977. 

4.  Neace,  K.  S.,  A  Computational  and  Experimental  Investigation  of  the  Propulsive 
and  Lifting  Characteristics  of  Oscillating  Airfoils  and  Airfoil  Combinations  in 
Incompressible  Flow,  Engineer’s  Thesis,  Naval  Postgraduate  School,  Monterey, 
CA,  September  1992. 

5.  Walker,  G.  J.,  Subroto,  P.  H.,  and  Platzer,  M.  F.,  "Transition  Modeling  Effects 
on  Viscous/Inviscid  Interaction  Analysis  of  Low  Reynolds  Number  Airfoil  Flows 
Involving  Laminar  Separation  Bubbles,"  ASME  Paper  No.  88-GT-32  presented  at 
the  Gas  Turbine  and  Aeroengine  Congress,  Amsterdam,  The  Netherlands,  6-9  June, 
1988. 

6.  California  State  University  Aerospace  Engineering  Department  Report  AE-90-2, 
An  Interactive  Boundary-Layer  Stability-Transition  Method  for  Calculating  Flow 
Over  Hydrofoils,  by  T.  Cebwi,  H.  H.  Chen,  and  H.  M.  Jang,  August  1990. 

7.  Carr,  L.  W.,  Chandrasekhara,  M.  S.,  Brock,  N.  J.,  and  Ahmed,  S.,  "A  Study  of 
Dynamic  Stall  Using  Real-Time  Interferometry,"  AIAA  Paper  No.  91-0007 
presented  at  the  29*  Aerospace  Sciences  Meeting,  Reno,  NV,  7-10  January,  1991. 

8.  Goldstein,  R.  J.,  Fluid  Mechanics  Measurements,  Hemisphere  Publishing 
Corporation,  1983. 

9.  Anderson,  D.  A.,  Tannehill,  J.  C.,  and  Pletcher,  R.  H.,  Computational  Fluid 
Mechanics  and  Heat  Transfer,  Hemisphere  Publishing  Corporation,  1984. 


102 


APPENDIX  A:  COMPUTER  PROGRAMS  AND  FILES 


The  following  programs  and  input/output  files  are  presented  in  logical  order  of  use.  The 
input  and  output  files  are  examples  only.  The  input  values  may  be  modified  as  described 
in  the  appropriate  User’s  Guide  section.  The  program  incompbl  comprises  90  pages  in 
length  and  was  not  modified  other  than  inserting  WRITE  statements  to  extract  data  for 
plotting;  thus,  the  source  code  is  not  included. 
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Input  file:  points.dat 


: .000000 

o.ooocco 

C. 999013 

"0.000141 

0.996057 

-0.000562 

0.991144 

-0.001258 

0.984292 

-0.002222 

0.975528 

-0.003443 

0.964888 

-0.004909 

0.952414 

-0.006603 

0.938153 

-0.008510 

0.922U4 

-0.010610 

C. 904508 

-0.012883 

0.88525? 

-0.015310 

0.864484 

-0,017868 

0.842274 

-0.020535 

C. 818712 

-0.023291 

0.793893 

-0.026111 

0.767913 

-0.028974 

0.74C877 

-0.031856 

0.712890 

-C. 034733 

0.684062 

-0.037582 

0.654508 

-0.04C378 

0.624345 

-0.043094 

0.593691 

-0.045705 

0.562667 

-0.048182 

0.531395 

-0.C5C499 

0.500000 

-0.052625 

0.468605 

-0.054534 

C.  437333 

-0,056195 

C. 406309 

-0.057581 

0.375655 

-0.C58666 

C. 345491 

-0.C59424 

0.315938 

-0.059834 

0.287110 

-C.C598'’6 

0.259123 

-0.059535 

0.23208? 

-0,058799 

0.206107 

-C. 057661 

0,191288 

-0,056119 

C.  1  57726 

-0,054176 

0,  :  35516 

-0.051039 

0.114743 

-0.049121 

C.C9S491 

-C.C46037 

C. 077836 

-0,042608 

0,061847 

-0,038854 

0,047586 

-0.034800 

0.035112 

-0,030471 

C.C24472 

-0.025093 

0.015708 

-0.02:008 

O.OCS356 

-0.016078 

C.0C3943 

-C. 010884 

c.ooooe’ 

-0,005521 

C  .OOCCCC 

c.oococo 

0.000987 

C.0C5521 

C.OC3943 

C. 010804 

C.CO0856 

O.OlfC’B 

C.C157C8 

0.021088 

C. 024472 

C. 025893 

C. 0351  12 

0,030471 

0.047586 

0,034800 

C.06184-» 

0.038854 

0.077038 

0.042608 

C. 295491 

0.046037 

0.114743 

0.049121 

C. 1 35516 

C. 051839 

C.  1  5  7  7;4 

0.0541  :6 

C.  If  1218 

0.056119 

0.2061C7 

C. 05766* 

C. 23208  7 

0.058 ’99 

C. 259123 

0.059535 

c.2f ic 

0.C59876 

c  315938 

0.059834 

r.  M54«2 

^ .  05'»424 

C.3’5655 

0.058666 

0. 406309 

Oi  AS  '  C  0  2 

c!43-333 

c]c56i95 

C.  468605 

0.054:34 

C.5CC7CC 

0.052625 

0.531395 

0.050499 

C. 562667 

0.048:82 

0.59369! 

0.:45-'C5 

0.624345 

0.043094 

0.65450? 

0.040370 

0.684C62 

0.037582 

0.712890 

0.034733 

0.  ’4C8’7 

0.031856 

C .767?: 3 

0 .:?9974 

0. 7?3B?3 

0.026111 

0.8»8M2 

0.02379! 

0.8422:4 

C.02C535 

0.864484 

0,01 '968 

C.  88525  : 

C.01531C 

0.904509 

3. 012003 

0.922164 

C. 010610 

0.938153 

0.008510 

0.952414 

C. 006603 

0. 964ff? 

0.CC4909 

0. 975529 

C. 003443 

C.984;*: 

0.032222 

:  991:44 

0.001758 

0 . 994'‘5  • 

0.^0:562 

^  f***:  •  > 

: . :  c :  1 4 1 

1 . 

»• 

Source  code:  panel.f 


Page  i 


Fpyr.aivM  fanei 

AL'T'C'R;  L.  M.  NOMAK 

da:?.;  t  NOV  i\  (nodutadi  ?1AV,  Aur, 

FUAfUEF.!  CAlCUDATf  THE  VELOCITIEJ  ON  AM  AIRrojI.  OaJN'J  A  PAfif.L  METHOD. 

t  !M!  Array*  currently  dtn>«>n*lon*d  Jot  maKlnaff  at  N»JC0  pinel* 

'nDJ?  data  fllr  Falnta.dat  will  have  N*1  points 

Ojtpur  valocltlaa  ara  referenced  to  freest rese,  Je.  V/Vlnf 

MEtHE”'!  ELCWriELO  CONSISTS  Of  THREE  SIMPLER  fLOMSs  TREESTREA-H,  SOURCE, 
AND  VOBTICITY.  SOURCE  DISTRIBUTIONS  qljl  VARY  FROM  PANEL  TO 
PANEL.  VORTICITY  STRENGTH  GAMMA  IS  THE  SAME  FOR  ALL  PANELS. 
boundary  CONDITIONS  INCLUDE  FLOW  TAMCF.NCY  AT  CONTROL  POINTS  AND 
KUTTA  CONDITION  FOR  FIRST  AND  LAST  PANELS.  INFLUENCE 
COEFFICIENTS  COMBINED  TO  FORM  NEW  COEFFICIENTS  IN  LINEAR  SYSTEM 
OF  n*!  E8JATIONS,  n»l  UNKNOWNS  (q ( 1 I . . .q ( n) ,  GAMMA).  VELOCITIES 
AT  CGNTROI.  POINTS  CVALUATCO  FROM  ql)|  ANO  OAMMA. 

REA!  xn:?02),Y(ll202),XM{l;20?).YM)l:20J), 

:  At (1 : 202, l:202),Bt 11:202, 1:202), 

:  a(:!2C2, 1:202), b(l:202), 

:  q(l :202) ,Vt  <1 :2C2) ,ALPHA,V,VtC  11:202) , 

:  P1,GAMMA,TIIETAI1:202),NUM.DEN, 

:  P  d  :202,  1  :2C2) ,  BETA (1 : 202, 1 ; 202) ,  NUMl ,  CENl,  NUM2,OEN2, 

;  AAi'Gl!  :  202,  1  :202),An(  1:202, 1:202 ),  Bn  II  :202, 1  :202) 


•  HUMBER  or  Noors  ON  AIRFOI!.  SURFACE: 

•  PRINT",  '  INPUT  NO.  OF  PANELS  (1  leaf  t)ian  lllnej  In  oolnt! 

•  READ  ".N 

FI«AC0S(-: .) 

OPEN  lUNIT-eS.FlIE  '  rt  1 '*.  » .dat ' ,  STATUS* '  "JNENOW*;' ) 

OPEN  (UN IT* 8?,  FILE-'  vo .-a: .dat STATUS* ' UNKNOWN' ) 

OPEN  (UN:t-91 , FILF ’*  cp.dat ' , STATUS*'  UNKNOWN’ ) 

OPEN  (UNIT  «C,FtLP-'bl2d.dat',STATUS*’UNKNOWN') 

print  •,'TNPUI  REYNOtnS  NUMBER:' 

REAP  ‘.RI. 

print  ','EN*rP  r  IP  TRANSITION  LOCATIONS  ’USKNOWN’ 

PRINT  ','  ’  TOA'iriTION  LOCATIONS  KNO)i'N:' 

READ  '.IAN: 

IF  (lANS.EC. ’. )  THEN 

PRINT  INPUT  X/C  TRANSITION  LCCATION  FOR  UPPER  S’.'RFAC 
READ  *  TRANSUPPER 

PRINT  INPUT  X/C  TRANSITION  LOCATION  FOR  tCKEP  SURFAC 
READ  TRANSLOWER 
ELSE 

•••These  are  arbitrary  values  Intended  to  be  dewnstresn-  c!  the 
•••  actual  trar.siticr  points,  for  use  with  Michel’s  criterion  ; 
TPANSUrFER*.B 
TRANSLOWER*.  99® 

ENDIF 

WRITE  (90,50)  R!,,TP.ANS'JPPER,TFANSICWER 
WRITE  199,501  RL, TRANSUPPER, TRANSLOWER 
50  FORMAT  (F10.0,ciC.<,  FIO.91 

PRINT  •,' INPUT  ANGLE  Of  ATTACK  IN  DEGREES:’ 

READ  •, ALPHA 
ALPIIA»Al.rHA*PI/l  BO.O 


in  BL2D 


00  30  1*1, N*1 

READ  (88,251  X(1),Y(II 
rCBMAT  (P(r8.«,2*)! 
CONTINUE 
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Paae2 

u 

u 

li 

•T*>1»  «»'CT!on  dn!!n«»  th*  influencd  co«f  f  Iclent  s: 

DO  Iil,N 

H 

XM(1)'0,5MX(I).X(I«1H 

n 

Kf<(n*0.5*(T(ll*lf(l*lM 

J8 

8(1, 11  •  tXMdl  -X(l)  I  ‘n.*  1  ••2. 

7» 

DO  J«!,N 

80 

NUM”Y(Jdl-Y(J) 

81 

0EN*X(J.1)-X(JI 

82 

THr.TX  ( Jl  •  ATAN2  (NUM,  OENJ 

8) 

N0K!-YM(I)-Y(J*1) 

84 

DENl<XM(ll*X(JM) 

88 

NUM2»YM(n'Y(JI 

8« 

DEN2*XH(1)-X(JI 

87 

BETAd,  Jl  •ATAN2  (  (MUMl  ‘OEW-OENl  •HUH2)  ,  (DENI  •DrN2-N'iW,  I 

88 

R(l,J*lt«(XM(I)-X(J*lll‘*2.»(YMdl-Y(j*ll)**2. 

88 

THETADir*THETA(II -THETA(J) 

80 

ir  (l.EO.JI 

81 

;  THEN 

82 

And,  J)  *0.5 

83 

Bnd,  Jl*0.0 

84 

ELSE 

88 

An  (1,  Jl  •  (1/(2‘P1I  1  •  (SINdHEIADin  *ALCC  (P  ( I .  J‘ 1 !  /P  (1,  .') ) 

88 

:  •  .5*005 (THETAOin ‘BETAd,  Jl  1 

87 

Bn  ( 1 .  Jl  •  ( 1  /  (2  * P 1 1 1  *  (COS  (THETADID  •  ALOG  ( 8  (1 ,  J*  1 1  /  P  ( * ,  Jl  1 

88 

;  •  .5-SIN(TMETAOiri ‘BETAd,  J)  1 

88 

ESC  ir 

100 

At  d,  Jl--8r>.  d,  Jl 

101 

8td,JI-And,JI 

102 

END 

103 

END  DC 

1C4 

105 

•  5'n*.rlx  coefficients  llnenr  tys'er  defined  Is'i  »nd  b'slt 

108 

107 

•  (N‘1,H‘11«0.0 

108 

00  :»1,N 

108 

•  d,N‘ll*C.0 

no 

00  J-i'.,N 

111 

•  (l,Jl«And,J) 

;!7 

*(l,N‘ll*»(l,N‘ll*Snd,Jl 

113 

ENO  DC 

114 

bdl  ’-1 .0*SlN(ALPHA-TMETA(in 

115 

*  (N‘l,  11  ‘At  d,  11  ‘At  (N,  :i 

lit 

*  (K‘1,N*11  "a  (Nd  ,N*1I  ‘Bt  (1, 1)  *Bt  (N,  11 

117 

END  DO 

111 

b(N*ll •-! ,0* (CCS (ALPHA-THETA (1 1 1 -COS (ALPHA-THETA (Nl 1 1 

118 

,  , 

120 

•  Define  augmented  mattl*  fcr  Input  to  linear  solver  lubroj'Ine  GA'Jsr- 

121 

122 

DO  1*1, N*1 

123 

DC  J»1,N*: 

124 

AA'v'Cd,JI  T8(l,  Jl 

125 

END  DO 

128 

AAUad,N<2l  'bdl 

127 

END  DO 

128 

128 

CALI.  GAUSS  (Nd.AAUCI 

no 

131 

•  Define  source  snd  vortlcity  strengths: 

132 

133 

DO  I*1,N 

134 

qdl  •AA'JC(I,N«2l 

135 

END  DO 

138 

137 

CAMH(A'AAUC(N‘l,N‘2t 

138 

138 

*  Calculate  velocity  on  each  panel  at  control  point 

140 

141 

NSTAGELAC^C 

142 

JSTAG-'O 
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:  13 
HI 
Hi 
H( 
HI 
HS 
H» 

110 

111 
is: 
1S3 
151 
ISS 
1S« 
1S1 

158 

159 
HO 
HI 
H2 
H3 
HI 
H5 
H€ 

H8 
H9 
1  7? 
PI 
1  72 
P3 
P4 
PS 
!  74 
;  77 
P8 
P9 
18C 
181 
182 

193 

194 
las 
!S« 

187 

188 
189 
3  9C 
1  91 

192 

193 
:  94 

195 

194 

197 

198 

199 

200 
201 
202 

203 

204 

205 
204 

207 

208 

209 

210 
211 
212 
213 


nti  1 » 1 ,  N 
VI (11 '0.0 
00  JM,N 

Vt  (II  'At  (I,  J|  'qlJl  'CAWIAHt  (I,  J1  'Vt  (II 
END  DO 

VI m • Vt ( I I 'COS (ALPHA-THETA (III 
Cp-l.O-Vt (II ••2 

<•  WRITE  (25,  *1  KMdI.Vtdl 

ir  ((VtdI.GT.Ol  .AND.  (NSTACrLAG.EO.Ol  1  THEN 
1  STAG' I 
NSTACrLAGd 
ENDir 

IF  (VtdI.LT.OI  Vtdl--Vtdl 
WRITE  (89,  451  XM(ll,Vtdl 
WRITE  (91,451  XKdl,-Cp 
END  DO 

45  FORMAT  (2 (FI 0 . 5, 2X1 1 

48  FORMAT  (3(110.511 

49  FORMAT  (3151 

WRITE  (90,491  N, ISTAG, IANS 
DO  I -UN 

WRITE  190,  481  XMdl  ,  YKdl  ,  Vt  dl 
END  CO 

print  •, 'CAt.CULAT! CHS  COMPLETE' 

PRINT  ♦, 'OUTPUT  FILES  ARE  v«l  o- .  dnt ,  cp.dat,  H/d.flat' 
END 


•  Gauss  e!  Imlr.at  (-'  procedure  obtained  from  Numerical  Methods  text  Ch.s 

SUBROUTINE  GAUSS (N,ZI 
INTEGER  PV 

REA’..  2  (',  ;’C2,  1  :2031  ,  E 
C«1  .0 

ID  IF  d  .O'E.G'-.  1  .01  THEN 
E--E/2.0 
GOTO  10 
END  1  F 
E*E'2 
EPS2»2*r 

FPl'.'T  ','  MACHINE  EPSILON^' ,  F, 

1005  DET'  1 

DO  lOlC  1  -1 ,N-1 
PV'- 1 

DO  102C  J'I*1,N 

IF  (ARS(Z  (PV,  II  1  .LT.  ASS(2(J,  nil  PV'U 
1020  CONTINUE 

IF  (PV.EO.II  GOTO  1050 

DO  1040  JC*l,N'l 
TM'Zd,  JCI 
Zd,  JCl  'ZIPV,  JCI 
ZIPV, JC) "TM 
1040  CONTINUE 

1045  OET'-l'DET 


3050  IF  (Zd,  II  .EO.Ol  THEN 
GOTO  1200 
FND  IF 
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Output/input  file:  bI2d.dat 


540000 . 

100  4  9 

0.99950 
0. 99150 
0.99555 
0.98165 
0.91985 
0.91015 
0.95860 
0.94525 
0.93010 
0.91330 
0.89485 
0.87480 
0.85330 
0.83045 
0.80625 
0.18085 
0,75435 
0.12680 
0.69840 
0.66925 
0.63940 
0.60895 
0.51810 
0,54695 
C. 51565 
0.48430 
0.45295 
0.42180  - 

0.39095  - 

0. 36050  - 

0.32065  ■ 

0.30150  - 

Q.2T310  • 

0.24555 
0.21905 
0.19365 
0.16945 
0.14660 
0.12510 
-1. 10505 
0.08660 
0.06980 
0.05465 
0.04130 
0 .0r9’5 
02005 
0 0 1 2  2  i 
0.00635 
0 .00240 
0.00045 

f, 

c!0O24O 

0.00635 

0.01225 

0.02005 

0.02915 

0.04130 

0.05465 

0.06980 

0.08660 

0.10505 

0.12510 

0.14660 

0.16945 

0.19365 

0.21905 

0.24555 

0.21310 

0.10150 


0.  380 


0.162 


0.00005 
0.00030 
0.00085 
•0.00110 
.0.00280 
.0.00415 
-0.00515 
-0.00155 
-0.00955 
-0.01110 
-C. 01405 
-0.01655 
-0.01915 
-0.02185 
-0.02465 
-0.02750 
-0.03035 
-0.03325 
-0.03610 
-0.03890 
-0,04165 
-0.04435 
-0.04690 
-C. 04925 
-0.05150 
-0.05355 
-0.05530 
-0.05680 
-0.05805 
-0.05900 
-0.05960 
-0.05980 
-0.05965 
-0.05910 
-0.05815 
-0.05685 
-0.05510 
-0.05295 
-0.05045 
-0.04155 
-0.04430 
-0.04070 
-0.03680 
-0,03260 
-0.02810 
■C. 02340 
•0.01850 
■0.01340 
-0.00815 
-0.00215 
0.00215 
0.00815 
0.01340 
0.01850 
0.02340 
0,02810 
C. 03260 
0.03680 
C.04C10 
0.04430 
0.04155 
0.05045 
0.05295 
0.05510 
0.05685 
0.05815 
0.05910 
0.05965 
0.05980 


0.15382 
C. 80956 
0.84584 
0.81525 
0.89440 
0.91199 
0.92801 
0.94058 
C. 95235 
0.96041 
0.91159 
0.98018 
0.98850 
0.99553 
I .00410 
1.01099 
1.01100 
1 .02456 
1 .03023 
1 .03612 
1,04113 
1 .04948 
1.05502 
1 .05902 
1.06629 
1 .01290 
1  .  07664 
1.08126 
1 .08640 
1 .09161 
1 .09584 
I .39’56 
1 .10065 
1.10082 
1.10060 
1. 10136 
1.09665 
1.09134 
1  .08489 
1  .01358 
1  .0595' 
1.03695 
1.012’! 

0. 969’0 
0.913!'. 
0.82914 
C. 10489 
0.49993 
0 .19469 
0.24808 
0.61814 
1 .03123 
1.20922 
1 .29430 
1.32604 
1.33560 
1.31815 
1.33195 
1.32113 
1.31302 
1.30239 
1.29118 
1.21918 
1.21004 
1.25919 
1.24515 
1.23599 
1  .22431 
1.21308 


0.33065 

0.36050 

0.39095 

o!4219C 

0.45295 

0.48430 

0.51565 

0.54695 

0.51810 

0.60895 

C. 63940 

0.66925 

0.69840 

0.12680 

0.15435 

0.19085 

0,90625 

0,93045 

0.85330 

0.81480 

0.89485 

0.91330 

0,93010 

0.94525 

0.95860 

0.91015 

0.91985 

0.98165 

0.99355 

0. 9»15C 

0.999SC 


0.05160 

!  .2r2'e 

0.0590. 

; .190:1 

C . 05  80  5 

i  .  1 ' 9  4  ■ 

o!c569: 

1  .  1  6662 

C .05530 

1.11685 

0. 05355 

1 . 1 4  - :  3 

0.05150 

1.13455 

0.04925 

1 .12330 

0.04690 

1.1 1552 

0.044)5 

1 .10402 

0.04165 

1. 09349 

0.03890 

1.08354 

C. 03610 

1.01480 

0.03325 

1.06496 

0.03035 

1.05419 

0.02150 

1.04603 

0.02465 

1 .03433 

0.02185 

1.02466 

0.01915 

1.01330 

0.01655 

1 .00401 

0.01405 

0.99028 

O.Oll’O 

0.9’9I» 

0.00955 

0.96712 

O.OC155 

0.95331 

0.00515 

0.93161 

0.00415 

0.91836 

0.00280 

0.89968 

0.00110 

0.81407 

0.00085 

0.83918 

0.00030 

0.79199 

0.C0005 

0.1S392 
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Source  code:  bl2d.r 


KodUlr«tlon*i  L.  M.  NOHAK 
V€r .  ? 

It  July  1»9J: 

Rdded  wtll#  itattititntt (30*40l  In  th*  “00  Pi’  loop  to 
output  boundary  layar  proftlai  for  plotting 
<100  panal  airfoil  cnlyl 

20  July  1992:  Addad  wrlta  atatamanta420-21)  to  output 
Cf'lskln  friction!  and  OLSIdelta  star!  for  plotting 
ver.  3 

SI  Aug  1992;  Added  calculation  for  oniet  of  tranaltlon 
bnied  on  Michel* a  criterion,  added  Input  ITAANS 
ver .  4 

3  Sept  1992:  Modified  boundary  layer  profile  output 
to  be  compatible  with  airfoil  of  any  number  polnta 
ver,  i 

3  Sept  1992:  Aedlmeni ioned  all  arrays  to  accept 
airfoil  up  to  200  panels 


SUBROUT  INC  BL 

COMMON  /BLC2/  NX, NXT, NP, NPT,  NTB,  IT,  ISC 

COMMON  /BLC3/  X  (2001  , UC (2001  , Pi <2001 , P2 12001 , CMTR (2001 

COMMON  /BLC7/  ETA (2C1 1 , DETA (201 1 , A (2011 

COMMON  /Bice/  r(2Cl,2),U(201,2l,V(201,2l,B(201,21 

COMMON  /BLC6/  DELF ( 20 1 1 , OELU (201 1  .  DELV (201 1 


NX  -0 

ITMAX  -  10 
ICROHT  •  2 
EPSL  =  0.0001 
EPST  »  0.01 
NPT  •  101 

ETA-CBIO 
ETAE  ■  8.0 
VCP  ■■  ”1.10 

DETA(l)  ”0.01 

NF  •  AL0G(  (ETAE/DETAdl  I  •  (VGF-1 .01  *1 .01 /ALOGIVGPI  ♦!  ,00! 
ETA (11  «  0.0 
DO  10  J-2,NPT 

ETA(J1  •  ETA(J-11  *  CrTA(J-ll 
r!ETA(Jl  »  VGP*DETA(J-3 1 
A(JI  •  O.S*OCTA(J-11 
CONTINUE 


INITIAI  LAMINAR  VELOCITY  PROFILE 
DO  20  J’l,Nr 

ETAB  •  ETA(JI /ETA(NPI 
ETAB2  •  ETAB* *2 


F(J,2I 
U(J,2) 
V(J,21 
8(J,  21 
CONTINUE 


0.25*ETA(NPI 'ETAaF* (3.0  - 
O.i'ETAB* (3,0  -  ETAB2I 
l.i'd.O  -  ETAB21 /ETA(NPI 
1  .0 


C.i'ETAB?! 


IT  •  :t*i 

IF  (IT  .GT.  ITMAXI  CO  TO  101 
IF  (NX  .CE,  NTRI  CALL  EDDY 
CALL  COEF 
CALL  SOLV3 


CHECK  FOB  CONVERGENCE 
IF  (NX  ,LT,  NTRI  THEM 
IF(ABS(DELV(1II  ,GT. 
ELSF. 


EPSL I  CO  TO  5 


•»2 

:f<absioelv(1i /v(:,2>  >  .gt. 

f3 

ENCIF 

T4 

c 

75 

n 

FRCFI’.  ES  TCP  GROWTH 

76 

DC  30  J*NP*1,NPT 

11 

r(J.2)  *  F(J-i.2)  •  DETA(J- 

79 

U(J,2}  =  UfJ-1.2) 

79 

VU\23  =  0.0 

9C 

B<J,2)  =  BlJ~l,2) 

ei 

30 

CONTINUE 

82 

c 

83 

0 

CHECK  FOR  CRCWTH 

84 

IF  (AP<(V(NF,2) )  .GT.  C.CC05  . 

85 

•  .GT.  C.005)  THEN 

86 

N?  J  NP*? 

87 

IGPCW  =  IGRCW*! 

88 

IF  (NF  .LE.  NPT  .AND.  IGPOW 

89 

IT  =  0 

90 

GO  TO  5 

91 

END  IF 

92 

FNDIF 

?3 

c 

Q4 

'' ' 

ca:l  c-’?ut 

!f  (NX  .IT.  NXT)  GO  TO  1 

96 

r 

9 

PET'RN 

99 

END 

99 

SVPRCC'ISE  CCEF 

COHIMON  aici-  NX,  SXT,  S?.  S7  T,  ST 

;  *: ; 

CCr-WOS  '9ic3/  X  (2003  ,vt  (200)  ,  r 

>  < 

rOMMCS  .BIC-'  ETA(?0:),-O-TA(;0 

-  3 

r:«“CN  .  eoc9  r»2c:.2i /.co:,? 

:  :4 

co'^cN  ,'9!,c9-  sn2r: I , s;  co:  • . 

;  T  i 

.  S*  (201 5  ,  G°  , 

I  C  f 
1C7 

- 

;  09 

F>  *  0.5  •  P:  (NX) 

I  0  9 

:*•  INX  .fC-  ‘3  TMES 

I)  -  c.c 
’  1  -  C 
1 }  s  c !  0 


':cs’:nj£ 

r:sr 

cr:  ^  C.5  •  <x(sx) ‘XiNx-:? '  jx (nx- -x»sx-: * j 
C':‘'=  c.^  •  cii 

r-  :::  ?.np 

C'Scfvt  STATICS 

F3  »  o.5MF(:.2j  *  n:-:,:}'. 

i'?  =  -  V(j-:,:n 

FVB  *  0.5*  (F (j. 2)  *r <j-i, ‘v*  .:m 

VB  *  C.5MV(J,  2)  •*  V(J-:,2M 

•jsg  *  C.s*  -’.2;  *T  . 

rrpgv  'v<c,?)  -  p<j-i,23 


PBEVICVS  STATION 

era  =  c  .iMF  ( j, !)  *  F  ::  i 

CUB  *  c.5M'j<-MJ  ♦ 

CV9  =  c.5*(v(j,i)  ♦ 

C'v’SB  =  C.5MU(J,n**2  •* 

ervB  ^  c.  5*  <r  ^ j, ! Mv^j,  n  *F < J-i .  I ;  •vfj-i . :  n 

CCER9V  =  (B ( J, 1) 'VJ J, ! »  -  a ( J-1 , i } 'Vf J-I , I J )  TETA f T-1 ) 


S-  CCEFFICIFNTS 

51  IT)  ^  CEIH*  if  '  -  CfSl 

52  (/:  ^  ce:.“* 


P:hT(J,2)  *  a(J,?>  /TFTAJ 
r  :ht  »  j>i  ,  ?)  -b ( ~r: 
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:r(M»5(orLvm/v(i,?)i  .gt.  cpst)  ca  to  s 
tNDIf 


fRCi:;  t5  TOP  CBOWTP 
nc  30  ,?*NPM,NPT 
F(J,3I  ■  r(j-l,2l 
U(J,2)  .  U(J-1,?) 
v(j,j)  •  0.0 
B(J,21  * 

30  eONTINUr 


nfTA(J-ll *0(3-3, ?l 


CHTCK  rcR  cp^wth 

ir  (ABS(V(NP,  J))  ,CT.  O.OOOi  .OB.  ABS  (1 .0-0  (fJP-J ,  7 1  (UP,  ? !  I 
•  .GT.  O.OOA)  THTN 

NP  MP.? 

joBow  ■  !cncw*i 

!r  (KP  .’..F.  F'PT  .AMD.  JCPOW  .OF.  IGHOHT)  TPFM 
!T  •  0 
GO  TO  5 
EKDTF 
FNDIF 


CA!,0_  0I‘  rOL'T 

c.'x  .OT.  go  to 

BETl’RN 

END 

SOBRCUTINE  COEF 

COMMON  /BLC7'  NX,  NXT, NP, NPT, NTP, IT, ISF 

COMMON  /B0C3/  X(?00|,'JE(200),P1  (2001,P7  (200),SMTB(?0C» 

COMMON  /BLC7/  ETA (203 ) , DETA (201 » , A (201 ) 

COMMON  7BI.C8/  E  (201 , 2 1 , 'J  (201 , 2 1 ,  V  (201 , 2) ,  B  (201 , 2) 

COMMON  /BLC9/  Si (201 1 , S2 (201 1 , S3 (701 ), S4 (201 », S5  (201 1 , S« (20! ) , 
•  ST (2011,5?  (70!), PI  (201),R2(70!I,R3(701).R4(;C1I 


PIH  «  O.S  •  PI (NX) 

IF  (NX  .EO.  1)  THEN 
CEL  »  0.0 
CELH*  0.0 
DO  5  J*1,NP 
r(J,ii  •  0.0 
U(J,1I  •  0.0 
V(J,1I  •  0.0 
B(-M)  •  C.O 
CONTINUE 
ELSE 

CEL  •  0.5  *  (X(NX)-X(NX-1I1/(X(NXI--X(NX-1II 
CELH»  0.5  *  CEL 
END  IF 

DO  IOC  J*  2,NP 
CURRENT  STATION 

FB  -  0.5*(F(J,2)  ♦  F(J-l,2n 

UB  ■  0.5*(U(J,2I  •  U(J-1,2)I 

FVB  •  0.5*(F(J,2l*V(J,2(‘r(J-l,2)*V(J-l,2)) 

VB  •  0.5* (V(J,2)  *  V(J-1,2) I 

USB  •  0.5* 'U(J,2) •*?  •  U(J-1,2)**2I 

DEBSV  •(B(J,7) *7(3,7)  -  B (3-1 , 7 ) *7 (3-1 , 7) ) /DETA (3-1 ) 

PRE7I0U5  STATION 

CFB  •  0.5*(F(3,1)  «  F(3-l,l)) 

CUB  •  C.5*(U(3,1)  ♦  U(3-l,l)l 

C7B  ■  0.5*(7(3,1)  ♦  7(3-1,11) 

CUSB  •  0.5* (U(3, 1) **2  «  0(3-1, 1)**2) 

Cr7B  »  0.5* )F(3,  1) *7(3, 1) *r(3-l, 1) *7(3-1, 11 ) 

CDERB7  .  (B(3,l) *7(3,11  -  B (3-1 , 1 1 *7(3- 1 , 1 1 1 /DETA (3-1 1 


S-  COEFFICIENTS 

Sl(3)  •  CELH*(F(3,2)  -  CFB)  ♦  P1H*F(J,2)  •  B (3, 7 ) /DETA (3-1 ) 

S7(31  •  CELH* (F'3-l,2)-CFe)  ♦  PIH *F (3-1 , 7) -B (3-1 , 7 ) -OETA 13- 1 ) 
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M3 

51('!  •  CELM*  (V|..',3)  •  C'lBI  •  MM*'.'(J,2t 

M4 

MU)  .  CEl.K' (VU-1,  2)  •  rVBl  .  P1H*VU-1,21 

MS 

E5(.!)  .  .(cr.i  *92  (NXl  1  'UU,  2) 

14( 

56U)  .  -(CF*  -r,’ l!)X)  )  VU-l ,  2) 

Ml 

c 

MS 

r 

p-  ccEmcirNTs 

MS 

IP  (NX  .FO.  11  THEN 

!S3 

CRB  •  *P2fNX) 

15! 

92  (J)  •  CRP  -  (OERBV  .  PI  (NXl  *rVB  -  P2(NXML'JRt 

1  52 

K1.5F, 

153 

COB  •  CPERR'2  .  91  (NX-ll ‘crVB  -  92  (NX-1 ) -CUEB  •  P.’CIM-ll 

!  54 

-CRB  •  -COB  -  CEOT’JSB  -  92  (NX) 

155 

P2U)  >  CRB  -  (DEPBV  •  91  (NX)  •  F'.'B-  (CEI  »P2  (NX  1  1  "J'.B  •  CEO 

156 

•  (FVB  •  CVBTB  -  VB*CFB  -  CFVB)  1 

157 

ENOIF 

158 

B1  (J1  •  F(J-1,2)  -  F(J,2|  •  OETAU-lMOB 

159 

R3(J-1M  U(J-1,2)  -  UU,2I  ♦  DETA(J-'.)  n’B 

160 

100 

CONTINUE 

Ml 

C 

1  62 

r* 

BOUNPABY  CCHniTlnNS 

•6  3 

f 

164 

91(11  '  0.3 

165 

R2(:)  c.c 

1  66 

R3(NP|  •  C.r 

HI 

c 

'.68 

RETURN 

169 

ENP 

1  !C 

COMMON  /BOOC.  RO.  NBO  ( 2 ! ,  XTT9I  (2) ,  f  lag,  •- raBsi*- (? ) .  ’  ' 

COMMON  /BOCi/  !rR,XCT9,XC(2C01,YC(2C0l 

1*12 

COMMON  /BIC2/  NX, NXT, NF , NPT, NTP , IT, ISF 

1  ?3 

COMMON  '91.C3,  X  ( 2 30) ,  UE  (2001 , 9 1  (2001 , 92  (200) ,  SMTP  (.733! 

1  14 

COMMON  /0I.C4'  d:,S(2301  ,VW(200),Cr(230),THT(2CC) 

1  15 

DIMENSION  Nxrsf  (2l , Xl (200) , YI (200) ,VEI (200) 

116 

c 

1  n 

c 

??? 

1 18 

CFF.N  (UNIT  -  9,  Fi:  f  'b:2d.dat’ ,  STATUS-' 'UNKNOWN' ) 

1  J9 

CFF.N  (UNITES,  FI'.F  .  'tl 2d. out ' ,  STAT'JS*' UNKNOWN' ) 

1  BC 

OPEN  ( UN !  T •  2  0 ,  r  1 E  ’ '  C  f  . da t '  ,  STATUS » '  'JHKNOWN' ) 

181 

OPEN  (UNIT";  1,  FILE-  'dla.dat ' ,  ST.XTUS*' UNKNOWN' ) 

162 

1P3 

WRITF.(6,  •)  'REA0IN3  THE  DATA,..' 

184 

READ  (  9,15  )  RU.XCTRI (11  ,XCTRI (2) 

1  85 

READ  (  9,10)  Nl, IS, ITRANS 

:  96 

READ  (  9,H  1  (XHIl.YHII.VEKD.Ml.NII 

191 

WRITE(6,'I  'INPUT  OF  DATA  COMPLETE.' 

1  98 

c 

1  9  9 

r 

j «: 

WRITE  (8,  90  PL,  XOTR!  l! )  ,  XCTPl  (21 

191 

NXTSFdM  NI  -  IS  • 

i  92 

NXTSr(2)«  IS 

1  93 

z 

194 

c 

BATA  FOR  each  S-URfacE 

1  95 

DO  20C  ISF  •  1,2 

1  96 

ntflad*0 

191 

NXT  •  NXTSFdSFI 

I  98 

CO  TO  (201,2021  ,  !S>- 

1  99 

c 

2C3 

r 

UPPER  SURFACE 

2C1 

23! 

II  •  IS-1 

2C2 

DO  211  Ml,  NXT 

203 

n  «  iMi 

204 

XC(I)  •  XI  (II.) 

205 

YC(I)  '  YI  (II) 

206 

UEC)  «  VEI  (11) 

207 

211 

CONTINUE 

208 

GO  TO  300 

209 

f 

210 

f 

LOWER  SURFACE 

211 

202 

11  •  1S*1 

212 

DO  212  Ml, NXT 

213 

11  •  II-l 

112 


;i4 

xt'i;!  xn!!! 

?1S 

YCi  1 1  T I  u :  1 

:u 

ufii)  -  vridii 

717 

212 

roNTiNur 

71B 

c 

7:<t 

300 

XU1  »  0.0 

720 

00  30]  i^2,Nxr 

721 

101 

X(I1  •  X(!-ll  •SORT!  (XCItl-XCd-ll  1  ••2*(Y'--(n -Y' <1-11  )  ••21 

727 

c 

773 

c 

77< 

c 

TRANSITION  l.CCRTION 

77i 

DO  320  I-.1,NXT 

776 

CMTRin  •  0.0 

777 

ir  ixc(i)  .at.  xCTRidsm  co  to  121 

276 

120 

CONTlNUf 

779 

321 

NTR  •  I 

730 

PCAMTP  ’  1200. 

731 

RXNTR  '  X(NTP-11*  UE(NTR-1I  •  RL 

337 

CCrT  •  RL*  *2/RXNTR“  1 .  34*UK(NTR-!  )  **3 

733 

UEINTG  •  0.0 

734 

III  ^  0.5/UE(HTR-l)  /  fGAMTR 

23i 

no  322  T  •  NTP.NXT 

236 

U2  .  n.5/L'E(I) /FCAMTP 

231 

lIF.INTC  ’  L’EINTO*  CJl -Ol)  *  <X<  ! ! -XI '-!  )  ) 

236 

L’l  .  U2 

234 

CC  "  CGET-OEIKTC* (X 1 :i -X(NTR-ll  1 

24C 

iriGG  .GT.  IC.O)  oo  TO  323 

741 

GMTRd)-"  1. 0-EXP  l-GG( 

747 

322 

CONTINOE 

743 

323 

DO  324  1I«I,NXT 

344 

324 

nMTPdll  •  l.C 

34i 

c 

246 

c 

PRESSURE  GRADIENT  PAPAM'-'FPG 

241 

DX  «  X(2)-Xf:i 

246 

DUE  .  UE(2)-UEC.t 

249 

ANG2  •  ATAN2(DUE,rX) 

250 

DL2  •  OX 

25! 

no  331  I  •  2,NXT-; 

252 

AKGl  •  ANG2 

253 

on  ■  DL2 

254 

OX  »  Xddl-.Xdi 

255 

DUE  •  UE  d  ♦ ; )  -'jf  ( : ! 

256 

ANG2  •  ATAN2 (DUE.DXI 

257 

OL2  »  CX 

259 

ANG  ’  (CL2*A!.Gt -D:  1 'ANGli  1011-012: 

259 

P2(II  «  TAN' (ANG) 

260 

331 

CONTI NUF 

26! 

P2(NXT)  ’  2.*D'JE/DI2  -  P2(NXT-n 

262 

DC  3  30  I  ■<  2,NXT 

263 

P2dl  •  Xdl  •  PPdl  /L'Ed) 

264 

Pldl  *0.5*d.0-P2dH 

265 

330 

CONTINUE 

266 

P2d)  >  l.O 

26-’ 

Pi  d!  «  0.5  •  (1 .0  ♦  P2  (1 1  1 

268 

c 

269 

c 

210 

c 

BOUNDARY  LAYER  CALCULATION 

27! 

WRITE (6, *1  'BOUNDARY  LAYER  COMPUTATIONS  IN  PROGRESS..’ 

2  72 

CALL  BL 

273 

WRITE  (8,  910)  ISF,  d  ,  XO  (1 1 ,  X  d ) ,  VW  (I ) ,  CP  (I) ,  DLS  d ) ,  THT  (1 1 ,  Id.NXT'. 

2'4 

IfdSr.EO.l)  then 

215 

wrlte(2C,  9C5)  (XCd)  ,CPd) ,  :  =  2,NXT) 

216 

write (21, 905)  (XCd)  ,DLSdi .  !»P,NXT) 

277 

end  1 1 

278 

905 

rCRMAT(PB.4,4X,EI1.4) 

279 

200 

CIJNTINUE 

280 

c 

281 

• 

•  ir 

ADA  Is  0  deg.,  iraXe  trans.  Iocs,  eoual: 

282 

lf(vel(2).eq.v«l(nl-l))  trsnsnew(l)«transnew(2) 

283 

294 

1  f  dTPANS.eq.O)  then 
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prl"!  *,  *  Cfl  tirat*  for  urp«f  t  rs»!!  15  st-)  ’ ,  t  ( i ) 

rrlnt  for  lowar  t  rantU  icr  ,  t  rar.fippw  I2 1 

791 

t 

799 

n0Sf,<UHfT<81 

J90 

CL0Sr.aiNlT«9) 

J91 

STOP 

797 

r 

791 

10 

rORMAT(3l5) 

799 

15 

rORMAT(3«'10.0> 

;95 

9C 

rCRMAT(//5X,'RL<-*  ,  Sll 2 . 5,  5X,  '  XCTRI  (11  ,  7  8. 3, 5x, '  XCT8  <2) 

:9* 

910 

rORMAT(///2X,  ‘  SUHMABY  OP  BOUHDARY  LAYER  SOL'JTICNS  OP 

791 

♦//2X,  'HX'  ,  4X,  •  XO*  ,  9X,  '  S'  ,6X, '  VH*  ,  8X,  'CP'  ,  BX,  'OLS'  ,  SX, '  T« 

2  98 

♦  /(:5,2pR.4,4FS1  ,411 

299 

CNU 

300 

SUBROUTIKE  CCOY 

301 

COMMON  /BLCC  •  RL.  NBL  (21 . XCTRl  (21 ,  nt  f  la<j,  t  ransnrw  ( 21 ,  N t 

302 

COMMON  /BLC2'  NX, MXT, NP, NPT, NTR,  IT,  ISP 

3CJ 

COMMON  .'BLCS'  X  ( 2001  ,  L'C  ( 2001  ,  PI  (700) ,  P2  (2001 ,  GMTR  (2001 

304 

COMMON  /BLC7'  ETA (201 ) , OETA (23! 1 , A (201 1 

305 

rCKMON  .'B1.CS/  '■(201,21  ,0(201, 7). V(701, 2). "1201, 21 

3C« 

O'.MPNSICN  F.DyH20n 

307 

c 

3C8 

r 

309 

RL2  •  SQP:(Bl"L'F.(NXI ‘XlNXll 

310 

B14  '  S0BT(R!.21 

311 

RLPlf  ’  0.16  •  RL2 

312 

c 

313 

ALFA  »  0.0168 

314 

EO'/O  '  AIP4*Rtr*0MT»(NX1  •(’v’(NP,2)*ETA(NP1-PtNP,2l1 

ilS 

EOVKll*  0.0 

318 

YBAJ  •  BL4*S0®T(ABS(V(l,2n  1 /76.0 

317 

00  7C  j'P.NP 

318 

Jw  *2 

31? 

YBA  •  YBA."ETA(J) 

320 

EL  •  l.C 

321 

IP(YBA  .IT.  10.01  EL  •  1.0  -  EXP«-YBA) 

322 

EOVKJl  '  Rt,216*OMTR(NX)  •  (Pt.'ETA(01  1  '  APS(V(;,2ll 

323 

IF(E0V1(J1  .OT.  r.DVO)  00  TO  90 

324 

IP  (EDVK.'l  .LE.  E3'/:(J-1)1  EOVI  (J)  •  POVKJ-ll 

325 

8(J,2)  •  1.0  •  EOVKJl 

32« 

70 

CCNTIN'JE 

327 

90 

CO  ICO  .7J"J,HFT 

328 

100 

2(J.’,21  '  l.C  •  PDVO 

329 

B(l,7t  '  1.0 

330 

c 

331 

RET'JRN 

332 

ENO 

333 

SUBPOUTlf.T  CVTPUT 

334 

COPWCN  /»LCC/  RL,HBL(2I  ,XCTRK2)  ,nt  f  laq,' ra-sr-o. 

335 

COMMON  /BLC2/  NX, NXT, NP, NPT, NTR, IT,  3 SP 

336 

COMMON  /D1.C3/  X  (200) ,  UE  (2001 ,  PK200I , F2  (2"0I , GV"  1 2 '  " 

337 

COlWON  /BLC7/  ETA(201),DETA(201),A(2e'.) 

338 

COMMON  /BLC8/  P  (201 , 2)  ,  IJ  ( 201 , 2 ) ,  V(2C1 , 2 1 ,  B  (231 . 2 1 

339 

COMMON  /P'.CS/  'JL5(2001  ,VH(200)  ,CF1200),THT(20CI 

340 

dlnenslon  rdl  I  K2C1) ,  rdlow  (201 ) 

341 

c 

312 

c 

343 

IP(NX.E0.1  )  THEN 

344 

DLSINXI*  0.0 

315 

THTINXI-  0,0 

346 

CF(NX)  ■  0,0 

347 

VW(NX)  •  V(1,2I 

348 

tdlKlowlOOO 

349 

nstop’O 

350 

ELSE 

351 

c 

352 

SOBX  •  SOBTIL'E(NX) ‘XINXI 'RD 

353 

CP1NX)  •  2.0  •  VC., 7)  •  B(1,2)  /SORX 

354 

VW(NX)  '  V(l,7) 

355 

CLSINXI*  X(NXI/SCRX  •  (ETA(NP)-r(NP,21) 

114 


Pages 


ui  '  u(i,?i  •  (1.0  -un, 211 

351 

SOM  •  0.0 

35t 

DO  20  .7'?,NP 

359 

U2  •  U(J,2l  •  11.0  -U(J.2)  1 

3iO 

SUM  .  SUM  *8(01  •  R'l  *  02) 

3il 

01  >  02 

3(2 

20 

CONTINUE 

3(3 

TMT(NX)*  XtNXl/SORX  •  SOM 

3(4 

r*li>UC(NX)  *X(NX)  *81 

365 

rth»t«*Uf (NXl *THT(NX) ‘BL 

3(( 

tt r«ni"l . 174* (1 .0<2?400.0/reK) ‘to** •0.46 

361 

rdl!f(nx)>«bs (rtheta-rtrani) 

3(8 

If  nNX.qt..2)  .and.  (rdiow(nx-l)  .eq.rdlowtnx-?) )  I  than 

J69 

If  ( rdlow  ( nn-?  >  .  cq.  rrtlow  (rK-3) )  nstepi'l 

endl  f 

371 

If  (3Sr.aq.2)  than 

1'? 

If  ( (nt f lag .eq. 1)  .and.  (natop.eq.Ol )  than 

373 

ransncwdSFI  •reK/(Rb*UE1NXl) 

374 

nt ;  laq*0 

37S 

endl  f 

37« 

end!  f 

377 

1 f ( (rdlf f (n*l  .LT.rd! fflow)  .and.  (natop.aq.O) I  than 

378 

transnewlISri  Tex/  (RfJElNX) ) 

379 

rdlf  f  low-'rdl  (1  (rxl 

38r 

nt  flag*! 

381 

cndif 

382 

rdlow (nxt "rd! f  f lew 

383 

END  IF 

384 

c 

print  • , natop, nx, rdl f f (  nx-1 ) , rdl £f (nx-?) 

385 

c 

print  *,  1 sf, rex, c ransnawliaf ) 

386 

C 

print  *, rt beta, rtrans,rdiff(nx),rdlfflcw 

387 

c 

388 

c 

389 

c 

SHIFT  PROFtLES  FOR  THE  NEXT  STATION 

390 

c 

391 

yiMrXx.OOOS 

392 

DO  175  J-l.NPT 

393 

IflISF.CO.n  then 

394 

lf(U(J,n.LT. '0.995)  1  then 

395 

laaty‘1 

396 

yp!ot*ETA(J) •SORT»X(NX)/(RL*UE(NX)) ) 

397 

do  nxloop»5, Nl/2*!, 5 

398 

1 f (NX.CO.nxlocp)  then 

399 

91 

warXxwSX/S 

40C 

numw*marxx*  30 

421 

write  (nu">w,*)  U(J,  1)  •marXx, yplot 

402 

write  (60,*)  U(J,  1 )  •r.arXx,  yplet 

403 

I f fyplot .qt .ypark)  then 

104 

write{55,'l  maikx.ymarlc 

405 

ydi f  f»yplot-yplotoid 

406 

udiff»U(J,l)-U<J-l,l) 

407 

xva j  ue’U ( J-1, 1) ‘udlf f • (ynark-yplotoldl /yd; f  f 

408 

write(55,*)  Kvalue*narkx, ymark 

409 

wrlte(55,921 

410 

9? 

format  (/) 

111 

ymarkwymark* .0005 

412 

if  (yplot.CT. ymark)  goto  91 

413 

endl  ( 

414 

endl  f 

415 

end  do 

416 

•  Isc 

417 

if  (laity. CO. 1)  then 

418 

lastywO 

419 

do  m'1,3 

420 

do  nxloop*5,HI/2-1, 5 

421 

1 f (NX.CO.nxloop)  then 

42? 

markx'NX/5 

423 

nunw>markx*30 

424 

write  (numw, •)  markx,yplet 

425 

write  ((0,*)  marks, yplol 

42( 

end!  f 
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427 

end  do 

428 

yplot -0.0 

429 

end  do 

430 

end!  f 

431 

endl  i 

432 

endl  f 

433 

yDlotold=yplot 

434 

F(J, 1)  =  F(J,2)  j 

435 

U(J, 1)  *  U(J,2)  1 

436 

VU.l?  V<J,2) 

43? 

B(J, 1)  »  B(J,2) 

438 

175 

CONTINUE 

439 

C 

440 

RETURN 

441 

END 

442 

SUBROUTINE  S0LV3 

443 

COMMON  /BLC2/  NX, NXT, NP, NPT, NTR, IT.  I SF 

444 

COMMON  /  BLC?/  ETA <201 7 , DETA 1203  ) ,  A (203  ) 

445 

COMMON  /BLC8/  F ( 201 , 2 1 , U  (201 , 2) , V (201 , 2 » ,  fi 1 201 , 2) 

446 

COMMON  /BLC9/  Si (20 1 ) , S2 (201) , S3 (2011 . S4 (201 1 .  S5 ( 20 1 )  ,  S6 (201 )  , 

447 

♦  S7 (201) ,sa (201) , R1 (201) , R2 (201) . R3 (201 J , R4  (201) 

448 

COM.MON  /BLC6/  DELF  ( 201 )  ,  DELC  (201 ) .  OELV  (2C1 ) 

449 

DIMENSION  All  (2C1) , A12(2C1) , A13(2C1)  .  A14 (20!)  , 

450 

♦  A21  (201 )  ,  A22  (201) ,  A23  (2Cn  .  A24  (2C1  > 

451 

c 

452 

c 

453 

All (1)  =  1.0 

454 

A12(ll ^  C.C 

455 

A13(l) *  C.O 

456 

A21 (1) *  0.0 

457 

A22 (1) =  1  .0 

4  58 

A2  3 ( 1 ) »  C.C 

459 

Gil  ^-1.0 

460 

g:2  =-a(2, 

461 

G13  *  0.0 

462 

G21  -  S4(2) 

463 

G23  *-s:(2)/A(2) 

464 

G22  *  G23-S6(2) 

465 

All (2) «  1  .C 

466 

A12  2)  ’-Ai:? -013 

467 

A13(2)  '  A(2;  'cn 

468 

A21  (2)  --  S3  (2J 

469 

A22(r» -  S5(2)-G23 

470 

A23  (2j  s  SI  (2)  •A(2) *023 

471 

91(2)  »  r: (2) - (Gl 1 'RI (I ) •G12'R2 (1) *G: 3*R3 » f 

4  72 

92(2)  •-  P/ (2)  -  (G21  *R1  (1 )  •G2?*R2  (1)  *G2.*  *9^  C  5  i 

473 

c 

4  74 

c 

FCRWAPP  SWEEP 

<73 

- 

476 

DC  50C  2-2, NF 

f.  77 

DEN  •  (An ( j-n  •a2i  (2-i)-A2? (j-i) ‘a: :  <j-i : -AC)  • 

478 

*  (A:2  (J-l )  •A21  (2-1)  -A22  ( J-15  ‘A*,  n/- :  )  J  J 

479 

DEN!  *  A22 (J-I) *A(2) -A23(2-l) 

480 

GU  '  (A2  3  (2-1)  ‘AfO)  •  (A{J)  'A?!  (2-1 5  -A22  (  :-V.  )  )  /CEN 

491 

G12  ^  -  (Al  J1  ‘AiJ)  'Gin  (Al2  ( J-l)  •A(2)  -Al  3(2-1  ’■  i  )  .'DENl 

492 

G13  *  (G:nA13(J-l)*G12*A23(2-in /A(2) 

493 

G21  =  (S? (2)  •A21  (J-l) -S4 (J) •A23(2-:) •A(2)  • (S4  <21  • 

484 

*  A22  (J-U  -S6(J)  -A?!  (J-l)  )  }  /DEN 

495 

G22  -  (-S2 (2) ♦S6(2) 'A (J)-G21* (A (2) 'AI2 (2-l)-Al 3(2-111)  'CENl 

496 

G23  ^  G21*Al2(J-mG22*A22(2-ll-S6(Jl 

487 

A]1(J)=  1,0 

488 

a:2 ( J) =-A(J} -Cl  3 

489 

A13(J) =  A(J) *013 

490 

A21(J)=  S3(J) 

491 

A22(J)s  SS(J)-G23 

4  92 

A23 ( J) =  S3  (J) ♦A ( J) *C23 

493 

PI (J)  »  R1 ( J) - (Gll 'Rl (J-l) *312 •R2 (2-1 1  -G!  3*R3 ( J-!) ) 

494 

R2 ( J)  =  R2 (J) - (G21 'Rl (J-l) *022 ‘R? (J-l) *02 3*93(2-1) 1 

495 

500 

CONTINUE 

496 

C 

497 

C 

BACKWARD  SWEEP 
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Output  file'  bi2d  out 


• 

1- 

f  ‘  e.)«: 

j  L  ,  ’  • : 

... 

Of  •CVlit 

ACT  lAlf  »?tVT!r%3  Zt 

!Sf  -  1 

Ml 

1C 

s 

Vli 

Cf 

?LS 

7»«T 

» 

0  CC24 

C.12)2t*51 

o.oocct*;o 

C'  -CC 

: . 

0.0CC5 

0 , C05  7 

0.845iC  CC 

C.4C98t-C: 

c. 

) 

0.0005 

0.0112 

C. 1242r>01 

S. J87)8-C! 

0,12158-:! 

C,54^4£*C4 

4 

:.CC24 

O.Ol '0 

s.:j8)c*C) 

0. 22288-01 

o.i2:nz} 

C  .  5i99t-C4 

3.0564 

0.02)6 

0.8CltC*CC 

0.  :292t*:; 

C  l454t'C5 

C  ,  4995C-C4 

4 

3.0123 

0.0)14 

0.652Ct«Ce 

0.88068-02 

C.3U4C03 

C.f 94fe-C4 

7 

0.0701 

0. 04C4 

0.5MU*OC 

0.6C)38-:2 

0.28286-0) 

0. :  UU-C3 

• 

0.0268 

0.051) 

0.443U*C0 

0.68211-0; 

C  )4’9€-0) 

e.:  )77t-c) 

9 

C.041) 

0.0437 

0.40231.00 

C.3’»49C-02 

c.4;:iC"C3 

:  *  u  3 :  c  - : ) 

10 

0.054  7 

0.0777 

0.3660e.00 

0.30968-02 

C.4944t-03 

c. :9c:e-c3 

11 

0.0498 

C.09)4 

0.32268.00 

C.;6988-02 

C.5920t-0) 

C.21858-C) 

I? 

C.0844 

0. 1106 

0. 3097C.0C 

C.:2I2£-02 

C.66!3£-0) 

C.2<4'E-C3 

:  3 

c.r5o 

C. 1293 

C.2804C*30 

0.  :e4CE-02 

C  .  '49CE-:3 

C.2’55£-03 

14 

0.1. j! 

0.1496 

0.2644E.CO 

0.161 '£-o; 

C.8  38  ‘£-23 

0. 3C51E-C3 

15 

:  .  1  444 

C.1712 

0.24268*00 

O.HllE-C; 

2. 9339E-D3 

C.  3  35  3E-C3 

U 

0.1495 

0.1941 

0.23938*00 

C.lBllE-C*. 

:  .I02Cr-02 

C. 3646E-03 

1  ^ 

0.10)7 

0.2104 

0.20918*00 

C.  10B5E*C: 

c.:‘7'ro7 

0 . 3954E-C3 

18 

C  .  2  !  9  ! 

0.2438 

C.!926t9-00 

0.9513E-03 

o.;24;*:-:; 

:.4270E-O3 

19 

0.2456 

C.27C4 

C.1784t*C0 

C.8401E-C3 

0.: 349E-C2 

:.4513£-C3 

2C 

C,2731 

0.2979 

0.16088*00 

C. 72698-03 

C.:469E-C2 

0.43-H-03 

21 

c.3c;5 

0.3263 

0.1398E*CC 

0.6O46E-0J 

C. ',  6028-02 

C.:<< '--03 

22 

C. 3307 

0.3555 

0.12028*00 

C.5001E-C3 

C.  : 746E-02 

c .  559:r-:3 

23 

r 

0.3853 

0.74108-01 

0.2978E-C3 

C.1965E-C2 

0.59408-03 

24 

0.39:9 

0.4158 

C. 4473e-Cl 

0.17398-03 

C.2!',E-:2 

0.63158-03 

25 

0.4218 

0.4467 

0.11778*00 

0.64398-03 

994E-:2 

:.6692E-C3 

26 

0.4530 

0.4778 

0.28238*00 

0.!C33E-C? 

C. : 746E-C2 

0.70358-03 

27 

C,4843 

0.5092 

0.43088*00 

0.1534E-C2 

C.:655E-C2 

0.75208-03 

28 

C.5156 

0.5407 

0.56888*00 

0.:976E-C2 

C. : 64rE-C2 

0.8127E-03 

29 

0.5469 

0.5720 

0.72468*00 

0.2460E-C2 

C.1656E-02 

0.87948-03 

3C 

0.5791 

0.6033 

0.85748*00 

0.28458-02 

C.l  ■’CCE-C2 

C.9503E-03 

31 

0.6090 

0.6342 

0.9414£*OC 

0.30628-02 

C . !  '92E-C2 

0. 10318-02 

32 

0,6394 

0.6648 

0.10178*01 

0.32478-02 

C.1991E-C2 

0. 11198-02 

33 

0.6693 

C.6948 

0.10778*01 

0.33798-02 

0.19878-02 

0. 12068-02 

34 

0.6994 

0.7241 

0.11198*01 

0.34548-02 

0.21008-02 

C.1296E-02 

35 

0.'»269 

C.7S26 

0.11338-01 

0.34438-02 

0.22098-02 

0.13908-02 

36 

0 . 76:44 

0.7803 

0.11468*01 

0.34388-02 

0.21648-02 

0.14968-02 

37 

C.78C8 

0.8070 

0.11528*01 

0. 34118-02 

C. 24998-02 

0.15918-02 

39 

0.9C62 

0.8325 

0.11318*01 

0.33178-02 

0.26558-02 

0.16948-02 

39 

0.8304 

0.8S69 

0.11078*01 

0.3216E-C2 

:.28! 3E-02 

0.17878-02 

40 

0.8533 

0.8799 

O.IC09E*O1 

o.3:39E-c: 

0.29:58-02 

0.18918-02 

41 

0.8748 

0.9016 

0.10458*01 

0.2989E-C2 

C. 3:49E-C2 

C.  :999E-C2 

42 

0.8949 

0. 9218 

0.9844E*0C 

0.28048-02 

0. 33518-02 

C.2119E-C2 

43 

C.9133 

0.9404 

0.9361E*00 

0.2655E-0: 

0.35518-02 

C.2239E-C2 

44 

0.93C1 

0,9573 

O.0514E*OO 

0,2408E-C2 

0. 376’8-02 

0.2366E-02 

45 

0,9452 

0.9726 

0.7083E*00 

0.2CC2E-02 

C. 40 '9E-C2 

0.2513E-02 

46 

0,9586 

0.9860 

0.45048*00 

0.1275E-02 

C.4490E-02 

C.26B8E-02 

47 

C.9701 

0.9977 

0.2796E*00 

0.7951E-C3 

C . 4  909E-02 

0.28778-02 

48 

C .9798 

1,0075' 

'0,n93E*01 

-0. 3410E-C2 

C,''C82E-02 

0.31538-02 

49 

0,9876 

1 .0154' 

-0,9638E*03 

-0.2784E*0: 

C.17!?e-C3' 

-0.44038*07 

50 

0.9936 

1.0213 

0.24088*12 

0.7C77E.6C9- 

•C'.5:44E-C8 

-0.3973e*19 

51 

C, 9975 

1.0253 

0.3259E*12 

0.980€E*C9- 

■C.2:88E*C9 

-C.4494E*20 

52 

0,9995 

1,0273 

0.9400E+12 

C.2907E‘1C- 

•C. 3:2iE*09 

-C.5061E-20 

•  •  • 

SUMMARY 

OF  BOUNDARY  LAYER  : 

SOLUTIONS  OF 

•  ISF  =  2 

NX 

XC 

S 

vw 

CF 

DLS 

THT 

1 

C.0024 

0,0000 

0.1232E*01 

0.0OC0E*0C 

O.00CCE*00 

0.0O00E*C0 

2 

0.0064 

0.0066 

0.8911E+00 

0.4236E-C1 

0.1316E-C3 

0.5759E-04 

3 

0,0123 

0.0144 

0.7976E*00 

0.2157E-01 

0,  l814E-:3 

0. 7B6'’E-04 

4 

0.0201 

0.0236 

C.7397E*00 

0.1440E-02 

0.2262F-03 

0. 9T51E-C4 
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5 

0.0298 

0.0344 

0.6680£«00 

O.lOJtE-Ol 

0.2790E-03 

0.11848-03 

t 

0.0413 

0.0468 

0.625U*00 

0. 7990E-02 

0.3328E-03 

0. 14018-03 

7 

0.054“^ 

0.0601 

0.5123£*00 

0.62S0E-02 

0.3921E-03 

0.16298-03 

s 

0.0(96 

0.0164 

0.5379E*00 

0.S2O2E-O2 

0.45b0E-03 

0.18148-03 

i 

0.0866 

0.0936 

0.512ie«00 

0.4426E-02 

0.5ie0E-03 

0.21128-03 

10 

0.1050 

0.2123 

0.4810C«00 

0.2B1SE-02 

0.5802E-03 

0.23588-03 

11 

0.1251 

0.1326 

0.4621E*0C 

0.3321E-02 

0.6472E-03 

0.26098-03 

12 

0.1466 

0.1S42 

0.44198*00 

0.2932E-02 

0.7162E-03 

0.28618-03 

13 

0.1695 

0.1112 

0.42838*00 

0.2644E-02 

0.7834E-03 

0.31218-03 

H 

0,1931 

0.2014 

0,40908*00 

0.23S4E-02 

0.B539E-03 

0.33188-03 

15 

C.2191 

0.2269 

0. 38428*00 

0.2092E-02 

0.933eE-03 

0.36538-03 

16 

0.2456 

0.2534 

0.31428*00 

0.1928E-02 

O.lOOeE-02 

0.39238-03 

n 

0.2131 

0.2809 

0.35828*00 

0.!753E-02 

0.1084E-02 

0.41918-03 

18 

C.3015 

0.3093 

C. 34118*00 

0.1593E-02 

0.n66E-02 

0.44698-03 

19 

0.3301 

0.3385 

0.32128*00 

0.1462E-02 

0.1247E-02 

0.41438-03 

20 

0.3605 

0.3683 

0.30328*00 

0.1301E-02 

0.1341E-02 

0.50328-03 

21 

0,3909 

0.3988 

0.28108*00 

0.1U2E-02 

0.1442E-02 

0.53338-03 

22 

C.4218 

0,4297 

0.26918*00 

0.1075E-02 

0.1535E-02 

0.56298-03 

23 

0.4530 

Q.46C9 

0.26438*00 

C.1021E-02 

0.1618E-02 

0.59098-03 

29 

0.4843 

0.4923 

0.24768*00 

0.9272E-03 

0.nilE-C2 

0.61878-03 

25 

0.5156 

0.5237 

0.21028*00 

0.7655E-03 

0.1844E-02 

0.64978-03 

26 

0.5469 

0.5551 

0.19838*00 

C.7C38E-03 

0.1962E-02 

0.68088-03 

2  ’ 

0.5181 

0.5863 

0.19978*00 

0.6912E-03 

0.2043E-C2 

C. 10828-03 

28 

C.6090 

0.6173 

0.16418*00 

C.5549E-03 

0.2182E-02 

0. 13158-03 

29 

0.6394 

C.6478 

C. 12198*00 

C.4237E-03 

0.2366E-02 

0.76978-03 

30 

0.6693 

0.6118 

0.11738*00 

0.3810E-03 

0.2504E-C2 

0. 19968-03 

31 

0.6984 

0.7071 

0,97148-01 

0.3098E-03 

0.2658E-02 

0. 82858-03 

32 

0.1268 

0,7356 

.,39928-01 

0.1220E-03 

0.2973E-02 

0.85928-03 

33 

0.7544 

0.1633 

0,16048-01 

C.2349E-03 

0.2607E-02 

0.76638-03 

34 

0.1808 

0.1 90C 

0.5406E-02 

0.1646E-04 

0.3036E-02 

0.79698-03 

35 

C.8062 

0.8155- 

0,82018-01- 

-0.2493E-03 

0.3390E-C2 

0.76398-03 

36 

0.9304 

0.8399- 

0.90418-01- 

■0.269ie-03 

C.33e4E-02 

0.79978-03 

37 

0.8533 

C.  9629- 

0.15628*00- 

-0.4602E-C3 

0.3787E-C2 

0. 53658-03 

38 

C.8148 

0.8846- 

C. 25118*00- 

■0.7339E-03 

C.446eE-02 

0.5^588-03 

39 

0.0949 

C.9C48- 

0.33608*00- 

-0.97;3E-03 

0.5324E-02 

0. 9352E-C3 

40 

0.9133 

0.9234- 

0.14388*01- 

-0.2’.50E-01 

0.2542E-C'. 

•0.5421800 

41 

0,9301 

0.9403- 

0.29898*07- 

-0.e597E.04 

0.7C34E-03 

•0.2021809 

42 

0.9452 

C.9556- 

0,32298*07- 

-0.927:e.04- 

-0,1513E.07 

43 

0.9586 

C.9691- 

0.19358*12- 

-0.5553E-09 

0.1354E-C9 

•:.2^468 -22 

44 

0.9101 

0,98C1- 

0.80498*14- 

-0.2316E.I2 

0.21181-10 

-0.55028*23 

45 

0,9198 

0.9905 

0.14128*17 

0.4082E‘14 

0.307CE-;4 

-0. ■’0588-30 

46 

0.9876 

0.9984 

NaN 

NaN 

NaN 

-NaN 

47 

0.9936 

1.0C44 

NaN 

NaN 

NaN 

-NaN 

48 

0.9975 

1,0083 

NaN 

NaN 

NaN 

-NaN 

49 

C.9995 

1,0104 

NaN 

NaN 

NaN 

-NaN 
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Gnuplot  command  file:  profile 


set 

terminal  tek40xx 

set 

nogr id 

set 

nclabel 

set 

si2e  1, 

1 

set 

data  St 

yle  lines 

set 

noxtics 

set 

ytics 

set 

title  - 

Veloc 1 

ty  Profiles  A0A=!0  deq.-  0,0 

set 

nokey 

set 

xlabel 

■Airfoil  Upper  Surface  Statloo"  0,0 

set 

xranqe 

10  :  101 

set 

y label 

■y/c" 

0,  .5 

set 

yrange 

(0  :  . 

010) 

set 

size  .€ 

,  .6 

set 

label  * 

Station  0  is“  at  .S,  .0C*8 

set 

label  * 

stagnation**  at  .8. -CO’’ 

set 

label  • 

point** 

at  ,8,  .0C6 

set 

label 

5”  at 

1,-.''CC3  center 

set 

label  * 

10-  at 

2,-.C0C3  center 

set 

label  •• 

15-  at 

3, -.0003  center 

set 

label  •’ 

20-  at 

4. -.0003  center 

set 

label  " 

25-  at 

5, - . CCC 3  center 

set 

label  •' 

30“  at 

. 0003  center 

set 

label  '* 

35-  at 

■»p«.:003  center 

act 

label  ’* 

40“  at 

8, -. :C03  center 

set 

la  be.  •• 

45-  at 

9, -.0003  center 

plot 

‘■pro!  0 

“lin 

1  o»* 
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Input  file;  FOR001.DAT  (Stardenl)  or  fort.l  (IRIS) 


3 


NACA  0 

•'!?  AIRFOIL 

ALP! 

PIVOT 

f.COOCOC 

0.250000 

SLOWER 

SUPPER 

50 

50 

C 

1.00000 

0.98000 

0.96000 

0.94000 

0.92CCC 

0.9CCC0 

0.68000 

0.86000 

0.84000 

0.82003 

0 .  eococ 

0. 780CC 

0.76000 

0.74000 

0.72000 

0.70000 

C. 68000 

0.66000 

0.64C00 

0.62000 

0.60000 

0.58000 

O.56C0C 

C.540CC 

0.52000 

0,50000 

0.48000 

0.46000 

0.44C0C 

0.42000 

0.40000 

0.38000 

0.36000 

0.340CC 

0.32CCC 

0 . 30000 

0.28000 

C.26C00 

0.24000 

O.22C0C 

C.2C0CC 

: .  isooc 

0.16000 

0.14000 

0. 1200C 

O.IOOOC 

0.C90C0 

0.C6OCO 

C.C4000 

C.C20C0 

0.00000 

0.02000 

C.C4COO 

0 .06000 

0.08000 

0. ICOOO 

C. 12000 

0.14000 

C.  16000 

0. 18000 

0.20000 

0.22COC 

0.24000 

0.26000 

C.28CCC 

C.  30000 

0.32000 

C. 34000 

0.36000 

0.36000 

0. 40000 

0.42000 

0.44000 

0.46000 

0.48000 

0.50000 

C.520CC 

0.540C0 

0.56000 

0.58COO 

0.60000 

0.62000 

0,64000 

0.66000 

0.60000 

0. 7CC00 

C. 72000 

C. 74000 

0.76000 

0 . 78000 

o.aoooo 

C.82C00 

C.840CC 

C.B60C0 

0. 89000 

0 .  oooco 

0.92C0C 

C . 94CCC 

C.9600C 

0.9800C 

1  .OOOOC 

“C.OC126 

-C.OC4C3 

-C.0C674 

-0.00939 

-O.C1196 

-0.01449 

-0.01694 

-n  n’  R  ■»- 

-O.C21"C 

-0 .02399 

-0.02623 

0 

80 

40 

-O.C3C56 

-0,33264 

-0.03467 

-O.C3664 

-0.03856 

04042 

-0. 04222 

-0.04396 

-0.04563 

-0.04723 

-0.0487S 

-0 . 0;026 

-0,05:65 

-C.C5294 

-0.05415 

-0.05530 

-0.05634 

-0.0:  2  fc 

rs  qn  1 

-0.05869 

-C. 05923 

-0.05966 

-C . 05995 

-0.C6CO6 

-0.0:997 

-0. 05966 

-0.05911 

-C. 05938 

-C.O573 • 

-0.05607 

-0.05444 

-C. 052  36 

-0.04990 

-0.04683 

-C. 04309 

-0,03912 

-0.0323: 

-C.C23S2 

O.CCOOC 

0.02392 

C.C3231 

0.03942 

0.04309 

r 

0.04990 

0.05236 

0.05444 

0.15607 

0.05737 

0.0:939 

0.05911 

0.C5966 

0.0599‘> 

O.C6C:6 

0.05995 

0.05966 

0.05923 

C  C5868 

0.05803 

0.:;'26 

0.05634 

0.C553O 

0.05415 

0.05294 

0.C5'.6: 

:. 0:026 

0.04978 

C.C4  ’?:• 

0.04563 

0.04396 

0.04222 

:.C<C42 

0.03856 

C.C3664 

0.0346' 

0.03264 

C.0^05<^ 

02842 

0.02623 

0.02399 

c.ori'c 

0.01935 

0.01694 

'.01448 

0.07 196 

0,00939 

0.006‘’4 

C.CC4C3 

c-.oo:2t 

Input  file:  incompbl.dat 


TWAKE  NXT 

1  161 

ITP(l)  ITR(2) 
C  C 

IP 
1 


NW  ITRENO 

37  40 

rSWPMX  ®L  XCTPtl) 

1  54C0CC,C  C.j:C?C 
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Output  file;  incompbl.out 


C  NACA  0012  AIRFOIL 

INPL’T  DATA  FOR  IMVISCID-FLOW  CALCULATIONS 

ALPI^  6.0000  PIVOT-  0.250C 

NL0WER=  50  NUPFER=  SC 


COORDINATES  OF  THE  BOOT 
X/C 


1  .cooooo 

0.980000 

0.960000 

0.94CC00 

0. 920CC0 

C. 900000 

0.880C0C 

C.86000C 

C.84CC0C 

C.820CCC 

C . 80C0C0 

0. 79C0CC 

0.760000 

C.740CC0 

0.720000 

C. 7CC000 

C.690CCC 

0.660000 

O.64OCO0 

0.620CCC 

0.6C0C0C 

0. 580000 

O.560COO 

0. 540OCO 

0.520000 

o.socco: 

C.49CCOO 

C-46CC00 

r-  ^  40^0^ 

0. 420CC0 

0. 4C0C0C 

C.38000C 

0. 3600CC 

c.34c:oo 

o!3:oocc 

C.300COC 

0.280CC0 

0.26CCCC 

C. 240000 

:'.22C0CC 

o.occoc: 

0.190000 

0.160000 

0. 14C0CC 

0.12OC0O 

0. lOOOCC 

C.C600CC 

c*.o6:o:c 

C.C4000C 

0.C20C0C 

c . cccooo 

0 . C2C000 

C.C400CC 

0.0600C0 

Q.ceocoo 

C. ICCCCC 

C. I2C00C 

ca4cccc 

C. 160000 

0.180000 

0,200000 

0.220000 

0.240000 

C. 260000 

C. 280000 

0. 3C000C 

0.320000 

C. 3400CC 

C.3600CC 

0.3900CC 

0.40CCCC 

0.420000 

0.440CO0 

C.460CCC 

C,49C0C0 

C.5000CC 

C . : rocoo 

0.5400C0 

0.56000C 

C.59CCC" 

C.600CCC 

0.620000 

0.64000C 

0.660COO 

Q.680CC0 

0. 70000:^ 

0 . *2C0CC 

C, 740000 

0. *60000 

f\  -  ^  A  ^ 

C.8000CC 

o.erocoo 

0.840000 

0.960000 

0.880000 

**  ^  ^  ^  r-  p 

0,82CC0C 

C  ^4000** 

/s  ^  'ifel'OCC' 

T/C 


-0.OC126O  -0.C0403C 

0.006740 

-C.0C938C  -:.ci:96c 

0:448C 

-0.016940  019350 

<*  22  *  ^  '^C 

-0.02399C  -j. 0262  30  -0 

028420 

-C.03C560  -C.C’?<40 

:.0346  7C 

-0.036640  -0.039560  -0 

:4C4:: 

-0.042220  -0.C4396C 

C.C45630 

-0.047230  -C. 048780  -0 

05C260 

-0.051650  -0.052940 

0.054150 

-0.055300  -0.056340  -0 

05-260 

-0.058030  -0.058680 

0.059230 

-C. 059660  -0.059950  -0 

060C6C 

-C.0599''0  -0.059660 

0.059110 

-0.058380  -0.057370  -? 

0560 '0 

-0.054440  -0.052360 

C.C499CC 

-0.046830  -0.043090  -0 

039420 

-0.032310  -0.023820 

c.ccoccc 

0.023800  0.032310  0 

0  3942'^ 

0.043090  0.046830 

0.0499:: 

C.C5236C  C. 05444:  C 

C560': 

0.057370  0.058380 

0.059110 

0.059660  C.0599'0  0 

06C06C 

0.059950  0.059660 

0.059230 

0.058680  C.C58C3C  0 

057260 

0.056340  0.05530C 

0.054150 

O.C52940  0.'05:6'-0  0 

:5C260 

0.048780  C. 047230 

0.04563C 

C. 043960  0.042220  0 

-4:42'' 

0.03856C  C.03664C 

0.034670 

0.032640  :.0.’0560  0 

•:2?420 

0.026230  0.02399C 

C.021'’C0 

0.019350  C.;'.(9'.;  0 

0 «  4  a 

0.011  960  0.009.180 

C.00674C 

C.004C30  0.00:260 

ISVISCID  FLOW  PF.SULTS 

PANEL  XP 

TP 

CP 

1  0.99000E-OC 

-C, :4950E'C2  C.25966£*0: 

2  0.97000E-00 

-C.4481 

3E-C2  C,38882E‘C0 

3  C.95COOE*OC 

-C.744: 

5E-02  C.147C3E*C0 

4  0.93000E*CC 

-0.1032 

:e-01  0.1I'35E‘0C 

5  o.9ioooe*oo 

-0.130615-01  0.98371E-0! 

6  o.e9oooe«oo 

-C.15646E-01  0.87319E-:i 

7  0.87000C-00 

-O.ieinE'OJ  0,  79369E-0I 

8  C.85000E«00 

-0.205C2E-01  0.72019E-C; 

9  C.8300CE*00 

-0.22829E-0!  0.65’20E-OI 

10  C.81000E*0C 

-C.25104E-0:  0.58929E-C1 

11  0.79COOE-00 

-0.27325E-01  0.53452E-0! 

12  0.77000E*OC 

-C.2  9490E-'01  C.  4eiC9E-Cl 

13  C.750C0e*C0 

-0.31600E-C1  0.43^94r-c: 
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H 

li 

16 

n 

18 

19 

20 
:i 
22 

23 

24 

25 
28 
21 
28 
29 
3C 

31 

32 

33 

34 

35 

36 
31 

38 

39 
4C 

41 

42 

43 

44 

45 

46 
41 

48 

49 
5C 

51 

52 

53 

54 

55 

56 
51 

58 

59 
62 
61 
62 

63 

64 

65 

66 
61 
68 
69 
1C 

12 

13 

14 
75 
16 
11 
18 
19 
8C 
81 
82 

83 

84 


0. IIOCCE'OO 
0. ’ICCCE-CO 
o.69Qc;e-co 
0.610CSE*00 
0.65CrCE-00 
D.63OC0E-OO 
0.610COE*00 
0.59000e*00 
0.510008*00 
0.55000E-00 
0.530008*00 
0.510008*00 
0.490008*00 
0.410008*00 
0.450008*00 
0. 4.10008*00 
0.410008*00 
0.390008*00 
0.310008*00 
0.350008*00 
0. 330008*00 
0.310008*00 
0.290008*00 
0.210008*00 

0.2300CS*0C 
C.21CC0E*0O 
0. ! 9CCC£*0C 

c.  nooc£‘Cc 

C  . '  5CC0J1*  00 
C.'  i  2C0C£*0C 

r>  •  ^  z'rs"***  *  Q  J 

oi^occvE-c: 
C, -COCOE-Cl 
C.:CCC:£*01 
C,3CC:0£-C1 

c.icecoe-ci 

C. • COCOE-01 
0. 3C000E-C’. 


-0. 33655r-03 
“0.356S5E-0) 
-0.33€OOE-C1 
“0. 39490E-01 
-0.41320E-01 
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0.23786 

0.231375 

0.00295 

0.0C0B4 

0.97198 

C.C5526 

3 

110 

C. 25939 

0.297918 

0.0028! 

0.00088 

0.97777 

0.04396 

3 

111 

0.27127 

0.269797 

0.00267 

0.00092 

0.98337 

0.03298 

3 

112 

0.28846 

0.281986 

0.00252 

0.00095 

0.98811 

0.02363 

3 

113 

0.30593 

0.299969 

0.00237 

0.00100 

0.99165 

C. 01662 

2 

114 

0.32367 

0.317203 

0.00222 

0.00109 

0.99917 

0.01163 

2 

115 

0.39169 

0.335171 

0.00210 

0.00109 

0.99595 

0.00808 

2 

116 

0.35980 

0.353339 

0.00200 

0.00113 

0.99732 

0.00535 

2 

117 

0.37819 

0.371683 

0.00193 

0.00117 

0.99871 

0.00258 

2 

118 

0.39667 

0.390221 

0.00187 

0.00121 

1.00031 

-0.00062 

2 

119 

0.91535 

0.908916 

0.00181 

0.00129 

1.00197 

-0.00395 

2 

120 

0.93414 

0.927728 

0.00173 

0.00128 

1.00311 

-0.00624 

2 

121 

0.95303 

0.496636 

0.00163 

0.00133 

1 .00390 

-C. 00682 

3 

122 

0.97198 

0.965613 

0.00159 

0.00138 

!  .00306 

-0.00612 

3 

126 


:2j 

C.  49095 

C. 484624 

0.00149 

12< 

C.5C994 

0.503652 

0.00145 

]25 

C.  52891 

0.522663 

0.00140 

126 

C.  54783 

0.541631 

0.00134 

127 

C.  56668 

0.560530 

0.00128 

128 

0.58543 

0.579334 

0.00123 

129 

0.60405 

0.598017 

0.00120 

130 

0.62252 

0.616553 

0.00116 

131 

0.64081 

0.634916 

0.00112 

132 

C. 65891 

0.653083 

0.00108 

133 

0.67677 

0.671027 

0.00104 

134 

0.69439 

0.688723 

0.00100 

135 

0.71173 

0.706149 

0.00096 

136 

0.72878 

0.723279 

0.00092 

137 

C. 74550 

0.740090 

0.00088 

138 

0.76188 

0.756558 

0.00084 

139 

0.77789 

0.772662 

0.00079 

1  40 

0.79352 

0.788379 

0.00074 

141 

0.80873 

0.803688 

0.00069 

142 

0.82352 

0.818567 

0.00064 

143 

0.83785 

0.832996 

0.00059 

144 

0.85172 

0.846954 

0.00053 

145 

C.865C9 

0.860424 

0.00047 

146 

0.97796 

0.873387 

0.00041 

147 

0.89C3C 

0.885823 

C. 00034 

148 

C. 90210 

0.897717 

0.00026 

149 

0.91334 

C.  9C9052 

0.00018 

150 

0.92400 

0.919813 

0.00009 

151 

0. 93407 

0.929985 

0.00002 

152 

0.94354 

0.939553 

-0.00006 

153 

0.9524C 

0.948504 

-0.00016 

154 

0.96063 

0.956827 

-C. 00030 

155 

0.96923 

C.96451C 

-O.CC050 

156 

0,97518 

0.971542 

157 

0.98149 

0.9''79!3 

-0'.C0!06 

158 

0.98'^13 

0.9B3615 

-0.00144 

159 

0.9921C 

0. 988640 

-C. 00181 

16C 

0.99639 

0.992982 

-0.00208 

161 

1 .00000 

C. 996632 

-0.00223 

wake 


T 

xs 

UE 

CP 

U1 

1.04119 

0,92565 

0.14318 

U2 

1.04769 

0,9U81 

0.15945 

U3 

3 .05492 

0.92128 

0.15124 

U4 

1.06297 

0,92771 

0.13935 

165 

1.07193 

0,9319C 

0.13157 

166 

1.08191 

0,93520 

0.1254! 

167 

1.093C3 

0.93674 

0.12251 

166 

1.10541 

0,93908 

C. 11813 

169 

1.11919 

0,94130 

O.II394 

170 

1.13454. 

0,94431 

0.10828 

171 

1.15163 

C.94‘’76 

0.10175 

172 

1 .17066 

0.95176 

0.09415 

173 

1.19185 

C. 95609 

0.08589 

174 

1.21545 

0,96029 

0.07784 

175 

1.24173 

C. 96418 

0.07035 

176 

1.27099 

0,96780 

0.06337 

177 

1.30357 

0,97110 

0.05696 

178 

1 .33985 

0,97422 

0.05090 

179 

1.38025 

0,97700 

0.04547 

180 

1 .42523 

0.97960 

0.04038 

181 

1.47533 

0.98193 

0.03582 

182 

1 .53111 

0.98409 

0.03157 

183 

1.59322 

0.98602 

0.02777 

184 

1.66239 

0.98778 

0.02428 

185 

1.73941 

0.98935 

0.02118 

0.00143 

1 .0C262 

-O.OC524 

3 

0.00147 

1  .OC243 

-0.00486 

2 

0.00150 

1.00230 

-0.00460 

2 

0.00155 

1.0C17B 

-0.00357 

3 

0.00160 

1.00087 

-O.OC173 

3 

0.00164 

0.99985 

0.00030 

3 

0.00168 

0.99895 

0.CC210 

2 

0.00172 

0.99806 

0.00387 

2 

0.00177 

0.997C2 

0.00595 

3 

0.00181 

0.99581 

0.00836 

3 

0.00186 

0.99450 

0.01096 

3 

0.00190 

0.99310 

0.01375 

3 

0.00195 

0.99162 

0.01670 

3 

0.00199 

0. 99006 

0.C1978 

3 

0.00204 

0.98842 

0.02303 

3 

0.00209 

0.98668 

C. 02647 

3 

0.00215 

0.98483 

0.03011 

3 

0.00221 

0.98287 

0.03397 

3 

0.00227 

0.96080 

0.03804 

3 

0.00233 

0.97962 

0.04231 

3 

0.00240 

0.97637 

0.0467C 

3 

0.00248 

C. 97407 

0.05119 

3 

0.00255 

0.97173 

0.35573 

3 

0.00264 

0.96935 

C.06C37 

3 

0.00274 

0.96693 

0.06305 

3 

0.00281 

0.96448 

0.0&9?e 

3 

0.00302 

0.96204 

0.07448 

3 

0.00320 

0.95970 

3 

0.00341 

0.95757 

0.36306 

4 

0.0C366 

0.95569 

0.39665 

a 

0.00391 

0.95402 

0.08995 

8 

0.00417 

C. 95231 

0. 09310 

8 

C. 00445 

0.95021 

0 . 0  9  ’  1 1 

2 

0.00469 

0.94749 

o!  1022? 

1 

0.00489 

0.94382 

o.:c92: 

9 

0.00501 

C. 93889 

C.  :  1949 

8 

0.C05C9 

0.93267 

0.13013 

0.30512 

0.92576 

o!i4296 

6 

0.00511 

0.91935 

0.1549': 

4 

OT.S 

II 

0.01781 

2.33C14 

c 

0.016C4 

2.00843 

2 

C. 01449 

1 .93997 

3 

0.01336 

l!^5i91 

3 

C.C1264 

1.68424 

3 

0.01209 

1,63209 

3 

o,c:i69 

1. 598^6 

2 

0,01129 

1 .5^ '34 

2 

0.01C92 

1.50866 

2 

0.C1051 

1,47057 

2 

0.01012 

1.43422 

z 

0.00974 

1.39982 

2 

0.00937 

1.36790 

i 

O.OC9C3 

1.33017 

7 

0.00873 

1.31099 

2 

0.00846 

1.28600 

2 

0,00822 

1.26354 

2 

0.00800 

1,24320 

2 

0,00781 

1.22443 

2 

0,00763 

1.2C7C6 

2 

0,00748 

1.19108 

c 

0.00733 

1.17633 

2 

0,00720 

1.16290 

2 

0,00709 

1,15040 

2 

0.00699 

1.13904 

2 
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18€ 

1.82517 

0.99C78 

0.01835 

O.C0689 

1 .12864 

2 

1 .92067 

0. 992C4 

0.01585 

0.00681 

!. 11910 

2 

188 

2.02702 

0.99317 

0.01361 

C. 00673 

1.11030 

2 

189 

2.14543 

0.99416 

0. 01164 

C. 00666 

1  .  1C21  7 

2 

190 

2.27730 

0.99504 

0.C0990 

C .  C  C  6  6  0 

1 .09469 

2 

191 

2,42413 

0.99580 

C.0C839 

0.CC654 

1 .08789 

2 

192 

2.58763 

0.99647 

O.OC7C5 

C.CC649 

1 .081 71 

2 

193 

2,76969 

0,99705 

0.00590 

C.0C645 

1 .07607 

2 

19<1 

2.97243 

0.99755 

0.00490 

C.0C641 

1 .07090 

2 

195 

3.19818 

0.99797 

0.00405 

0.00639 

1  .06616 

2 

196 

3,44956 

0.99833 

O.OC333 

C.C0635 

1.0617" 

2 

197 

3,72949 

0.99863 

0.00273 

0.00632 

1.05778 

2 

198 

4.04119 

0.99859 

0.00281 

O.OCOCO 

O.OCCOO 

SUMMARY  OF  THE  DRAG,  LIFT  AND  PITCHING  MOMENT 
COEFFICIENTS  WITH  THE  CYCLE 


CD  IS  EVALUATED  FROM  FAR -WAKE  FORMULA 
Cl.  (  CM  FROM  INTEGRATION  OF  CP 


c*  r 

CD 

CL 

CM 

1 

C. 012592 

C. 011899 

0.665033 

2 

0.011748 

C. 012118 

0,6655 w 

3 

C.01C915 

C. 011455 

0.6659C6 

4 

0.012261 

0.011899 

C. 665108 

0.012459 

0.012117 

0.666036 

6 

0.0H949 

0.011452 

0.666010 

1 

C.012C65 

0.C11898 

0.665143 

8 

0.012272 

0.012116 

C. 666060 

9 

C.Ol 1606 

0.011459 

0.6660  ’f 

1 0 

0.011963 

0.011896 

0.6651:: 

1 1 

0.012167 

0.659153 

0.0C2397 

12 

0.011474 

0.6586C6 

0.002286 

13 

C. 011860 

C. 658523 

0.002802 

14 

C.012C45 

0.658816 

C.CC251 7 

15 

0.011453 

0,657755 

0.0C250C 

16 

C. 011853 

0.659142 

0.002906 

1  " 

0.012176 

0.658669 

0.002550 

19 

0.01160C 

0.659903 

0.'C0204: 

19 

0.011900 

0.661304 

C.0C22“3 

2C 

0.012125 

0.662937 

0.0C1650 

21 

0.011497 

0.663051 

0.001391 

22 

0.011987 

0.663583 

0.001815 

23 

0.012100 

0.664779 

o'oc:284 

24 

C.Cr.4  66 

0.664628 

D 

2  5 

C. 011872 

0.664564 

c.:o:6i  3 

2  6 

C. 012094 

0.665672 

O.OCH24 

2  7 

0.011457 

0.665504 

0 . 00091 0 

28 

0.011898 

0.664903 

0.0C152C 

29 

0.012119 

0.665869 

0.001081 

3C 

0,011455 

0.665745 

0.000863 

3; 

0.01  1899 

0.665033 

0.CC1489 

32 

0.012118 

0,665977 

0.001061 

33 

C. Cl  1455 

0.665906 

C.00C835 

3^ 

0.01 1898 

0.6651C8 

0.001471 

35 

C.012m 

0.666036 

C.CCI05C 

36 

C. 011452 

C. 666010 

C.0CC917 

37 

0.011898 

0.665143 

0,001461 

38 

0.012116 

0 . 666C6C 

0.001C46 

39 

0.C11458 

0. 666076 

C.CCC8C6 

4C 

0.011896 

0.665151 

0.00145? 
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Source  code:  pres.f 


2 

3 

4 

5 

6 
? 

8 

9 

1C 

n 

12 

13 

14 

15 
U 

n 

18 

19 

2C 

21 

22 

?3 

24 

25 
2f 

27 

28 
29 
3C 


program  pres 
REAL  cp(lOO) 

OPEN<unlc=20, fi le=*cp.ir* ,sta-us=' unknown' > 

print  ‘j'no.  points  before  bubble?* 
read  ‘.N 

print  ‘j'no.  points  after  bubble?' 
read  ’jM 

C  cp(l)  corresponds  to  stagnation  for  M=.3,  dark  fringe  5.5 

C  cp<21  corresponds  to  next  dark  fringe  4.5,  etc. 

do  1=1, N 

read (20, *  I  cp ( 1) 
end  do 


dc  1=1. N 

read ( 1, *  J  k, x, y 
wrlte{2,*l  X, -cp  < 1 ) 
end  do 


dc  i=N,N-M*l,-l 
read  (1 , *)  k, x, y 
write  (2.  •)  X, -cc  f i 1 
end  do 

end 


Page  1 
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APPENDIX  B:  BASIC  COMPUTER  COMMANDS 


The  following  collection  of  commands  is  intended  as  a  very  basic  user’s  guide  to 
the  various  support  programs  which  are  necessary  to  do  research  on  a  UNIX  system. 
It  is  by  no  means  a  comprehensive  list,  simply  enough  to  get  started  without  wasting  a 
great  deal  of  time  on  finding  elementary  procedures  and  syntax.  In  some  commands,  the 
arbitrary  word  filename  or  abbreviation  fii  is  used.  In  others,  specific  examples  are  used 
when  it  aids  the  clarity  of  the  explanation.  Both  methods,  however,  indicate  that  the  user 
may  substitute  an  appropriate  name. 


130 


BASIC  UNIX  COMMANDS 


man  filename 

obtain  on-line  help  information  for  a  program  or  command 
from  the  on-line  manual 

q)  fhl  fh2 

copy  filename  1  to  filename2,  both  in  current  directory 

q)  /alpha/ nowak/bl/fnl  . 

copy  fnl  from  another  directory,  path  specified,  to  the  as 
name  in  the  current  directory 

q)  ../fhl  . 

copy  fnl  from  the  directory  above  to  the  same  name 

mv  fhl  fn2 

move,  or  rename,  fnl  to  fn2  ;  fnl  will  no  longer  exist 

more  fn 

type  the  text  of  the  file  on  the  screen,  read  only 

!v 

repeat  the  last  command  that  started  with  v 

cd  paneldir 

change  directory  to  paneldir 

rm  fn 

delete  (remove)  fn 

mkdir  paneldir 
rmdir  paneldir 

make  directory  (example  name  paneldir) 
remove  (delete)  directory 

Is 

list  contents  of  directory  (like  dir  on  a  pc) 

pwd 

print  working  directory 

batch  <  fn 

execute  a  command  file  in  batch  mode  (runs  even  after 
logging  ofO  -  useful  for  long  run-time  programs 

telnet  131.120.254.92  (suzqt -Stardent) 

131.120.254.91  (madmax  -  IRIS) 


BASIC  EDITING  COMMANDS  FOR  THE  VI  EDITOR 


NOTE:  ALL  commands  in  VI  are  case  sensitive,  type  exactly  as  shown.  Check  the 
status  of  the  CAPS  LOCK  key  if  a  command  does  not  seem  to  work  properly. 


vi  fill  (invokes  editor,  calls  old  file  if  it  exists,  otherwise  creates  new 

file) 

To  get  started  into  text  mode  from  command  mode: 

a  (add  to  document,  cursor  moves  to  right  and  input  is  enabled) 

i  (insert,  cursor  does  not  move  and  input  is  enabled) 

o  (open  a  new  line  below  the  cursor,  input  is  enabled) 

To  return  to  command  mode: 

Esc  (disables  input,  enables  move,  write,  save,  etc.) 


NOTE:  All  of  the  following  commands  assume  command  mode 
To  move  around: 


ctrl-h 

move  1  left 

ctrl-1 

move  1  right 

ctrl-k 

move  1  up 

Ctrl-1 

move  1  down 

The  above  commands  will  always  work.  On  some  machines,  the  arrow  keys  MAY  work 
as  well.  Other  commands: 


l(shift)G 

(shift)G 

/bl2D 

/(Enter  key) 


go  to  first  line 
go  to  last  line 

search  for  the  next  text  string  "bl2D"  after  the  cursor  location, 
CASE  SENSITIVE 

search  for  another  occurrence  of  the  previous  search  string 
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When  done  editing  or  viewing: 

(Esc):w  write,  or  save,  but  do  not  exit 

(Esc):q  quit,  exits  only  if  no  modifications  were  made 

(Esc);q!  quit,  exits  without  saving  even  if  modifications  were  made 

(£sc);wq  write  quit,  exits  and  saves  all  modifications  to  the  original  fn 

(Esc):wq  newfn  write  quit,  exits  and  saves  all  modifications  to  a  new  fn 

To  delete: 

dd  delete  current  line 

2dd  delete  current  line  and  next  line 

lOdd  delete  current  line  and  next  9  lines 

To  cut  (copy)  and  paste: 

5yy  yanks  5  lines  to  buffer  (leaves  original  5  lines  also) 

move  cursor  to  desired  location 
p  pastes  the  5  lines 
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BASIC  GNUPLOT  COMMANDS 


NOTE:  The  following  is  a  list  of  some  example  commands.  Extensive  on-line  help  is 
available  by  typing  help  while  in  gnuplot. 

set  term  tek40xx  (Stardent) 

plot  "VEL.DAT" 

plot  "VEL.DAT"  with  Unes 

set  xrange  [0: 1] 

set  yrange  [0:5] 

plot  "cf.dat",  "dls.dat"  (two  different  data  files) 

plot  "cf.dat"  with  lines,  "dls.dat"  with  lines 

plot  "bl.dat"  using  2:5  (one  data  file  with  multiple  columns) 

plot  "bl.dat"  using  2:5,  "bl.dat"  using  2:6 

set  data  style  lines  (option  with  lines  will  then  not  be  needed  after  each  plot  command) 

set  key  (legend) 

set  nokey  (no  legend) 

set  grid 

set  nogrid 

set  nozeroaxis 

To  print  (these  commands  can  be  put  in  a  command  file): 
set  term  postscript 
set  output  "gnuout" 
replot 

set  term  tek40xx  (reset  to  terminal  being  used) 

set  output 


134 


FILE  TRANSFER  USING  FTP 


ftp  suzqt 
madmax 

get  fn 

get  oldfn  newfn 

put  fn 

quit 


(numerical  computer  address  may  be  used  instead) 
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INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  Virginia  22304-6145 

2.  Library,  Code  0142  2 

Naval  Postgraduate  School 

Monterey,  California  93943-5002 

3.  Dr,  M.  F,  Platzer  5 

Dept,  of  Aeronautics  and  Astronautics,  Code  AA/PL 

Naval  Postgraduate  School 
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4.  Professor  M.  S.  Chandrasekhara  1 

M.  S.  260-1 

NASA  Ames  Research  Center 
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5.  Dr.  L.  W.  Carr  1 

M.  S.  260-1 

NASA  Ames  Research  Center 
Moffett  Field,  California  94035 
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Chief,  Fluid  Dynamics  Division 

U.  S.  Army  Research  Office 
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